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Zusammenfassung 



Quantenchromodynamik [QCD] wird heute als die Theorie betrachtet, welche die Starke 
Wechselwirkung zwischen den fundamentalen Konstituenten der Hadronen, den Quarks 
und Gluonen, korrekt beschreibt [1-5]. Bei grof^en Energien (kleinen Abstanden) ver- 
schwindet die Kopplungskonstante - ein Phanomen, das gemeinhin als asymptotische 
Freiheit bezeichnet wird. In dem entsprechenden Impulsbereich wurde die Storungstheorie 
angewendet und in tiefinelastischen Streuexperimenten erfolgreich iiberpriift. Fiir mittlere 
und kleine Impulse hingegen wird die Kopplungskonstante so groi3, dass die Storungsthe- 
orie nicht mehr angewendet werden kann. Daher miissen andere Methoden angewendet 
werden, um das Confinement zu erklaren, also das Phanomen, dass im Experiment nur 
farblose Hadronenzustande beobachtet werden. In dieser Dissertation wird der Quark- 
Sektor der QCD in Coulomb-Eichung mit Hilfe der Dyson-Schwinger-Gleichungen unter- 
sucht. In diesem Rahmen nutzen wir verschiedene Strategien, um die Eigenschaften der 
QCD sowohl bei grofien als auch bei kleinen Impulsen zu studieren. 

Im ersten Kapitel erinnern wir an einige grundsatzliche Eigenschaften der QCD. Aus- 
gehend von der Lagrangefunktion der QCD prasentieren wir die Eichfixierung und mo- 
tivieren unsere Wahl der Coulomb-Eichung. Verschiedene Aspekte des Confinements wer- 
den diskutiert, insbesondere der Gribov-Zwanziger-Confinement-Mechanismus und seine 
Relevanz in Coulomb-Eichung. 

Das zweite Kapitel beschaftigt sich mit der grundsatzlichen Ableitung der Dyson- 
Schwingen-Gleichungen. Funktionale Methoden werden eingefiihrt und der Quark-Pro- 
pagator sowie die Quark-Beitrage zur Gluon-Zweipunktfunktion und zur Quark-Gluon- 
Vertexfunktion werden formal abgeleitet. 

Im dritten Kapitel werden diese Funktionen in 1-Loop-Storungstheorie ausgearbeitet. 
Um die in den Gleichungen auftretenden nicht-kovarianten Loop-Integrale in Coulomb- 
Eichung zu behandeln, wird eine neue Methode basierend auf Differentialgleichungen und 
partieller Integration entwickelt. Physikalische Resultate werden verifiziert, so z.B. die 
Gltigkeit der analytischen Fortsetzung zwischen Minkowski- und Euklidischer Raum-Zeit 
und die Renormierung der Quarkmasse. Des weiteren wird der Quark Beitrag zum 1- 
Loop-Koeffizient der /3-Funktion berechnet. 

Das vierte Kapitel ist der Slavnov-Taylor-Identitat der Quark-Gluon-Vertexfunktion 
gewidmet. Insbesondere wird das Auftreten des sogenannten Quark-Geist-Streukerns un- 
tersucht. 

In Kapitel 5 nahern wir uns dem Confinement-Problem durch Betrachtung des Gren- 
zfalls schwerer Quarks. In diesem Limes nutzen wir den (vollstandig nicht-storungsthe- 
oretischen) funktionalen Formalismus kombiniert mit einer Entwicklung nach Potenzen 
des Inversen der schweren Quarkmasse. Durch Einschrankung auf die fiihrende Ordnung 
in dieser Massen-Entwicklung leiten wir eine strenge analytische Losung fiir den Propag- 
ator der schweren Quarks ab. AnschlieBend nutzen wir die Gleichungen fiir gebundene 
Zustande von Mesonen und Baryonen, um das linear wachsende Potential abzuleiten, das 
Quark-Confinement erklart. 



Kapitel 6 behandelt die Vier-Punkt Greenschen Funktionen der Theorie. Ausgehend 
vom in Kapitel 2 eingefiihrten Funktionalformalismus werden diese Funktionen explizit 
abgeleitet und ihre Beziehung zu den Gleichungen gebundener Zustiinde aus dem vor- 
angegangenen Kapitel diskutiert. 

Kapitel 7 enthalt die Zusammenfassung und das Fazit. Es folgen die Anhange in denen 
unter anderem die in Kapitel 3 abgeleiteten nicht-kovarianten Integrale explizit iiberpruft 
werden und auch einige in der Arbeit benotigte Zwei- und Drei-Punkt-Integrale berechnet 
sind. 



Abstract 



Quantum Chromodynamics [QCD] is widely believed to be the correct theory of strong 
interactions between the fundamental constituents of the hadrons, the quarks and gluons 
[1-5]. At high energies (small distances), the coupling between quarks and gluons tends 
to zero, a phenomenon known as asymptotic freedom. In this momentum region, per- 
turbation theory has been applied and successfully tested in deep inelastic processes. At 
intermediate and low momenta though, the coupling constant becomes strong enough 
to invalidate perturbation theory. Different methods must be employed to investigate 
color confinement, i.e. the phenomenon that only colorless hadronic states are observed 
in the experiment. This thesis deals with the quark sector of Coulomb gauge QCD and 
the method we employ is the Dyson-Schwinger equations. In this framework, we utilize 
different strategies to explore both the large and small momentum properties of QCD. 

In the first chapter we review some basic properties of strong QCD. Starting with 
the QCD Lagrangian, the gauge fixing procedure is presented, and our choice of using 
Coulomb gauge is motivated. Certain aspects of confinement are discussed, in particular 
the Gribov-Zwanziger confinement mechanism and its relevance in Coulomb gauge. 

The second chapter is concerned with the formal derivation of the Dyson-Schwinger 
equations. Functional methods are introduced and the quark propagator, along with the 
quark contribution to the gluon two-point functions and the quark-gluon vertex functions 
are formally derived. 

In the third chapter, these functions are examined at one-loop perturbative level. To 
handle the Coulomb gauge noncovariant loop integrals entering the equations, a new 
method based on differential equations and integration by parts technique is developed. 
Physical results are verified, such as the validity of the analytic continuation between 
Minkowski and Euclidian space, and the quark mass renormalization. The quark contri- 
bution to one-loop coefficient of the /3- function is also calculated. 

Chapter 4 is devoted to the Slavnov-Taylor identity for the quark-gluon vertices. In 
particular, the appearance of the so-called quark-ghost scattering kernels is explored. 

In Chapter 5, we address the problem of confinement by restricting ourselves to the 
heavy quark sector of the theory. In this limit, we employ the (full nonperturbative) 
functional formalism, combined with an expansion in the inverse of the heavy quark mass. 
Restricting to the leading order in the mass expansion, we derive an exact, analytical 
solution for the heavy quark propagator. We then consider the bound state equations for 
mesons and baryons and use them to derive the linearly rising potential which confines 
the quarks. 

Chapter 6 is devoted to the four-point Greens functions of the theory. Based on the 
functional formalism introduced in Chapter 2, these functions are explicitly derived and 
their connection to the bound state equations from the previous chapter is discussed. 

Chapter 7 includes the summary and conclusions. It is followed by appendices where, 
among others, the noncovariant integrals derived in Chapter 3 are checked explicitly, and 
also some standard two- and three-point integrals needed in this work are evaluated. 



Chapter 1. 

QCD in Coulomb gauge 



A deep understanding of QCD requires a whole toolbox of theoretical methods: analytical 
perturbative methods for weak coupling, numerical lattice gauge theory, along with the 
canonical approach [6], and functional methods [2]. The latter two methods have the ad- 
vantage that they are not restricted to weak coupling and can still be treated analytically. 
The functional methods most commonly used are the renormalization group equations 
(see, for example, Refs. [7, 8] for reviews) and the Dyson-Schwinger equations of motion 
for the Green's functions of the theory [9-11]. In this thesis, we employ the functional 
equation techniques and derive the Dyson-Schwinger equations. 

Since QCD is a non-abelian gauge theory, within the functional approach considered 
here it is necessary to fix the gauge. Thus, after briefly introducing the QCD Lagrangian, 
we will present the gauge fixing procedure and the problems related to it, with emphasis 
on Coulomb gauge. We will then motivate why among various gauges. Coulomb gauge has 
the advantage that it is "physical" : after converting to first order formalism, we will show 
that the number of dynamical variables reduces to the number of the physical degrees of 
freedom. Further, the Gribov-Zwanziger confinement scenario [12, 13] will be introduced 
and its relevance for Coulomb gauge will be put forward. The alternative approaches 
to QCD in Coulomb gauge will be also reviewed, and in particular, the results obtained 
within the Hamiltonian formalism will be outlined. 

Although Coulomb gauge seems to be more efficient in identifying the physical degrees 
of freedom, noncovariance introduces severe technical difficulties^ and moreover, the prob- 
lems related to renormalization have not yet been solved. From a practical point of view. 
Landau gauge has the advantage of being covariant and thus many infrared investigations 
have been undertaken, however studies are still in progress. A brief review of the results 
obtained in this gauge, also in correspondence to the confinement mechanism, will be 
presented at the end of this chapter. 

1.1. QCD as non-abelian gauge theory 

QCD is a non-abelian gauge theory whose matter constituents, the quarks, are spin 1/2 
fermions and obey the Dirac Lagrangian^ 

Cg = qo^ix) [i7^D^ - m]^^ qfsix), q = r^t/ (1.1) 

where the Dirac 7^ matrices satisfy the Clifford algebra, {7^,7^^} = 2^^^^ and the indices 
a,P . . . commonly denote the Dirac spinor, flavor and (fundamental) color. The quark 



'^The noncovariant Feynman integrals in Coulomb gauge will be examined in Chapter 3. 
■^Initially we use Minkowski metric defined in the Appendix A. 
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fields q, q transform in the fundamental representation of the gauge group SU{Nc)^, with 
Nc = 3 realized in QCD. The covariant derivative (in the fundamental representation) is 
given by 

D^ = d^- igAf,, (1.2) 

where g is the coupling constant of the theory and Afj^{x) = ^^(x)T". The non-abelian 
gauge field ^^(x) transforms according to the adjoint representation of the gauge group. 
Given an infinitesimal transform U{x) = l — i6°'(x)T°' the variation of the gauge and quark 
fields is 

M«(x) = -i^-(x)r(x) (1.3) 
6qa{x) = -ir"r(rE)g„(x) (1.4) 
where the covariant derivative in the adjoint representation reads 

^- = 5-a^+5r^M^. (1.5) 
Defining the field strength tensor 

F^u = T'^id^At - d,A1 + gr^'^A^Al) (1.6) 
we can construct a kinetic term for the non-abelian gauge field A^ 

2 ^ 4 

such that the total Lagrange density Cqcd = Cg + Cym is invariant under local gauge 
transformations. 

Let us now consider the functional integral 



Cym = -^TtF^uF'''' = --F^,F''''\ (1.7) 



Z = j V[Aqq]exp{iSQCD}, (li 



where denotes the functional integration measure for the Yang-Mills and quark 

fields and the QCD action is given by 

Sqcd = I d'x |g-„(x) lij^Df^ - m]„^ q^{x) - ^F^.F'^'^^j . (1.9) 

The difficulty with the functional integral Eq. (1.8) is that the measure runs over 

infinitely many gauge equivalent configurations (field configurations that are connected 
by gauge transformations), whereas the the action Sqcd is gauge invariant. Hence it 
introduces for every gauge orbit^ a divergent factor (the volume of the gauge group). 
The way to handle this problem is to use a method introduced by Faddeev and Popov 
[14]. The idea is to single out one representative from each orbit by imposing a gauge 
fixing condition to the functional integral Eq. (1.8). In this thesis we shall be considering 
only Coulomb gauge, but also other choices such as Landau gauge are possible. As a 
consequence, a new set of Grassmann fields, known as ghosts, are introduced. The gauge 
fixing procedure, the ghost fields and the problems associated with them will be discussed 
in the next section. 



■^Some relevant formulas for the group SU{Nc) are collected in the Appendix A. 

*A gauge orbit contains all configurations connected by gauge transformations and will be explicitly 
defined in the next section. 
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1.2 Gauge fixing and ghosts 




gauge orbit 



Figure 1.1.: Depiction of a gauge orbit containing and the gauge transformed field 

1.2. Gauge fixing and ghosts 

We start by defining the gauge orbit for some configuration as the set of all gauge 
equivalent configurations, i.e. each point on the gauge orbit is obtained by acting 
upon Afj, with the gauge transformation U (see also Fig. 1.1): 

[^r**] =: {< = UA^U^ - ^{d^U)U^ I U e SU{N,)}. (1.10) 

Then the integral over all the gauge fields J T>A can be written as an integral over a full 
set of gauge-inequivalent configurations / 'DA'^'^f , where A'''^-^' is a reference gauge field 
representative for the orbit (i.e., integral over all possible gauge orbits), and an integral 
around each gauge orbit f VO{A^^f) 

JvA = J VA'-^f j VOiA^'^f). (1.11) 

Since one integrates over infinitely many equivalent configurations related by a gauge 
transformation, the integral around the gauge orbit J T>0{A^^f) is infinite and must be 
eliminated. The strategy is to extract only one representative gauge-field configuration out 
of each orbit by imposing a gauge fixing condition x[^] = 0- Each gauge field configuration 
on the orbit 0{A^'^^ ) is a gauge transformation of A^^f , i.e. one can (up to a volume factor) 
rewrite the integral J T>0{A^'^^) as an integral over the gauge group J T>U . The gauge 
fixing condition is implemented by inserting the identity 

1 = A[^] j Pi7<5[x(A^)] (1.12) 

into the functional integral Eq. (1.8). Ideally, the gauge fixing condition x(^) = should 
be satisfied by only one A^^ of each gauge orbit. ^ The Faddeev-Popov determinant 
A [A] = A[74^] accounts for the functional determinant arising from the argument of the 
delta function. Using the invariance of the action under gauge transformations, one can 
rewrite the functional integral of the theory as 



Z= j VUZgf (1.13) 

where the gauge fixed amplitude is given by 

Zgf = I V[Aqq]A[A]5[xiA)]exp{iSQCD}, (1.14) 



'But, as will shortly be discussed, this is in practice impossible due to topological restrictions. 
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and the divergent measure J DU has been factorized as an overall constant which can be 
absorbed in the normalization. In the following, we will use the notation Z (instead of 
Zgf) for the gauge fixed functional integral, since no confusion can arise. 
In Coulomb gauge, the (noncovariant) gauge fixing condition is given by: 

x[A]=V-A- = (1.15) 

and hence the temporal and spatial components of the gauge field must be treated dif- 
ferently. The gauge fixing condition can be implemented by rewriting the delta function 
with the help of a Lagrange multiplier A'* 



6[x{A)] ^ J PAexp 

and the Faddeev-Popov determinant 

A[A] = Det [V-D""^] 



(1.16) 



(1.17) 



can be written as a functional integral over two new Grassmann valued fields c and c, the 
so-called "Faddeev-Popov ghosts"^: 



Vice] exp 



i / d'^x7fV-D'"'d 



ab „b 



(1.18) 



In the above, we have introduced the spatial component of the covariant derivative 
Eq. (1.5) (in the adjoint representation) 



(1.19) 



Due to their unusual spin-statistics (they obey Fermi-Dirac statistics, and at the same 
time are scalar), the ghost fields are only allowed to appear in closed loops in Feynman 
diagrams and never as initial or final states in a physical process. Putting all these 
together, we can write for the Coulomb gauge functional integral: 



Z = i V[Aqqcc\] exp {\Sqcd + i^Fp}, 



with 



Sqcd 
Sfp 



d^x 



i7°-Do + i7 • 



m 



al3 



4 1^^ 



(1.20) 

(1.21) 
(1.22) 



In the QCD action, Eq. (1.21), we have separated the covariant derivative into temporal 
and spatial components (implicitly in the fundamental color representation), given by: 



a° - i^TV", D = V + igT^A" 



(1.23) 



where the temporal component of the gauge field A^'^ has been renamed to a". Spp 
collects the terms originating from the gauge fixing condition Eq. (1-16) and the Faddev- 
Popov determinant Eq. (1.18). 



pedagogic introduction on this topic can be found in [4]. 
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r 



[A] 




Figure 1.2.: Illustration of the hypcrplane F in gauge field configuration space obtained by gauge 
fixing. Furthermore, the first Gribov region D, and the fundamental modular region A 
are shown. A gauge orbit [A] intersects the hypcrplane F several times thus generating 
Gribov copies. The fundamental modular region, by definition, is intersected only 
once. 

The gauge fixing procedure described above is not yet complete, in the sense that 
the simple Faddeev-Popov trick is not sufficient to extract a single gauge configuration 
from each gauge orbit. The reason is the presence of the so-called Gribov copies [13], i.e. 
configurations Afj_ connected by a gauge transformation that produce multiple intersection 
points of a gauge orbit with the hypcrplane T generated by the gauge fixing condition (see 
also Fig. 1.2). In order to avoid the Gribov copies, it is necessary to restrict the hyperplane 
F to the so-called Gribov region Q,. This is obtained by minimizing the L^-norm of the 
vector potential 

Fa[U] = I d^xtT[AY{xf (1.24) 

along the gauge orbit [15]. In this region, any local minimum of this norm implements the 
gauge fixing condition (in this case. Coulomb gauge) and the Faddeev-Popov operator is 
restricted to positive eigenvalues: 

= jl: V-i'= 0;-V-l5 > o| . (1.25) 

Importantly, the Gribov region J7 contains the trivial configuration gA = (i.e., it contains 
all configurations relevant for perturbation theory) and therefore the ultraviolet regime is 
not influenced by this restriction. Moreover, the Faddeev-Popov determinant vanishes on 
the boundary dfl of the first Gribov horizon [16]. 

In general, is still not free of Gribov copies^, thus in principle one has to restrict 
further the gauge field configurations to the fundamental modular region A - the region 
of global minima of the norm defined above: 

A= [A:Fa[1]<Fa[U] yu}. (1.26) 

This gauge condition is also known as "minimal Coulomb gauge" [12]. In practise, these 
configurations are extremely complicated to identify®. However, in the continuum Zwan- 
ziger showed by means of stochastic quantization that Gribov copies inside the Gribov 

^A proof of the existence of Gribov copies inside f2 was given in [17]. 
*For recent lattice studies see, for example, [18-21]. 
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region do not affect the Green's functions of the theory [16]. Precisely, the dominant 
configurations he on the common boundary of Q and A and hence, in practise, restriction 
to the Gribov region is sufficient. This restriction might eventually generate nontrivial 
boundary terms which could influence the derivation of the field equations of motion (and 
implicitly the Dyson-Schwinger equations), but since by definition the Faddeev-Popov 
operator vanishes on the boundary of fi, these boundary terms are identically zero. 

1.3. Approaches to Coulomb gauge QCD 
1.3.1. First order functional formalism 

As already mentioned, the gauge of choice in this work is Coulomb gauge. Although this 
gauge does not have so many practical advantages as Landau gauge for example, there 
are several reasons motivating our choice: in Coulomb gauge there is a natural picture 
of confinement, Gauss law is naturally built in (such that in principle gauge invariance 
is fully accounted for) and the total color charge is conserved and vanishing [12, 22]. 
However, except for the original work of Khriplovich [23], only recently Dyson-Schwinger 
studies in Coulomb gauge have been undertaken. The technical barrier stems from the 
so-called energy divergence problem - the unregulated divergences generated by the ghost 
loops [24-26]. The way to circumvent this problem is to use the first order functional 
formalism^. Within this formalism these divergences cancel exactly and moreover, the 
system reduces automatically to "would-be-physical" degrees of freedom [12]. All these 
aspects will be discussed in detail in the course of this section. 

Motivation and main idea 

The main advantage of working in Coulomb gauge stems from the fact that in this gauge 
the Gribov- Zwanziger mechanism of confinement^" becomes particularly transparent. In 
this picture, the long range confining force is provided by the instantaneous Coulomb 
interaction, which appears to be enhanced for small three-momenta — )■ 0, whereas the 
transverse (colored) gluon is suppressed, reflecting the absence of the colored states in 
the spectrum. In covariant gauges such as Landau gauge, a quantity that leads to the 
confinement potential has not been identified so far. 

The conversion to the first order (or phase-space) formalism is motivated by the fact 
that within this formalism the famous Coulomb gauge energy divergence [24-26] can 
be avoided. Energy divergence means that the functional integral Eq. (1.20) gives an 
ill-defined integration, stemming from the energy-independent ghost loops^^ which are 
integrated over both 3-momentum and energy. As an example, consider the following 
one- loop integral [30]: 

J dko j d^k[{k - p)^P]~\ (1.27) 

^However, since we are mainly concerned with the quark sector of QCD, in the course of our investigations 

we will come across few points where restriction to second order formalism [27, 28] will be sufficient. 
^°The Gribov- Zwanziger confinement scenario will be presented in more detail in Section 1.4. 

The energy independence of the ghost propagator follows from the fact that the Faddeev-Popov operator 
involves only spatial derivatives and spatial components of the gauge fields (see Refs. [27, 29] for a 
complete derivation and discussion). 
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This divergence appears in any number of dimensions and cannot be regularized by usual 
dimensional regularization, although when taken as a full set in the Dyson-Schwinger 
equations, all such loops cancel. To handle this type of divergences, Leibbrandt in- 
troduced a modified form of the dimensional regularization, the so-called split dimen- 
sional regularization [24], in which two complex parameters u and a are introduced, 
d^q — (P'^q and dqo — )• d'^'^qo with the limits uj ^ 3/2 and cr — t- 1/2 to be taken after all 
the integrations have been completed. An alternative approach is the so-called negative 
integration method (NDIM) [31, 32], where a "Feynman-like" integral is solved, i.e. a 
loop integral in negative D-dimensional space with propagators raised to positive powers 
in the numerator. With these two methods the Coulomb gauge integrals were studied 
up to one and two-loop level perturbatively and results for the divergent part for several 
of them have been achieved. These divergences do in principle cancel order by order 
in perturbation theory (tested up to two-loops [33]), but this cancellation is difficult to 
isolate. 

Furthermore, within the first order formalism we are able to cancel the unphysical ghost 
fields, i.e. the Faddeev-Popov term. This means that we reduce the functional integral 
Eq. (1.20) to "physical" degrees of freedom, the transverse gluon and transverse vf fields, 
which in classical mechanics would be the configuration variables and their momentum 
conjugates (see below). 

We keep the term "physical" into quotations marks because it is realized that the true 
physical objects are the color singlet states, their observables being the mass spectrum 
and the decay widths. 

The presentation of the first order formalism, together with the reduction to the "phys- 
ical" degrees of freedom, follows [29] and is discussed in some detail. We start by express- 
ing the field strength tensor Eq. (1.6) in terms of the chromo-electric and -magnetic fields 
(recall that the temporal component A^°' has been renamed to cr°) 



ijk 



1 



VjAt - -gr^A^Al 



(1.28) 



such that the Yang-Mills action can be split into chromoelectric and -magnetic terms 



d^x 



—E"" E"" —B'^ B"' 

2 " 2 " 



(1.29) 



Next, we consider the chromoelectric term in the action. We linearize this term and hence 
convert to the first order formalism by introducing an auxiliary field vf via the identity 
[12] 



exp<; i / d^x]^E''-E'' 



Vtt exp < i / d'^x 



1 



-TT -TT 



(1.30) 



Classically, the vf field is interpreted as the momentum conjugate to A. The tt field is then 
split into transverse and longitudinal components using the identity 



const 



V {(P,t} exp 



(1.31) 
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After making the change of variables vf — )• vf — V0 and cohecting together all the parts 
that contain tt, we can write our full functional integral as (<I> denotes the collection of all 
fields): 



Z 



j P$ exp {iSg + ISym + iSpp} = / exp {iSg + iS^ + iSs + iSpp} (1.32) 



with 

Sb 
Sfp 



d X Qa 



d^x 



i7 Dq + \'y ■ D — m 



1 



Q/3 



2 



(1.33) 



Having derived the functional integral and the full QCD action in the first order formal- 
ism, we are now in the position to show that the ghost loops indeed cancel and the system 
reduces to the "physical" degrees of freedom. In the next section, we will demonstrate 
that the above QCD action simplifies to an expression where only the transverse A and 
TT fields appear. 



Reduction to "physical" degrees of freedom 

Given the functional integral, Eq. (1.32), and Yang-Mills part of the action (the quark 
field is not considered in this discussion), Eq. (1.33), we start by rewriting the Lagrange 
multiplier terms as 5- function constraints, which automatically eliminates the V-A and 
V • vf terms in the action. In addition, we also rewrite the ghost terms as the original 
Faddeev-Popov determinant. This has clearly the drawback that the local formulation 
and the BRST invariance of the theory^^ are no longer manifest. The functional integral 
now reads 



/ 



P$Det 



V-Dd'^ix-y) 6 (V-A) S (V-n)exp{iS} 



(1.34) 



with 



d'^x 



2 2 2^^ 



1.35) 

where we have defined the effective color-charge of the gluons = f"''^^ A'^ -t:^ . Next, we 
use the fact that the action Eq. (1.33) has become linear in a (after introducing the field 
tt), and write the integral over cr as a 5- function constraint. This enforces the chromo- 
dynamical equivalent of Gauss' law giving 



Z = / P$Det 



V-DS^ix -y) 5[ V-A] 5 ( V-vf ) 5 ( - gp"" ] exp {15} 



(1.36) 



^^The standard BRST and time- dependent Gauss-BRST transformations used in this work will be intro- 
duced in Chapter 4. 
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with 



d X 



(1.37) 



The implementation of the Gauss' law is very important because this essentially ensures 
the gauge invariance of the system. Defining the inverse Faddeev-Popov operator M: 



-V-D 



ab 



M 



be 



we can factorize the Gaufi law 5-function constraint as 



6 { -V -D"^ ((>''- gp" 



Det 



V-Dd'^ix-y) 6 {(t)" - W'^gp^) . 



(1.38) 



(1.39) 



The inverse functional determinant cancels the original Faddeev-Popov determinant from 
Eq. (1.36), leaving us with 



Z 



V^d ( V-A) 5 ( V-7f ) (5 ( - W^gp^ ) exp {15}. 



(1.40) 



We now use the (5- function constraint to eliminate the field. Since the inverse Faddeev- 
Popov operator M is Hermitian, we can reorder the operators in the action to give us 



Z 



V<^6 (V-A) 6 V-7f exp{icS} 



with 



-gp''M'"'{-V'')M'"'gp'' + tt^-O^A 



,0 Xa 



(1.41) 



(1.42) 



As promised, the action Eq. (1.42) contains only transverse A and vf fields. All other 
fields, especially the unphysical ghosts, have been formally eliminated. However, the 
appearance of the functional 5-functions and the inverse Faddeev-Popov operator M have 
led to a non-local formalism. It is not known how to do practical calculations within this 
formulation, but the non-local nature of the above result certainly serves as a guide to 
the local formulation. 

Before we close this section, few more remarks are in order. Quite generally, one can 
argue that first order formalism in Coulomb gauge is better suited to describe physical 
phenomena then other gauges such as Landau gauge. Indeed, the natural decomposition 
of degrees of freedom, both physical and unphysical, inherent to the first order formalism, 
automatically leads to the cancellations of the unphysical components. The subtlety is 
then to identify how these cancellations arise, and also it is very important to ensure that 
approximation schemes employed respect such cancellations. For example, the unphysical 
ghost loop of the gluon polarization should be cancelled in the Dyson-Schwinger equations. 
Given that the ghost propagator is energy independent, the temporal component of the 
gluon propagator must itself have a part that is independent of energy [34] in order to 
cancel this divergence. Later on in Chapter 5 devoted to heavy quarks, this observation 
will be used to show that only color singlet quark-antiquark bound states are physically 
allowed, regardless of the specific form of the energy independent part of the temporal 
gluon propagator. 
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On the other hand, one unpleasant feature of this approach is the large number of 
fields, but as it turns out this does not have serious implications for the Dyson-Schwinger 
equations [29, 35]. More important though, the issue of renormalisability remains unclear, 
in the sense that a complete proof of full multiplicative renormalization is still missing 
and one does not have a Ward identity in the usual sense [12]. A brief overview of the 
attempts to renormalize Coulomb gauge will be given at the end of this chapter. 

1.3.2. Alternative methods 

Currently, the most popular continuum formalism to QCD in Coulomb gauge is the 
Hamiltonian formalism [36-42].^^ In this approach, one starts by imposing Weyl gauge 
Aq{x) = 0, and subsequently, the Coulomb gauge is fixed with the help of the Faddeev- 
Popov method, whereas the Gauss law is imposed as constraint. The Yang-Mills Schro- 
dinger equation is then solved using the variational principle for the vacuum state, with 
a Gaussian ansatz for the wave functional. Then the vacuum energy is minimized and 
this leads to a coupled system of non-linear Dyson-Schwinger equations for the gluon en- 
ergy, the ghost and the Coulomb form factor and for the curvature in configuration space 
[36]. These have been solved analytically in the infrared and numerically in the whole 
momentum regime [40]. Similar to Landau gauge (see [43] for a recent review), it has 
been found that in Coulomb gauge there are two different infrared powerlaw exponents 
for the gluon and the ghost propagator [41, 42]. The favored solution is the most singu- 
lar - with the ghost propagator dressing function diverging as 1/|A;| - which generates a 
linearly rising heavy quark potential at large distances [41].^^ Also the gluon energy is 
divergent in the infrared, reflecting the absence of the gluons in the physical spectrum at 
low energies, which is again a signal of confinement. Recently, the static potential between 
infinitely heavy color sources has been also studied with a method based on the Dyson 
equation for the Wilson loop. In [44], the authors considered the temporal Wilson loop 
with the instantaneous part of the gluon propagator and the spatial Wilson loop with the 
static gluon propagator, and solved the corresponding Dyson equation in Coulomb gauge. 
In both linearly rising potential has been found. 

A second approach to Coulomb gauge is the lattice QCD (see, for example, [45] for a 
review of the results obtained on the lattice). The first lattice calculations have concen- 
trated on the infrared behavior of the ghost and gluon propagators. In particular, in Ref. 
[34] it has been shown that the static transverse gluon propagator is suppressed in the IR 
limit, while the time-time component is enhanced. Moreover, in the infinite-volume limit, 
it has been found that the transverse gluon is well described by the Gribov's formula [13] 
but unfortunately this study was not conclusive in the ultraviolet [34, 46]. Recently, in 
[47] the residual temporal gauge has been fixed and the renormalization of the full gluon 
propagator has been studied. It has been found that the static propagator is renormaliz- 
able only in the limit of continuous time, i.e. the lattice Hamiltonian formulation, and the 
resulting static propagator satisfies the Gribov's formula at all momenta. For the ghost 
propagator, lattice results have been reported in [48, 49], and more recently in [50]. An 
infrared divergence stronger than has been found, in agreement with the horizon con- 
dition necessary in the Zwanziger confinement criterion. Very recently, lattice studies of 

^■'in Ref. [39], the results obtained within the Hamiltonian approach to Yang-Mills theory are reviewed, 
^''in the next section, the connection to Landau gauge as well as the infrared ghost dominance will be 
discussed in more detail. 
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the Coulomb gauge quark propagator have been also undertaken and preliminary results 
have been presented [50]. The ambiguities due to the Gribov copies have been analyzed 
and it was found that their influence on the quark propagator is small. Moreover, the 
residual gauge has been fixed, and it appears that the impact of the residual gauge fixing 
is only reflected in the time-dependence of the propagator. 



1.4. Connection to Landau gauge and aspects of confinement 

Due to its covariance. Landau gauge d^Alj^ = has been for a long time the preferred 
gauge for non-perturbative Dyson-Schwinger studies. It preserves Lorenz invariance and 
it has a distinct property which makes it very attractive for practical approximations, 
namely that the the ghost-gluon vertex remains bare in the infrared^^. 

The first non-perturbative calculations in Landau gauge date back to the late seventies, 
when the infrared gluon propagator was first studied by Mandelstam and Bar-Gadda 
[52, 53]. In this calculation, the ghost loop and the four-gluon vertex, which do not 
contribute at first level in perturbation theory, have been neglected and only the (non- 
abelian) triple gluon vertex has been considered. It was found that the gluon propagator 
is infrared divergent and moreover, assuming a single gluon exchange, a linearly rising 
potential between heavy quarks has been derived. Today this picture is known as infrared 
slavery. In Landau gauge this picture has been rather misleading, in the sense that only 
very late it has been shown that in this gauge it is not the gluons, but the ghosts that drive 
the infrared properties of the theory More precisely, the gluon propagator vanishes at 
zero momentum, whereas the ghosts provide for a long range correlation. The gluon and 
ghost propagators are characterized by the so-called infrared exponents, which are related 
by a scaling relation, hence the name scaling solution. We briefly mention that there exists 
a second solution, the decoupling solution, which possesses quite different characteristics: 
the ghost propagator is not infrared enhanced but remains bare in the infrared and the 
gluon propagator becomes finite instead of going to zero [56, 57]. The connection between 
the two solutions (scaling and decoupling), different only in the deep infrared, still remains 
to be understood. 

Let us now analyze the consequences of restricting the configuration space to a compact 
region (Gribov or fundamental modular region), for the infrared properties of the theory. 
As discussed in the beginning of this chapter (Sec. 1.2), the ultraviolet regime is not 
affected by this restriction (when the coupling becomes small, all relevant configurations 
lie in the vicinity of gA = 0), whereas the infrared can be in principle governed by any 
domain within the Gribov region. Zwanziger argued that only the behavior of the gauge 
field on the Gribov horizon (i.e., the boundary of the Gribov region) is important for 
the infrared properties [58]. Pictorially, the situation is similar to a compact sphere of 
radius r in high N dimensions, where the probability distribution is concentrated on the 
boundary due to the volume measure r^~^dr. In Coulomb gauge, the Coulomb potential 



Or at least, from a semiperturbative analysis of the Dyson-Schwinger equation for the ghost-gluon 
vertex (with nonperturbative ghost and gluon propagators and a bare ghost-gluon vertex), it follows 
that deviations from the tree-level vertex are very small [51]. 
We refer to the original works [54, 55], and for a review, see Ref. [43]. 
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Chapter 1. QCD in Coulomb gauge 



between two external color charges^'^ has been derived [41], and it is related to 



<M(-V2)M>"* {x,y) 



(1.43) 



where M is the inverse Faddeev-Popov operator, defined in Eq. (1.38). Since on the 
Gribov horizon the Faddeev-Popov operator develops a zero eigenvalue, the Coulomb 
energy Eq. (1.43) becomes very large in the vicinity of the Gribov horizon, leading directly 
to an asymptotically linear potential. This mechanism is known as the Gribov-Zwanziger 
scenario of confinement. One can also inspect the relation between confinement and 
the restriction to the Gribov region by examining the infrared behavior of the ghost 
propagator. This propagator is given by to the expectation value of the inverse Faddeev- 
Popov operator and therefore is strongly divergent, due to the vanishing eigenvalues of the 
Faddeev-Popov operator on the Gribov horizon - this is known as the horizon condition 
[13]. Hence, the ghost propagator becomes infrared enhanced from the effects of the 
Gribov horizon. 

As outlined above, in Landau gauge the infrared ghost dominance has been established, 
but in Coulomb gauge, the Gribov-Zwanziger scenario has a somewhat different realiza- 
tion, depending on the specific formalism. As discussed in the previous section, in the 
Hamiltonian formalism the ghost propagator is infrared enhanced, whereas in the func- 
tional formalism, apart from the ghost and the transversal spatial gluon propagator, there 
is a third propagator, the temporal gluon propagator (in the canonical formalism, this 
would correspond to the non-abelian color Coulomb potential, Eq. (1.43)). Assuming that 
this propagator is largely independent of energy and diverges like 1/k^ in the infrared (as 
indicated by the lattice results [49]) one directly obtains a linearly rising potential for 
color singlet quark-antiquark states, leading back to the old infrared slavery picture^^. 
We also mention that very recently, a relation between the ghost and temporal gluon 
propagators has been found [60]. Based on the Gribov-Zwanziger scenario, this result, 
together with the Slavnov-Taylor identities presented in Ref. [27], provide an important 
element towards connecting the infrared slavery with the ghost dominance picture. 

1.5. Issue of renormalizability 

As already mentioned, the renormalizability of Coulomb gauge in the continuum has not 
been proven yet. In the following, we briefly review the efforts made in this direction. 
Among the various attempts, the most sophisticated approach has been pursued in Ref. 
[61]. There, the authors define the so-called "interpolating gauge" 



and recover Coulomb gauge in the limit a — )• 0. A linear shift in the field variables is 
performed in order to exhibit a symmetry (called r-symmetry) between the Fermi and Bose 
unphysical degrees of freedom. Individual closed Fermi-ghost loops and closed unphysical 
Bose loops diverge like 1/y/a, but they cancel in pairs at every order in perturbation 
theory by virtue of the r-symmetry. Thus in the Coulomb gauge limit the correlation 

^'''This potential is renormalization group invariant and is an upper bound for the gauge invariant potential 

from the Wilson loop - a statement known as "no confinement without Coulomb confinement" [59]. 
^*We postpone the detailed presentation of this mechanism to Chapter 5, concerning the heavy quarks. 



adoAo + V -1=0 
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1.5 Issue of renormalizability 



functions are finite, and this remains true also for the renormahzed correlation functions. 
However, there have been also identified one-loop graphs that vanish like \/a, which do 
not exist in the formal Coulomb gauge (i.e. for a = 0). These graphs can not be neglected 
since they give a finite contribution at two-loop order, when inserted into the graphs 
that diverge like Ij^fa. One possibility is that these graphs are merely gauge artifacts 
and decouple from the expectation values of all gauge-invariant quantities such as Wilson 
loop, but up to know this has not been explicitly shown. 

Renormalization of Coulomb gauge QCD has been studied also within the Lagrangian, 
second order formalism. In Ref. [62] , a proof of algebraic renormalizability of the theory 
has been given with the help of the Zinn-Justin equation. Through diagrammatic analysis 
the authors have shown that in the strict Coulomb gauge g^D^Q {Dqq is the time-time 
component of the gluon propagator) is invariant under renormalization, in accordance 
with a similar result obtained by Zwanziger [12]. In a covariant gauge, no component of 
the gluon field has this property. 
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Chapter 2. 

Dyson-Schwinger equations 



Dyson-Schwinger equations are the equations of motions in quantum field theory (analog- 
ous to the classical Euler-Lagrange equations) and they relate the various Green's func- 
tions of the theory. They are very powerful tools to treat nonperturbative phenomena, 
such as confinement and chiral symmetry breaking, whereas in the weak coupling regime 
the perturbative series is recovered. The most convenient way to derive Dyson-Schwinger 
equations (which will also be employed in this work) is to use the functional method, i.e. 
to derive these equations directly from the invariance of the generating functional under 
the variation of the field [10, 11]. An alternative method is the Dyson resummation [9], 
which reorganizes perturbative corrections into subdiagrams. 

Dyson-Schwinger equations built an infinite tower of coupled non-linear integral equa- 
tions, providing a complete description of the theory. Thus, in order to solve the theory 
it would be in principle necessary to solve the whole set of equations. However, this is in 
practice impossible and hence the main question is how to truncate the system, i.e. how 
to find a way to reduce these equations to a smaller subset of simpler equations which can 
be solved. The only systematic truncation relies on perturbation theory, otherwise one is 
forced to make an ansatz for the unknown higher order Green's functions. Importantly, 
the ansatz must respect the symmetry properties of the theory, i.e. it must obey the 
Ward-Takahashi identities in QED or the Slavnov-Taylor identities in QCD. 

In this Chapter we present the formal derivation of the Dyson-Schwinger equations of 
the two- and three-point quark Green's functions (i.e., the quark gap equation and the 
quark-gluon vertex functions). In Chapter 3 we will then present results in the perturb- 
ative limit, and in the second part of this thesis we will derive the four-point Green's 
functions and analyze them in the limit of the heavy quark mass. 

2.1. Field equations of motion 

The full generating functional of the theory is constructed from the functional integral, 
Eq. (1.32), by adding the corresponding source terms. Explicitly, we have (recall that P$ 
denotes the integration over all fields): 




(2.1) 



with the action Eq. (1.33) and the sources defined by 




(2.2) 
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In the derivation of the quark field equation of motion (from which the Dyson-Schwinger 
equations will be derived) the gauge- fixing term in the action, Spp, and the terms arising 
from the conversion to the first order formalism, discussed in the previous Chapter, are 
unimportant because the quarks are not connected by a primitive vertex to any of the 
corresponding fields, including the ghosts (i.e., there is no direct coupling term in the 
quark Lagrange density). What is however important later on is that these extra fields 
will formally enter the discussion of the Legendre transform (through partial functional 
derivatives) which, in principle, gives additional terms but which will turn to be vanishing 
at one-loop order perturbatively. 

Also, it is important to note that the generating functional, Eq. (2.1), is restricted to 
the Gribov region. This restriction might generate complications due to the presence of 
the Gribov copies inside the Gribov region. However, as discussed in Section 1.3.1, the 
Gribov copies do not influence the derivation of the Dyson-Schwinger equations in the 
continuum. Moreover, the boundary terms that may in principle appear are identically 
zero due to the fact that the Faddeev-Popov operator vanishes on the boundary of the 
Gribov region. 

The quark equation of motion follows from the generating functional Eq. (2.1) and from 
the observation that the integral of a total derivative vanishes: 

J ^^'^YM + i j d'^x[ Qxa {il^Do + ij ■ D - rnj ^q^p 

+Xxaqxa + QxaXxa ] + ...>= 0. (2.3) 



In the above, we have inserted the explicit expression for the quark contribution to the 
QCD action. Also, we have written the quark sources explicitly and denoted the rest with 
dots. Using the expression for the components of the covariant derivative, Eq. (1.23), it 
follows that 

y P$ I ^i-f'>do, + i7 • V, + 9T"7°f^x - gT'i ■ K - m] q.p + Xxaj exp {15} = 0, 

(2.4) 

where S is the full action plus source terms. This expression can be rewritten in terms of 
derivatives of the generating functional Z: 



SZ r m 6^Z 



i7%x + i7-V,-m -;^+[gT^jO] , , y re y - 



6'Z 



+ XxaZ = 



(2.5) 

Th equation Eq. (2.4) is the starting point for the derivation of the quark Dyson-Schwinger 
equations. Before we proceed to explicitly derive them, let us first introduce some nota- 
tions and briefiy review the various Green's functions of the theory. 

In general, the vacuum expectation values of the time-ordered products of field operators 
- the full n-point Green's functions of the theory (both connected and disconnected) - 
are obtained by functional differentiation of the generating functional with respect to the 
sources: 

G„(xi,...,x„)= y (2-6) 

dlJ[Xi) . . . dlJ{Xn) 

However, in practice we work with connected and one-particle irreducible n-point 
Green's functions. In order to eliminate the disconnected vacuum to vacuum diagrams. 
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2.1 Field equations of motion 



we use the generating functional of the connected Green's functions W, defined as 

Z[J]=e^W, (2.7) 
such that the connected Green's functions are given by 

dlJ[Xi) . . . OlJ[Xn) 

We now introduce a bracket notation for the functional derivatives of W with respect to 
the sources, such that for a generic source Jq, 

<iJa>= (2.9) 
OlJn 



Explicitly, we have: 

6Z[J] 
5^Z[J] 



Z[J] <iXax>, (2.10) 

ax- ^[-^l i<'Px^Xxa> + <Vx><mxa>] ■ (2-11) 

^^Pi^^Xxa 

Using the above equations, we convert Eq. (2.5) into derivatives of W[J] and obtain : 
i7°5ox + if • - m <iXxi3> +9T'' { 7° [<ip><ixxa> + <ipliXxa>] 

- afi 

-1^ [<^Jkx><^Xxa> + <U^Mxc>] } + Xxa = 0. (2.12) 

We define the generic classical field (since no confusion can arise, we use the same notation 
for the quantum fields which are integrated over and for the resulting classical fields) to 
be: 

\ f 1 sz 

^»=Z J ^^'^'-explicS} = (2.13) 

Explicitly, the classical quark and antiquark fields are given by: 

qaix) = ^ j V<^q^ix)exp{[S} =^31^ =<iXaix)> 

qa{x) = ^ Jv^q,{x)ew{iS} = -\^M^^ = - <ixaix)> . (2.14) 

Furthermore, we define the effective action (function of the classical fields) via the Le- 
gendre transform of W[J] with respect to the fields: 

^<t>,q,q] = W[J,X,X] - '^Ja(t>a - yXaQa - iqaXa- (2.15) 

In the above, we have explicitly separated the quark and Yang-Mills components, such 
that Jq and (pa denote generic gluonic (Yang-Mills ) sources and classical fields, respect- 
ively, and we also use the common convention that summation over all discrete indices 
and integration over continuous arguments is implicit. The generating functional r[$] 
then yields the n-point proper or one-particle irreducible ( IPI) Green's functions, which 
are those Green's functions that are still connected after one internal line is cut: 

r„(xi, . . . , xn) = . (2.16) 

di$(xi) . . . dl^{Xn) 
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In the following, we introduce a bracket to denote derivatives of T with respect to fields 
- although the notation is similar to the derivatives of W with respect to the sources, no 
confusion can arise since we never mix derivatives with respect to sources and fields. This 
gives: 



<ijQ,>= 



5W 

5\Jn 



and <i$rK>= 



5T 
5ii: 



(2.17) 



Note that care must be taken to observe the correct minus signs associated with the quark 
(Grassmann) fields and sources. Explicitly, this reads: 



Qaix) =<iXaix)>, Xa{x) = - <iq^{x)>, 
laix) = - <iXa{x)>, Xai^) =<^qa{x)> ■ 



(2.18) 



Having defined the proper (IPI) Green's function, we can now rewrite the equation of 
motion, Eq. (2.12), (and from which the Dyson-Schwinger equations will be derived) in 
terms of proper functions as: 



<nxa> 



a/3 



+ ffTV <wl\Xxa>] - gT^l^ [A%^qxa+ <iJLiXxa>] • (2.19) 



In a similar fashion, one can derive the gluon field equations of motion. They are given 
by^ (the trace is over Dirac and fundamental color indices): 



r5°Aia/3(i9^")(i9x/3) - 5Tr {[r"7*]„;3 <iXx/3iXxa>} 



dSd'^z-r'^^l''^^{z,x,y) 
Ac^^2^rfjJ(x,y,z) 



+ 



(2.20) 



<ifj"> 



.(0)a 
qqaa, 

|4„. j4 



(\qxa)iSqxp) +5Tr{[r"7 <iXx/3iXxa>} 
1 

2 



<UjyiJ^^> —lA'jyiAl^ 



d-yd''zT'^X^{x,y,z) <iJ^yip> -iA^yia] 



in''- 



1^ . 



+ 



(2.21) 



where T^'^^" r^^^'^ 
wnere i i gg^^ 



and T^^\^a'^-,T^^I^^!' , T^}'^'^^ are the tree- level quark-gluon and triple- 



crAAij^ aAai ' HAijk 

gluon vertices, respectively. The first two terms in the above gluon equations of motion, 
containing the quark-gluon interactions, are needed for the derivation of the quark con- 
tributions to the gluon proper two-point functions. The rest of the terms represent the 
Yang-Mills self-interaction, and are required for the derivation of the Dyson-Schwinger 
equations for the quark-gluon vertex functions (the dots represent the terms which are not 
important for the quark sector of the theory and have been left aside). We also mention 



^The complete derivation, carried on in the context of the Yang-Mills studies, has been presented in 
Ref. [29] and will not be repeated here. 
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2.1 Field equations of motion 



that in principle one can derive another two equations of motion, for the vf and (j) fields, 
but since the quarks do not directly couple to any of these fields, the quark field will 
not give a contribution to the corresponding proper two-point functions (at least at one 
loop-perturbative order). 

At this stage, it is useful to introduce multiple functional derivatives with respect to 
quark fields and sources, which will be later on used to derive the Dyson-Schwinger 
equations. Consider the following partial differentiations, both with respect to sources 
and fields: 

<X{J)> = -i5[7] <\^fii^^><U^X{J)> (2.22a) 



(5i$_a 

5 
5\Jr 



<y($)> = i5[7] <iJp\J^><\^^Y{^)> (2.22b) 



where ^[7] = ±1 accounts for the fact that the quark fields are Grassmann- valued, i.e. a 
minus sign appears when the index 7 refers to the following combinations 

<. ..ig^><iX7...>, <. . . ix^Xig^ . . .> . (2.23) 
For X{J) = iJa, we have that 

± i^T^ <iJa>= ±5[7] <i$^i$7><iJ7iJc,>= 5^/3 (2.24) 

(the overall sign is negative for $a = Qa)- Taking the functional derivative of this with 
respect to the source iJs and using the relation Eq. (2.22b), we find 

-S['j] i3i^^><U^iJa> 



5Us 

= i5[7, k] <UsUf^><i^K^i3i^^><iJ.yiJa> +'nsi3S[j] <i^i3i^^><U.yUsUa>= 0, 

(2.25) 

where the factor risjs = —1 if the fields 5, /3 anticommute. For X{J) = iJs'iJa, Eq. (2.22a) 
becomes 

— — <UsUa>= -i5'[7] gi^ ^><iJ^UsUa> (2.26) 
or, with the help of Eq. (2.25) 

<UsUa>= -mfsSh, k] <iJ5iJK><i$K$/3i$7><iJ7iJa> • (2.27) 



61^ IS 

With this notation multiple functional derivatives of arbitrary order can be efficiently 
constructed. In particular, this will be used in the next Section for the derivation of the 
quark gap equation. Also, later on in Chapter 6 the above formula will be employed in 
the derivation of the 4-point quark Green's functions. 
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Figure 2.1.: Full nonperturbative diagram for the quark self-energy. Filled circles denote dressed 
propagators and empty circles denote dressed vertices. Springs denote connected 
(propagator) functions, solid lines denote quark propagators and wavy lines denote 
the external legs of the proper functions. 



2.2. Quark gap equation 

We start the derivation of the gap equation by taking the functional derivative of the 
quark field equation of motion, Eq. (2.19), with respect to the quark field ig^, and omit 
those terms which will eventually vanish when the sources are set to zero. We arrive at: 



d yd z 5{x — y)5{x 



(2.28) 

where the (configuration space) tree-level quark proper two-point function T^^^ (x) and 

quark-gluon vertices T^^j^, ^^^g^j obtained from the quark equation of motion Eq. (2.19). 
In this thesis we will consider the gap equation both at one-loop perturbative order and 
in the heavy quark limit. Consequently, the explicit form of the corresponding tree- level 
quantities will be given separately in Chapter 3, which deals with perturbation theory, 
and in Chapter 5, where the heavy quark limit is investigated. 

We now use the formula Eq. (2.27) to calculate the functional derivatives appearing 
in the bracket. As already mentioned, simply because of the presence of the vf and (j) 
fields arising in the first order formalism, we must allow for the additional terms to be 
generated. For completeness, we keep all these terms for the moment, bearing in mind 
that they vanish when we consider the one-loop order in perturbation theory: 

<i/)^iX27> = - J d^d'^u KiXzf^XvSX^Pppix^qvS^qwficrly 

-J d'^vd'^u <iXzyiXvS><ippi<''u><^qvsiqwi3^(l)i>, (2.29) 
<UjyiXz-f> = - j d'^vd'^u <iXz'yiXv5><^JyjiJtk><^WqwiB'iA^uk> 



6 



I 



d'^vd^u <iXz'r^XvS><'^Jy^K^k><'^qv5^qwi3^'n-ik> ■ (2-30) 



At this point, it is useful to introduce our conventions and notations for the Fourier 
transform. For a general two-point function (connected or proper) which obeys transla- 
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tional invariance we have: 

<iJa(y)iJ/3(x)> = / dkWa^{k)e-'''<y-''^ (2.31) 



<icl>,(y)i$^(x)> = J dkT^^ik)e-"'<y-^\ (2.32) 

where dk = d^k/{2TT)^. The propagator (connected 2-point function) Wa/3{y,x) and 
proper (IPI) two-point function Ta/siy, x) are related via the Legendre transform. Whereas 
in covariant gauges this is simply an inversion, in Coulomb gauge this may not always 
be the case. The relation between the connected and proper two-point functions follows 
from Eq. (2.24). For the quark propagator, we find (in momentum space) the standard 
relation 

Wc,^{k)r^p{k) = 6^p. (2.33) 

Returning to the gap equation, we insert the expressions Eq. (2.29), Eq. (2.30) into 
Eq. (2.28), Fourier transform into momentum space, and we obtain the quark Dyson- 
Schwinger (or gap) equation: 

+ J dL.rgt",^(fc,-c.,^-A;)1^5,^5(a;)r|^,5^(^,-A;,A;-o.)T^,t(^— ^) 
+ 1 dL.rgt%(/c,-a;,a;-A;)1^5,^5(a;)r|,^5^(o.,-A:,A:-o.)l^,"^(A:-a;) 

+ / duTf^2^a,(k,-u:,u;-k)W^,,s{u^)T'ggAJ5^^^ 

+ / durf^2ia,ik,-u,u-k)Wg,,s{u^)T'^,.^spi^ (2.34) 

The self-energy corrections are presented diagrammatically in Fig. 2.1. We see that the 
TT and (f> fields do make a contribution thanks to the existence of the mixed propagators 
WAnij and Wa(j) in the first order formalism. But, as emphasized, these contributions will 
drop out at one-loop order because of the absence of corresponding tree-level vertices, 
i.e., there exist no direct interaction terms in the action between the quark fields and the 
auxiliary fields of the first order formalism. However, for future studies one has to bear 
in mind that additional contributions may arise. 

In the second part of this thesis we will also consider the gap equation in the heavy 
quark limit. In this case we shall work in the standard, second order formalism, where 
the auxiliary fields do not appear from the first place. 



2.3. Quark contributions to the gluon propagators 

In order to understand the analytic structure of the gluon propagators, it is necessary to 
explore the quark contributions to the gluonic proper two-point functions.^ In contrast 

^Apart from the quark contribution, the gluon propagators contain ghost loops, as well as a collection 
of terms generated by the tree-level 4-gluon interactions (which give rise to tadpole and explicit 2-loop 
contributions to the gluon Dyson-Schwinger equations). All these will not be considered in this work 
(a complete derivation of the Yang-Mills terms appearing in the Coulomb gauge gluon propagators can 
be found in Ref. [29]). 
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A ,o 




Figure 2.2.: One-loop diagram for the quark contributions to the gluon proper two-point functions. 

Fihed circles denote dressed propagators and empty circles denote dressed vertices. 
Solid lines denote quark propagators and wavy lines denote the external legs of the 
proper functions. 



to covariant gauges, the various gluonic degrees of freedom (temporal and spatial) are 
being separated into three proper two-point functions, T^j^j, Taa and Tg-A- The Dyson- 
Schwinger equation for the third function, T^jA, can be derived, since the quark-gluon 
vertices entering the equation are defined, even though the function itself does not have 
a tree level component in the first order formalism. 

Starting with the a equation of motion Eq. (2.21) and following the same procedure as 
for the gap equation we derive the quark contribution to the proper two-point function 
Fo-cr (indicated by the index (g)) in configuration space: 

<icj"icr^>(g)=-Tr j d^y(fzd\d*vTfJ^^^{z,y,x) <iXyy'^Xui3><^qu/3^qvs'icrl^><^XvS^Xza> 

(2.35) 

In the above, the trace over Dirac and fundamental color indices is taken. Performing the 
Fourier transform, we get in momentum space: 



(k) = -Tr y" dLvVf^^^^iuj - k, -w, k)Wgg^p{u})T'^^^psi^, ^ - ^> -k)WggSai^ - k). 

(2.36) 

Similarly we obtain: 

KAiiq)ik) = -Tr J dLvVf^^^^iu - k, -w, k)Wqgyfi{uj)r''ggAifisi^^ k-iv, -k)WqqSai(^ - k), 

(2.37) 

(it is easy to check that T^Ai{q) = '^Aai{q)) and 



^AMmik) = -TrJ du;rf>2,^^{u;-k,-u;,k)Wqq,piuj)T'ggAM^^k 

(2.38) 

In Chapter 3 we will consider these loop contributions (shown collectively in Fig. 2.2) at 
one loop perturbative level, and compare the results with the calculations performed in 
covariant gauges. We can already anticipate that since at one loop perturbative level the 
quark loop cannot be different from its covariant analog (the difference is that the spatial 
and temporal degrees of freedom are separated into the corresponding proper two-point 
functions), the one loop results should equal the covariant gauge calculations. Later on 
we will see that this is indeed the case. 
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(a) (b) (c) (d) (e) 






(i) 



Figure 2.3.: Dyson-Schwinger equation for the (spatial or temporal) quark-gluon vertex, written 
in terms of proper (IPI) Green's functions. Bfobs denote dressed vertices and all 
internal propagators are fully dressed. Internal propagators denoted by springs may 
be either spatial (A) or temporal (cr) propagators, dashed lines represent the ghost 
propagator and solid lines represent the quark propagator. Symmetry factors and 
signs have been omitted. 



2.4. Quark-gluon vertex 

The quark-gluon vertex, as an important component that relates the Yang-Mills and the 
quark sector of QCD, has been intensively studied in the last years, mostly in covariant 
gauges. At perturbative level, it has been analyzed in Ref. [63], in arbitrary (covariant) 
gauge and dimension. Also, in Landau gauge, nonperturbative studies have been carried 
on (see, for example, [64] and references therein). 

In Coulomb gauge, due to the intrinsic noncovariance, there are two quark-gluon ver- 
tices, spatial and temporal.^ Just as for the quark contributions to the gluonic two-point 
functions, the Dyson-Schwinger equations for the spatial and temporal quark-gluon ver- 
tices are obtained by taking the functional derivative of the gluon field equations of motion 
Eq. (2.20), Eq. (2.21) with respect to iqvS, '^Qw-y and setting sources to zero. 

Defining the Fourier transform for the vertex functions (all momenta are incoming): 

T{x,y,z) = J dki dk2 dk3 {27r)^5iki + k2 + k3)e-''''-''-'''^-y-'''' 'T{ki,k2,k3), (2.39) 

we get for the spatial quark gluon vertex, in terms of proper Green's functions: 

■^In fact, within the first order formalism one can derive the Dyson-Schwinger equations for two more 
vertex functions, corresponding to the interaction of the quarks with the additional fields vr and cf>. 
However, since there are no direct interaction terms in the Lagrangian, these terms do not give a 
contribution at one loop order in perturbation theory. 
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(a) (b) (c) (d) 

Figure 2.4.: Momentum rooting for the quark-gluon vertex functions. Diagrams (a) and (c) rep- 
resent the QED-hke graphs, and diagrams (b) and (d) arc the non-abehan graphs. 



^Uia^{ki,k2,h) = -gT^pf 

^^\qAj^p{k^ + ^2, kl - Uj)W^j^Aji,{ki - w) 

+^qaa-y{^ - kl,ki, -Uj)Wgg^5{i^)rfgMSxik3, ^ , -ks - Uj)WggXn{k3 + w) 

xr|^,^^(/c3 + k2, kl - uj)Wt{ki - UJ) 

+r^g^fc„7(^ -ki,ki, -uj)Wggys{uj)Tlg^sp{uj, k2, -k2 - '^)WZ{k2 + ^) 

^G^aA%{k2 + ^, k^, kl - U)W'A^,{U - kl) 

+'^qqaa-yi^ -kl,ki, -U})Wgg^s{i^)^gAkSl3i^ , ^2, "^2 " U})WZ^j,,{k2 + Uj) 

x^'''A%ik2 + UJ, ks, kl - o.)l^i!^(a; - ki) 

+'^1qAla^iu^ -ki,ki, -u)Wgg^siu})T^ggAm.Spiu}, ^2, -k2 - Uj)W''^Amk{k2 + w) 

xrS?fc(^2 + ^, A:3, kl - u:)W'A,i{uJ " ki) 

+^q.a',{uj -kl,ki, -Uj)Wgg^s{uj)r1g,5^{u, k2, -k2 - Uj)WZ{k2 + Uj) 
XrSt? (^2 + ^, ^3, kl - U:)W'„1{UJ - kl) 

+Tf^,,^p{u + ^3, k2, ki,u)W^,l{uj)Tf^;^{u, ks, -u - k3)W^t{u; + ^^3) 
+ mqqaps-yiu; + k^, kiM, (a;)^^)^,^ (o., A^s, -o; - hi)Wll^^{u: + ^3)} 

+ . . . . (2.40) 
Similarly, the temporal quark gluon vertex is given by: 

v'gg^^pikiMM = gTipi"" 

{TlgAka-yiuJ - kl,ki, -Uj)Wgg^s{uj)Tf^^g^{k3,UJ, -ks - U})WggXrjik3 + w) 

x^qqAjv^iks + UJ, k2, kl - uj)W^^Ajkiki - uj) 

+^qaa'yiuj - kl, kl, -Uj)Wgg^s{uj)rf^^^^{k3,UJ, -k^ - Uj)WqgXr,{k3 + u) 

x^qqav^iks +u;,k2,ki- w)VF^:^(/ci - Uj) 

+^qAkajiuJ - kl, kl, -u;)Wgg^siuj)rig^si3iu}, k2, -k2 - uj)WZ{k2 + uj) 

XrS? (^2 + UJ, ks, kl - Uj)W'Ajk{u^ - kl) 
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+ ^qqqaps-yi^ + ^3, ^^1,^2, -u;)!^!;^, (a;)^;^',^ (u;, h, -u - fc3)^|,V(^ + /ca)} + . . . . 

(2.41) 

In the expressions Eq. (2.40), Eq. (2.41), the dots represent the interaction of the quarks 
with the vr and (p fields, which are not regarded here. Further, the (two-loop) Yang-Mills 
diagrams arising from the 4-gluon vertex are not considered in this work (see also footnote 
at the beginning of Section 2.3), but for completeness we have included these two-loop 
contributions in the graphical representation from Fig. 2.3. The momentum rooting for 
the vertex functions is shown diagrammatically in Fig. 2.4. 

In both equations, the external gluon leg is connected to a bare vertex in the loop dia- 
grams. Alternatively, one can start with the quark equation of motion and take functional 
derivatives with respect to the gluon fields. Then in the corresponding Dyson-Schwinger 
equations the external quark legs are attached to the bare internal vertex. In the full 
theory, both equations should give the same quark gluon vertex, however in a truncated 
theory one of these equations is eventually easier to solve. ^ In this work we will only use 
the first "version" of the quark-gluon vertex, based on Eqs. (2.40, 2.41). 

In the next chapter, we will explicitly evaluate the divergences of the vertex functions at 
one-loop perturbative level. Moreover, the corresponding one-loop expressions, combined 
with the so-called quark-ghost kernels, will be then used to show that the Slavnov-Taylor 
identity for the quark-gluon vertex is satisfied at leading order in perturbation theory, 
without explicitly evaluating the integrals. 



* See also Appendix B of Ref. [64] for an extended discussion of this topic. 
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Chapter 3. 



One-loop perturbative results 



This chapter is concerned with the hmit when the coupHng g between quarks and gluons 
is small, i.e. the whole formalism can be expanded in powers of giving rise to the 
perturbative expansion. 

We will first derive the Feynman rules for the quark sector of the theory (basically, these 
are the lowest perturbative order Green's functions), and we will establish the general form 
of the two-point functions. Then we will consider the proper two- and three-point Green's 
functions derived in the preceding chapter at one-loop order in perturbation theory, and 
evaluate the corresponding two-point dressing functions. Coming across the problems 
originating from the energy-divergence and the noncovariance inherent to Coulomb gauge, 
we will derive the required noncovariant massive integrals, using a technique based on 
differential equations and integration by parts. We then consider the renormalization 
of the quark mass and propagator and we verify that the corresponding renormalization 
factors agree with the results obtained in covariant gauges. We also evaluate the first 
coefficient of the perturbative /5-function, and again we find that our results agree with 
the covariant gauge calculations. Moreover, we will consider the divergent parts of the 
temporal and spatial quark-gluon vertices, and shortly discuss the implications for the 
Slavnov-Taylor identity for the quark-gluon vertices presented in Chapter 4. 

3.1. Feynman rules 

In this section, we derive the basic Feynman rules and collect all the tree-level quantities 
required for our one-loop calculations. In addition to tree-level propagators and proper 
vertices we also derive the proper two-point functions of the theory. 

With the field equation of motion written in the forms Eq. (2.4), Eq. (2.12), we can 
now derive the Feynman rules for the quark components of the theory. We first derive the 
quark propagator. From the quark equation of motion in terms of connected functions, 
Eq. (2.12), ignoring interaction terms and functionally differentiating we get the tree-level 
propagator in configuration space 



For the quark propagator, the Fourier transform Eq. (2.31) explicitly reads (recall that 
translational invariance is assumed): 



= i7°9os + i7 • V^. - m 



<iX7(^)iX/3(a^)>^ ' -i(^^Q(5(z - x). 



(3.1) 



J 0/3 




(3.2) 
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such that in momentum space, we get 



: 

The solution is: 



'ke-''^<^-y^ {w^X^k) [jX - ^% + m]^^ + i6p,] . (3.3) 

In a similar fashion, starting with the quark equation of motion in terms of proper 
functions, Eq. (2.19) the tree-level quark proper two-point function is derived. We obtain: 



r(°) fx) 



i7°5ox + if • Va; - m (3.5) 

- a/3 



or, after Fourier transforming to momentum space: 

= i b'ko - yfc, - m] 5^p. (3.6) 

Due to the noncovariance, we have written out explicitly the components of but later 
on (where appropriate), we will use the usual notation ft = j^ko — ^ki- 

The (spatial and temporal) tree-level gluon propagators needed in this work have been 
derived in [29] and are given by: 

wPM-ik) = S^'§^,, W!^r\k) = 5^'^^^ (3.7) 

where 

tij{k) = 5ij-kikj/}? (3.8) 

is the transverse spatial projector. It is understood that the denominator factors involving 
both temporal and spatial components implicitly carry the Feynman prescription, i.e., 

1 1 



{kl - k^) {kl -k^ + iO+) ' 



(3.9) 



such that the analytic continuation to the Euclidean space (/cq ~^ 1^4) can be performed. 
This will be explicitly verified at one-loop order in perturbation theory. 

There are two tree-level quark-gluon vertices, spatial and temporal, again obtained by 
taking the functional derivatives with respect to quark, antiquark and gluon field: 

r£l% = (3.10) 



r£l4%, = -[<7TV].«. (3.11) 



Later on, in the evaluation of the vertex functions, we will also need the tree-level 
gluonic vertices derived in [28] (all momenta are defined as incoming and momentum 
conservation is assumed): 

ri°]t.(pa,P.,Pc) = -^gr'^s.^ipl-p^), 

^SAiifc (PaiPfcjPc) = -ifi'/"'"' i^ijiPa - Ph)k + &jk{Ph " Pc)i + 5ki{Pc - Pa)j\ , 

^f^Mip-c^PcPA) = -igr'^Pci. (3.12) 
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3.2. Two-point functions 

In this section, we introduce the general decompositions of the two-point functions derived 
in the preceding chapter, and derive the one-loop perturbative expressions for the asso- 
ciated dressing functions. In order to evaluate the noncovariant massive integrals arising 
in the one-loop calculations, we employ a technique based on differential equations and 
integration by parts. We then give the perturbative results for the two-point functions. 

3.2.1. General decomposition 

In order to investigate the Dyson-Schwinger equations for the quark gap equation derived 
in the previous chapter, in addition to the tree-level forms given in the previous section, 
we will also require the general decompositions for the quark propagator and proper two- 
point function, and the relationship between them. Because we work in a noncovariant 
setting, the usual arguments must be modified to separately account for the temporal and 
spatial components. Starting with the quark propagator, we observe that (recall that 
denotes the integration over all fields) 

T^gga/3(fc°, k) ~ 5^^ j V^qq exp {15} (3.13) 

such that under both time-reversal and parity transforms, the propagator will remain 
unchanged - the bilinear combination qq is scalar. Since the propagator depends on both 
hP and k, it has thus four components, in distinction to the covariant case where there 
are only two (a dressing function multiplying ^ and a mass term). Hence we can write 

Wqq^p{k) = 6ap—^ {ko-f'^Ftik) - ka'Fsik) + M{k) + koka^^'Fd{k)] (3.14) 

k^ — k^ — 

where all dressing functions are functions of both /cq and k"^. At tree- level, one can 
trivially identify Ft = Fs = 1, F^ = Q and M = m. The last term with F^ has no 
covariant counterpart. The possible appearance of this term will be discussed below. 
For the proper two-point function, the same arguments apply and we write 

r^ga^ik) = i6a^ {koJ^Atik) - ka'Asik) - Bra{k) + k^ka"" Ad{k)] (3.15) 

and we will refer to At^Ag and Bm as the temporal, spatial and massive components, 
respectively. Again the last component A^ has no covariant counterpart. The relation- 
ship between the connected and proper two-point functions is supplied via the Legendre 
transform (introduced in Chapter 2) and we have 

^qq{k)Wgg{k) = 1. (3.16) 

Let us now discuss the possible appearance of the genuinely noncovariant term corres- 
ponding to the dressing function A^ (or equivalently Ffi). At one-loop order in perturba- 
tion theory, these terms are vanishing. This can be deduced from the Dirac structure of 
the self-energy loop, stemming from the quark propagator and tree- level vertices. Namely, 
the tree-level quark propagator does not contain a term with k^ki^y^^^ and from the tree- 
level vertices we only have either two or two 7* matrices together (the gluon propagator 
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is either purely temporal or spatial). This implies that there is no one- loop contribution 
that has the overall structure 7^7* and in turn this means that = at one-loop order 
perturbatively. Hence, the components and F^i will only appear (if at all) at two-loop 
order or beyond. In in fact, very recently lattice calculations have shown that = Ad = 
[50]. 

Using the definitions Eq. (3.14), Eq. (3.15), together with the relation Eq. (3.16), we 
then get the following set of relations for the remaining dressing functions: 

Ft = Fs = M = ^^.3.17 

i.2 42 _ 42 _ d2 1.2 42 _ 11,2 42 _ d2 p42_ 1.242 _r2 ^ ' 

We again emphasize that these relations only hold up to one-loop perturbatively — in 
possible future studies it must be recognized that the fourth Dirac structure, 7^7*, may 
enter in a nontrivial fashion. In this case, the set of relations Eq. (3.17) should be replaced 
with the following set (which includes the functions A^ and Fa): 



M 





{kl - P 


- 'm?)At 


^o^t 


- pAj - 


' + ^Ik'^A^ 




{kl - P 


- m'^)As 




- Pai - 


' + ^Ik'^A^ 




{kl - P - 


- m'^)Bm 




- Pai - 


' + ^Ik'^A'^ 




{kl - P 


- m'^)Ad 




- Pai - 


' + ^Ik'^A^ 



(3.18) 



3.2.2. One-loop perturbative expansions 

Let us now consider the one-loop perturbative expansions of the quark gap equation and 
the quark contributions to the gluon two-point functions. Although so far the formalism 
has been presented in 4-dimensional Minkowski space, in order to evaluate the resulting 
loop integrals we have to convert to Euclidean space. This means that we make the 
analytic continuation /cq — )• iA;4, where k4 denotes the temporal component of the Euclidean 
4-momentum, such that fe^ = k\+P. Additionally, to regularize the integrals, dimensional 
regularization is employed^ with the Euclidean space integration measure 

= (2^ ^^-^^^ 

where d = 3 — 2e is the spatial dimension. To preserve the dimension of the action we 
must assign a dimension to the coupling through the replacement 

5' ^ 9^^i', (3.20) 

where ^ is the square of a non- vanishing mass scale (which may be later associated with 
a renormalization scale). 



^An alternative method proposed by Leibbrandt is the so-called split dimensional regularization. This 
has been discussed in the introductory chapter of this thesis. 
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The perturbative expansion of a generic two-point dressing function is written as: 

r = r(°) + /r(^) (3.21) 

where the factor is included in F^^^ such that the new couphng and T^^) are dimen- 
sionless. 

Let us first consider the (fuh nonperturbative) gap equation, Eq. (2.34) (depicted in 
Fig. 2.1). We first insert the various tree-level vertices and propagators given by Eqs. (3.4), 
(3.7) (3.10), and (3.11). Then we insert the general decomposition of the proper two-point 
functions, Eq. (3.15), occurring on the left-hand side of the gap equation, take the Dirac 
projection, solve the color and tensor algebra and lastly perform the Wick rotation. The 
one-loop temporal, spatial and massive components of the quark gap equation in Euclidean 
space read (the Casimir invariant Cp = {N^ — l)/2Nc is listed in the Appendix A): 



2 g „ 1 /" J J k^ijJA k4U}4{d — 1) 



1 /■ r 2[k ■ {k - uj)][uj ■ {k - Co)] k-Lo 



Mk) = i-^vcF-i y rf^l 



(a;2 + m?){k - ujY{k - ujY (w^ + m?){k - 
I fe-^(3-d) ] 

B„W = ra + „..VC. / { ^ _ + } . (3.24) 

As mentioned earlier, the possible contribution corresponding to the genuinely noncovari- 
ant dressing function does not appear at one-loop. To evaluate the integrals occurring 
in Eq. (3.23), it is helpful to use the identity: 

k-oj = ]^[k^ +00"^ -{k-ojf]-kiOJA, (3.25) 
which enables us to rewrite Ag as a combination of more straightforward integrals: 

^2^2 2{k'^ + m'^)kAUA 



A -, 2 1 /■ J f 1 (A:2+m2)2 



2 (jj'^[{k — a;)2 -|- m?]uj'^ (jj'^[[k — a;)2 -|- m2]a;2 

2£fc2 + 2kl 2fc-a;(l-g) 1 m2 + 3A:2 \ 

^ [{k - w)2 + m2]tj2 + p _ a;)2 + ^2^2 2 [{k - + m?]^'^ I • ^ 



Let us now examine the quark contributions to the various proper two-point gluon 
dressing functions given by Eqs. (2.36-2.38) (presented in Fig. 2.2). Again, we insert the 
tree-level factors given by Eqs. (3. 4), (3. 7), (3.10) and (3.11), solve the color and tensor 
algebra and perform a Wick rotation. The one-loop integral expressions are (recall that 
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we have Nf flavors of identical quarks): 

k'r'X)^^) = ,^Nj2 I + + f , (3.27) 

<Ta{q)\ ' J J {uj^ + m?)[{k - ujY + m?] ^ ' 

hk^T^^ dk) = f,'Nf4 [ du Tr'tr ^'"^t ^, (3.28) 

Ptij{k)r^^^^^^^{k) + A:j/cjT^]^^^^(A;) 

= 2I\ffi / ctu 2u77 \T~, 21 ' (,3.29) 

where T^j^j.^^ and f ^j^^^^ are the transversal and longitudinal components of the proper 

two-point function r^]|^^.^^-| (given by Eq. (2.38)), respectively^. 

Having derived the one-loop perturbative expressions for the propagator dressing func- 
tions and the quark contributions to the gluonic two-point functions, we can now proceed 
to evaluate the corresponding loop integrals. There are two categories of integrals arising: 
those which can be solved using standard techniques (such as Schwinger parametrization 
or Mellin representation — the details of these techniques are presented in the Appen- 
dices B, C), and those which require a more complex approach. In the next section, we 
will concentrate on this later variety. Since this part is rather technical, the reader might 
skip this and go directly to Section 3.2.4, where the physical results are presented. 

3.2.3. Noncovariant massive loop integrals 

The noncovariant massive loop integrals appearing in the one-loop expansions from the 
previous section are studied by using a technique based on differential equations and 
integration by parts developed previously in Ref. [30]. We will consider the two integrals: 

A^iklk'') = J ^2p_^)2 + ^2p2' (3-30) 

At.iklk'') = j ^2p_^)2+^2p2- (3-31) 

3.2.3.1. Derivation of the differential equations 

Let us first write Eqs. (3.30) and (3.31) in the general form [n = 0, 1) 

/"(A;|,A;^)= / „ -^-^ „ . (3.32) 

J ^'^[{k — bj)'^ + m'^\uj'^ 

In this derivation k^ and k"^ are treated as variables whereas the mass, m, is treated as a 
parameter. The two first derivatives are: 

,91" _ [ du u2 f , k,ik,-u,) \ 

I . (3.34) 



" dk^ J uj'^[(k — Lv)"^ + m'^]uj'^ \ (k — oj)"^ + m'^ 
dr f du uj2 in k-{k-uj) 



dkk J uj'^[{k — uj)'^ + m'^]uj'^ | {k — oj)"^ + m"^ 



The decomposition of the spatial two-point function TAAij is explained in detail in Ref. [29] 
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There are also two integration by parts identities: 



d 



n+1 



duo uj2 



n + 1 -2 



/ 



.to 



(A; — a;)2 + 
(J- (tl; — k) 



up' {k — lS)^ + rrP 



(3.35) 



(3.36) 



Adding these two expressions gives 



bP\{k — u;)2 + rn^Jw^ 



d+n- 3- 2 



a; • (cj — A;) 
(fc — w)^ + 



(3.37) 



Expanding the numerator factor, we can rewrite Eq. (3.37) as 



dbj ufl 



d + n - 4 + 



k — LO + 
{k — uj)^ + m? 



(3.38) 



Combining this with Eq. (3.35) then yields: 



/ 



dio uj2 I 2 ^^^^^ ~ 

io'^[{k — oo)"^ + m^]cJ^ \ (fc — io)'^ + 

duj uj2 



Lo'^[{k — ujy + m'^]uj'^ 



2/C4 



doooo2 



k"^ + m? 
2 



{n + d-A)-n + l 



+ 



k^ + 
k"^ + m? 



uj^[{k - + rrP] J uo'^[{k - ooy + m'^]'^' 



dL0L02 



doooo'2 



[{k — Lo)"^ + m?Y(jo'^ 
(3.39) 



This leads to the temporal differential equations for Am and A. 



4 . 

rrf 



dA^ 

^ dki 



kS^ 
dkd 



1 + 2 



{d - A)kl 
k"^ + m"^ 
dco 



Am + 2 



dco 



k + m 

k"^ + J [{k — w)^ + m^j^w^ 
dio 



a;4[(fc-a;)2 + m2] 



(d-3)fe|^, ^ ^ P + m2 
k"^ + A;2 + 

duo UJ4 



uo^[{k — uo)"^ + m?] J Lo'^[{k — uo)'^ + m'^] 



Lo'^[{k — uo)"^ + m2]2 ' 

duo U04 
[{k — uo)"^ + m2]2a;2 
rfo; W4 

t212 ■ 



(3.40) 



(3.41) 
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In the same manner, we derive the differential equations involving the spatial components: 



dAn 

' dki 



* dki 



2-d + 2 



+ 2 



{d - 4)fc^ 

1^2 _|_ 



{d - 3)P 



duj 



d + 2- 



+ 2 



du 



uj^[{k - uY + m?] 



+ 2 



2fc^ /■ duj uii 



(3.42) 



fc2 + 771,2 



[(A; — + 7772]2d;2 



doj bJi 



(jj'^[{k — ojY + "T-' 



212 ■ 



(3.43) 



It is in fact possible to write down a mass differential equation, but in the light of the 
method presented here this would not bring any new information. However, this third dif- 
ferential equation will turn to be useful in checking our solutions - the detailed derivation 
will be presented in the Appendix C. 

At this point, it is instructive to show how the differential equations for the massless 
integrals considered in Ref. [30] are regained. In the massless limit, there are potential 
ambiguities arising in the integrals appearing in Eqs. (3.40 -3.43), because in part, the 
limits 777 — )• and e — t- do not interchange. Let us start by considering the following 
integral given by Eq. (B.16) (similar arguments apply to all the integrals appearing in the 
differential equations): 



duj 



up- \{k - w)2 r772]2 (47r)2- 



r(3/2-e) 



2F1 l,l + e;3/2-e; 



777^ 



(3.44) 

It is useful here to invert the argument of the hypergeometric with the help of the formula 
(see, for instance, Ref. [65]): 



2Fi(a,6; c; t) 



r(c)r(6 - g) 

r(6)r(c-a) 

r(c)r(a - h) 



{-ty^Fi ( a, 1 - c a; 1 - 6 a; - 



r(a)r(c-6) 



{-ty^iFx 6,l-c + 6;l-a + 6;- 



(3.45) 



Then we have: 



fc2 (47r)2-^ 



r(i-e)r(-e)r(l + e) 



(4vr) 



2-e 



2e) 



m 



1 77T,2 

2F1 ( l + e,- + 2e;l + e;— ^ 



(3.46) 



In the expression above, the problem of the non-interchangeable limits is seen explicitly 
in the first term. However, when all the integrals occurring in the various differential 
equations are put together, such terms explicitly cancel and only the second term of 
Eq. (3.46) (which leads to the correct massless limit) contributes. 

Returning to the differential equations, we evaluate the standard integrals in terms of 



34 



3.2 Two-point functions 



e (see Appendix B) and with the notation x = k^, y = k'^ we obtain for Am- 
dAm 



2x- 



dx 



1 - 2(1 + 2e) 



X 



+2 



m 



X + y + m 
y -\- m 



Ar, 



r(-e)r(l + e 



X2F1 l,l + e;3/2-e;-^ 



r(2 - e) 



-2F1 ( 2,l + e;2-e 



X + y 



Y2F1 1, 1 + e: 2 - e 



X + y 



(3.47) 



2y- 



dAr, 



dy 



-l + 2e-2(l + 2e) 



y 



X + y + 



Ar^ 



\m 



r(-e)r(l + e 



l,l + e;3/2-e;-^ 



r(2 - e) 



-2F1 ( 2,l + e;2-e 



x + y 



Y2F1 1, 1 + e: 2 - e 



x + y 



(3.48) 



and for the integral A"^ = k^A^ 
dAm. 



2x- 



-1 -4e- 



X 



Ar, 



dx [ X + y + m''' 

+2^7TATr— <! ^^2^1 l,l+e;3/2-e;-^ 



(4vr: 
2-e 



2-e 



X + y + 



2F1 ( 2,l + e;3-£;-^^) +(l-e)2Fi (l,l+e;3-e; 



(3.49) 



2?/- 



5y 



-l + 2e-4e- 



y 



X + y + 



^71 



\m 



2]-l-e 



y 



(47]-)2 e + y + m? 



X2F1 l,l+e;3/2-e; 



Y 



2-e 



2F1 ( 2, 1 + e; 3 - e; ) + (1 - e)2Fi ( 1, 1 + e; 3 - e; ^ ^ ^ 



(3.50) 



where 



r(i/2-£)r(i + £) 

r(3/2-e) 



Y 



r(i-£)r(i + £) 
r(2-£) 



(3.51) 



3.2.3.2. Solving the differential equations 

Let us first consider the integral Am- By the same method as in Ref. [30], we make the 
following ansatz: 

Amix, y) = FAm{x, y)GAm{x, y) (3.52) 
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such that 

dFA, 



2x 

FAm2x 



dx 

dFAm 

dy 

dGAm 



l-(2 + 4e)- 



X 



X + y + m 
-l + 2e-(2 + 4e) 



FAn 



y 



X + y + 



Faui, 



(3.53) 
(3.54) 



\m 



2l-l-e 



dx 



y + m^ / y 
tA \2~£ \ ~< \ o'^^-c'i I i, i + e; 3/2 - e; 2 



r(-g)r(i + £) 
r(2-£) 



-2F1 ( 2, 1 + e; 2 - e 



x + y 



+ Y2F1 [ l,l + e;2-e;-^^ 



(3.55) 



FAm'^y- 



dG 



Am 



\m 



2l-l-e 



dy 



y 



+- 



(47r)2 £ 1^ 2; _l_ y _l_ 77^2 
r(-e)r(l + e 



X2F1 l,l + e;3/2-e;-^ 



r(2-£) 



-2F1 ( 2, 1 + e; 2 - e 



X + y 



+ Y2F1 [ l,l + e;2-e;-^ 



(3.56) 



By inspection, it is simple to determine the solution for the two homogeneous equations, 
Eq. (3.53) and Eq. (3.54): 



(3.57) 



Since the mass m is treated as a parameter, the (dimensionful) solution, Eq. (3.57), may 
have an integration constant proportional to [m^]~^~^. However, returning to the original 
equations, Eq. (3.53) and Eq. (3.54), we see that the only consistent solution is the one 
for which this constant vanishes. 

Let us now make the following ansatz for the function GAm-, which will be verified below 
{z = x/y): 

GAm{x, y) = G\^{x, y) + GAm{z). (3.58) 

The component G\^{x,y) can be found by adding the differential equations Eq. (3.55) 
and Eq. (3.56), which lead to: 



9G^A dG^A 
X — TT^ + y- 



1 + ^ + 1 

y 



dx 



dy 



^ I m" V " ' y ' " 

(47r)2-- I ^xiy + m^) y^J^T^ - 1 



+ Oie) } (3.59) 



Because the function is multiplied by the function FAm (which does not have an e 
pole), the term of order 0{e) will not contribute. The solution of this equation is: 



^Am(^) y) 



(47r 



l2-£ 



' 1 + — In 



y 



Iny 



1 + ^-1 

y 



+ Oie) } +Ci. 

(3.60) 
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Before we proceed to determine GAm, we justify the ansatz for GAm, given by Eq. (3.58). 
First we observe that: 

2xFAm^^ = 2zFAm^ + 2xF^^^. (3.61) 
ox oz ox 

Then we subtract the above equation, Eq. (3.61), from Eq. (3.55). This gives: 

dGAn._ dG\^ 1 Mi^f^+^i^^ p . , ^.,/2 y 



9z -FAm (47r)2 ^ [ X + y + V ' ' ' m? 

( x + y\ r(-e)r(l + e) / x + y 

\ m"' / r(2 — e) \ 

(3.62) 

Evaluation to the first order in e is straightforward and we see that the right hand side 
of the above expression is only a function of the variable z. This allows us to write down 
a differential equation for GAmiz) in the form: 

dGAm 1 1 f 1 2 



from which we get immediately 

GAm{z) = -jv^TT^ | - - 7 + lnm2 + 0(e) j +C2. (3.64) 
{Airy ^ [ e J 

Returning to the original differential equations (3.53 - 3.56) with the function G{x, y) = 
G^^{x,y) + GAm{z), we see that the only consistent solution is the one for which the 
overall constant Ci + C2 vanishes. 

We may now put together the solutions Eqs. (3.57), (3.60) and (3.64) and write for the 
function Am- 

(x + ?/ + m2)-i-^ r 2 f x + y + 771^ 

Am[x,y) = TTZTw <^-- + 27 + 21n' 




+ 0{e)}. (3.65) 

V y VV^ + ^-V 

We see that for m? = we regain the result from Ref. [30] and that the singularities 
are located at x + y + m? = (with m?,y > 0). Two more useful checks arise from the 
study of the power expansion around x = and the mass differential equation (this will 
be explained in detail in Appendix C). 

We now proceed in the same fashion to determine the function Am{x, y) = F-^^{x^ y)^Am( 
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The resulting partial differential equations are in this case: 



2x- 



Am 



2y 



dx 

dFArn 

dy 



-1 -4e 



X 



-l + 2e-4e 



X + y + m? 

y 



^ Ami 



X + y + 



(3.66) 
(3.67) 



(4vr)2-e 
Y 



y + m^ / y 
:rX2Fi 1, 1 + e; 3/2 - e; ^ 



X + y + 



Fi 2, 1 + e; 3 - e; 



2-e 

+ (l-e)2Fi ( l,l+e;3-e;■ 



3; + y 
x + y 



(3.68) 



P-Am^y^ 



\m 



2i-l-e 



+ 



(47r)2-^ 
Y 



y 



X + y + 



X2F1 (i,i + e;3/2-e--A. 



Fi 2, 1 + e; 3 - e; 



2-e 

+ (l-e)2Fi ( l,l+e;3-e;■ 



x+_y 
x + y 



with X, y defined previously. The solution to the first pair is 

F-Ami^. y) = x"^l^y~^l^+'{x + y + m^r^\ 



(3.69) 



(3.70) 



For brevity, in the above expression and also in the derivation of the function = 
^\m Am (^^^ analogue of GAm) we omit the constants of integration - they vanish as 
in the case of the functions FAm and GAm- 
As before, for G^^{x,y) we make the ansatz: 

GAmi^iy)=G'jm(^iy) + GAmi^). (3.71) 

In the limit e — t- 0, the component G^^{x, y) is determined from the differential equation: 



dG^ dG^ 

^_Arn^y_ Am 



1 



m 



In 



y 



dx ' " dy (47r)2 = | x + y + m'^ + ^2 
The solution of this equation is: 

1 I., / Vl + 'f (\^z + \ 



y 



+ 0{e) 



G'-Aj-^y) 



(47r: 



i2-e 



iln 



In 



ix/i- 1 



-iLio 



l-i^/i 



1 + f -iv^ 



+ iLio 



1 + i\/i 



1 + f -iV^ 



+ 0{e) 

(3.72) 
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where Li2(z) is the dilogarithmic function [66]: 



U2{Z) 



ln(l - t) 



dt. 



(3.73) 



As before, we check that the ansatz for G^^^{x,y) given in Eq. (3.71) is correct and 
derive the differential equation for the function G^^, in the Hmit e — )■ 0: 



z 



1 



z+1 



dz (47r)2-^ 

The result we leave for the moment in the form: 

1 



[ln(l + z) - In z - 2 In 2] + 0{£] 



(3.74) 



„ { — 41n2arctan(-v/^) 

(47r)^ e 



dt 



Vt{i + t) 



ln(l + t) 



dt 



\nt + 0{e) }. 



r, , ^^^^ , ^v^; J • (3.75) 

With the solutions, Eqs. (3.70), (3.72) and (3.75), after some further manipulation we can 
write down the following simplified expression for the integral A^: 



-2 In 



+ y + ,„2)-l-.(l + ^ + mi 



(47r 



dt 



Vt{l + t) 



In 1 + t + 



m 



1 + ^ + 1 

y 



1 + ^-1 

y 







1 arctan 









with the integral 
dt 



In ( 1 + t + — 

y 



Viil + t) 

. l-i^/i^A 

+iln 



vr In 2 — i In 



1-i^ 
1 + i^/i 



2 In 2 



y 



u/z 



+iln 



In 



In 



1 + + 1 

y 



'1 + 



-iln {^/z — i) 



In 2 + In (1 + z + — ) - In (1 - - - In (^i - i) 



1 i 



iLi2 - - -\/z + iLi2 - + - a/z 



1 i 



iLio 



i + yz 

-'l+x/z 



+ iLi2 



-i + yz 

i + x/z 



+iln(V^+i) 



m 



In 2 + In (1 + z + ) - In (1 + i^i) - - In ( + i) 



-iLio 



iJl + ^ + Vi 
-i + A/i 



iLio 



i + V^ 



-iLi2 



-l + xfz 



+ iLi2 



i./l + ^ + ^/i 



(3.77) 
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We see that for = we get the correct hmit for the function A^. We also mention 
that the singularities are located at x + y + m? = and the apparent singularities at 
z = —1 (i.e., X + y = 0) in the expression Eq. (3.76) are canceling out. This can be easily 
seen by making a series expansion of Eq. (3.76) around z = —1: 




y/z{z + 1 + 



+ 0(e). (3.78) 



Again, the result Eq. (3.76) has been checked by performing an expansion around x = 
and by studying the mass differential equation (see Appendix C). 

3.2.4. Results in the limit e ^ 

Having derived the noncovariant massive loop integrals in the previous section, and col- 
lecting the results for the standard massive integrals from Appendix B, we are now in the 
position to write down the perturbative results for the dressing functions under consid- 
eration. For the temporal, spatial and massive components of the quark gap equation, 
Eqs. (3.22), (3.24) and (3.26), we find in the limit e 0: 



(3.79) 




(3.80) 



Bm{k) = m + m 



(4vr; 



2-e 




+10 - 2 Jl + ^ In '\ - 2 fl + ^ fl + S ) + 0{e) \ , (3.81) 



where the function fm{x, y) is given by (x = k1,y = k'^, z = x/y): 



1 + ^ + 1 



2 fm \ 2 / V w 
fmix, y) = —p^ In ( I arctan {\fz) H -= In — | arctan 




dt , / 

■ In 1 + + . (3.82) 



Vi{l + zt) \ y 
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The last integral has been rewritten using the identity: 



1 dt ( m? 
\ -Tr 7 In 1 + t + 



dt 



In 1 + zt + 



m 



(3.83) 



As a useful check, we can set m = and show that the results for the temporal and 
spatial components are in agreement with the calculation performed independently using 
the one- loop massless integrals derived in Ref. [30]. 

Another important point is related to the singularity structure of the above dressing 
functions. As has been shown in the previous section, in the noncovariant integrals the 
singularities appear at x + y + m? = 0. Further, it is easy to see that the standard 
integrals have the same singularity structure. As promised, since the singularities in both 
the Euclidean and spacelike Minkowski regions are absent, we find that the validity of the 
Wick rotation is justified. In addition, the results for At,As and Bm can be compared 
(allowing for the color factors) and agree with those of Quantum Electrodynamics [67, 68]. 

Having calculated the dressing functions for the quark proper two-point Green's func- 
tion, we are now able to discuss the structure of the quark propagator. In Eq. (3.17) we 
first analyze the denominator factor. Let us denote (in Euclidean space): 



D{k) = klAl{k) + k^Ai{k) + BlXh). 



"2 a2i 



(3.84) 



Inserting the expressions from Eqs. (3.79), (3.81) and (3.81) into the above equation, we 
have: 



D{k) 



k^ + {l + Q 



+(A;2 + m'^] 



(4vr)2-£ 



(4vr; 



2-e 



1 , m"' 4 

--7 - In h- 

e fj, 6 

2 

m 



+ 3 



m 
fc2 



In 1 + 



7 - In h9 




^^2 ^^)| 



(3.85) 



Defining the renormalized mass, m/j, via: 

= Z^m% with = 1 - 6^£^^ |^ - 7 - In ^ ^| , (3.86) 

we see that the expression for D{k), Eq. (3.85), then contains explicitly the overall factor 
fc^ -|- m|j. This means that the simple pole mass of the quark emerges, just as it does 
in covariant gauges. For the remaining part, the singularity structure is such that non- 
analytic structures do not appear for spacelike or Euclidean momenta. Moreover, we see 
that the renormalization factor, Z^, which defines the physical perturbative pole mass and 
hence should be a gauge invariant quantity, agrees with the result obtained in covariant 
gauges [3]. 

Because of the Dirac structure, it is more convenient to write the quark propagator 
in Minkowski space, i.e. to analytically continue k^ — ?• —kQ. We have shown that the 
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analytic continuation of the functions At, Ag, Bm (and implicitly, of the function D{k)) 
back into the Minkowski space is allowed and this enables us simply to write: 

Wqq^p{k) = idc.^ {7°A:oAt(A:) - ^'hA^ik) + Bm{k)] D'\k). (3.87) 

Inserting the denominator factor, Eq. (3.85), in the limit e — )• and replacing the mass 
with its renormalized counterpart, the above expression gives: 

'-^^(4^0-^ 



Wqqo,p{k) = -Sal3— \ (^ + ma) 





+ finite terms 

(3.88) 



We can thus write down for the quark propagator: 



Wqqapik) = {-i)5a[Sj:^^^^rP^Z2 + &mte terms (3.89) 



1.2 _ 1,2 



and identify the renormalization constant (omitting the prescription dependent constants) 

Turning to the quark loop contributions to the gluon two-point proper functions, in 
evaluating the integral structure of Eqs. (3.27), (3.28) and (3.29) we observe the following 
relations (in Euclidean space): 



rS(,)(^) = ^%ik) = -^r« (^)(fe) = -pZqik) = i{klP), (3.91) 



where the integral I{k1, k"^) reads (using the results of Appendix B), as e — ?• 



2 P^ _ _Jlf 



Nf r 2 



(4^)2-e \ 3 



1 , k^ 

7 — In — 

e /i 



2"! 10 2 / in^ 

¥ + 3 V + 



I{ki 1 k 



The relations Eq. (3.91) are similar to the Slavnov-Taylor identities for the Yang-Mills 
part of the theory, derived in [28]. Also, the above integral Eq. (3.92) agrees with the 
results obtained in covariant gauges (see for instance [3]). This was in fact expected, since 
at one-loop level the quark loop as a whole is identical with its covariant counterpart — 
the only difference is that the various degrees of freedom (temporal and spatial) are being 
separated into the corresponding proper two-point functions, i.e., Taa-, ^Aa and ro-o-- 

Let us now consider the one- loop gluon propagator dressing functions, in connection to 
the first coefficient of the perturbative /5-function. Analogously to the two-point proper 
functions, these are constructed by writing D = D^^'^ + g'^D^^\ As mentioned previously, 
in the first order formalism we have to account for the presence of the additional vf , (p and 
ghost fields, and implicitly the corresponding propagators (for example Datt)- We have 
already seen that at one-loop the quarks only contribute to three of the gluon proper two- 
point functions {Taa, Tylo- and To-a)- But since these gluon two-point functions are related 
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to the various connected (propagator) two-point functions arising from the first order 
formahsm, there will be (quark loop) contributions to many more of these propagators. 
The relationship between the connected and proper gluon two-point functions in the first 
order formalism is derived in Ref. [30] and will not be reproduced here. The full set of 
quark contributions to these gluonic type propagators is given by: 

D%,)ik) = I?il\,)(A:) = T^X^k) = liklP), (3.93) 
= = = = ^(^4, k'). (3.94) 

We are now able to identify the first coefficient of the /3- function. As is well known in 
Landau gauge, a renormalization group invariant running coupling can be defined through 
the following perturbative combination of gluon and ghost propagator dressing functions 
[43]: 



1 fllN, 2Nf 



(3.95) 



At one- loop in perturbation theory, the coefficient of the 1/e pole above is simply minus 
the first coefficient of the /3-function (/3o = — llA''c/3 + 2A^j-/3). By inspecting the relations 
Eq. (3.93), containing the quark contribution to the propagator Daa , and those obtained 
in Ref. [30] for the Yang-Mills part of the propagator Daa and the propagator Dc, we 
see that the same result is achieved in Coulomb gauge. Moreover, in Coulomb gauge, a 
second renormalization group invariant combination of propagators appears and is given 
by g'^Dfjfj [12]. Again, combining our results Eq. (3.93) and those obtained in Ref. [30] 
we see that indeed the coefficient of 1/e agrees with this. 



3.3. Quark-gluon vertices 

We now evaluate the divergent components of the quark-gluon vertex functions at one 
loop perturbative level. Just as for the two-point functions discussed in the preceding 
section, we write the the perturbative expansion of the vertex function as 

r = r(°) + ^^r^i) (3.96) 

such that the new coupling and F^^^ are dimensionless. 

We first consider the temporal part of the quark-gluon vertex, Eq. (2.41). By the same 
method as for the two-point functions, we insert the appropriate tree-level vertices and 
propagators from Eq. (3.4), Eq. (3.7), Eq. (3.12), Eq. (3.10) and Eq. (3.11) and solve the 
color algebra. For the one-loop expansion of the temporal part of the quark-gluon vertex 
we find, in Minkowski space (the factor Cf — Nc/2 arising from the color algebra is derived 
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in the Appendix A): 




^ [7°(A;i -a;)o - 7*(^i - w)i + m] q [-7°(^2 + w)o + 7*(^2 + w)i + m] jtikioj) 
{ki — Lo)^ — m? (k2 + w)^ — u"^ 

[7°(fci - a;)o - fjki - uj)i + m] p [-7°(fc2 + uj)o + 7'(^2 + + m] oJ_l 

-uj)o-Y(h -cu)i + m] ^ (2^3 tjfc(a;) 
{ki - a;)2 - m2 ^ (w + k-i)^ 

_^ [7°(^ - ^2)0 - I'juJ - k2)i + m] fc (-2A:3 -a;)j tjkji^) 

; [7°(fci - a;)o - y (fci - (j), + m] ^ + ^s) t^K^) 1 

Since we are only interested in the divergent part of the above expression, we first 
make use of the identity Eq. (3.25) to simphfy the non-abehan terms containing two 
noncovariant denominator factors, and then we employ a simple power analysis to separate 
the convergent and divergent integrals (at leading order in perturbation theory). After 
eliminating the convergent terms and using the expression Eq. (3.8) for the transversal 
projector tij(k), we are left with the following divergent part of the temporal quark-gluon 
vertex: 



r£%(^i'^2,fc3) = (^-'-T^N^ + iT^Cj^j 

7''(7°a;g + 7*7V^i^^07^ fx ^j^k 



f f 7^(7°a;g + 7^7V^^a;;)7J / 
J \a;2pi-w)2-m2][(A:2+a;)2-m2] V^^ 



w2 



+ - 



a;2 [(^ki — a;)2 — m"^] [(/c2 + a;)2 — m'^] 



[(A;i - w)2 - m2](u; + A;3)2a;2 



^ _ UJkUJi _ {uj + ^3)^0;; 
2(^2 2(cJ + 4)2 



(3.98) 



Collecting the results for the divergent factors of the three-point loop integrals listed 
in the Appendix D, it is straightforward to calculate the divergent part of the temporal 
quark gluon vertex, in the limit e — >• 0. It is given by (in Minkowski space): 

.(i)<i _ rrd^fl Cf 1 
(47r)2 e 



r^i>. = T'l'jj^- + finite terms. (3.99) 



Repeating this calculation for the spatial component of the quark-gluon vertex, Eq. (2.40), 
we obtain: 

= (cf + ^iVc) (4^^ + finite terms. (3.100) 

The results Eq. (3.99), Eq. (3.100) will then provide a useful check for the Slavnov- 
Taylor identity for the quark-gluon vertices. We will return to this subject at the end of 
the next chapter, when we analyze the divergent structure of the Slavnov-Taylor identity. 
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Chapter 4. 

Slavnov-Taylor identity for the quark-gluon 
vertices 



The Slavnov-Taylor identities relate Green's functions with a different number of external 
legs, analogous to the Ward-Takahashi identities in QED. They play an important role in 
connection to nonperturbative Dyson-Schwinger studies, since they provide constraints 
for the higher n-point functions that enter the Dyson-Schwinger equations. 

In order to derive these identities, we start with the observation that the BRST invari- 
ance is an unbroken symmetry of the gauge fixed theory, at the level of the Lagrangian. 
The BRST symmetry (and its variation in Coulomb gauge) will be presented in detail 
in the first Section of this chapter. Next, we will present the formal derivation of the 
Slavnov-Taylor identities for the quark-gluon vertex functions, based on the BRST in- 
variance of the QCD Green's functions, and in particular, we will examine the so-called 
quark- ghost scattering kernels. We will then demonstrate that these identities are satis- 
fied at one-loop perturbative level. The proof is based on the translational invariance of 
the loop integrals and it does not require the explicit evaluation of the loop integrals. 

This chapter concerns the Slavnov-Taylor for the quark-gluon vertices, but similar 
identities can be derived for Green's functions of arbitrary order. Some of the higher 
order identities will be explored in the second part of this thesis, in the limit of heavy 
quark mass. Moreover, their implications for the Dyson-Schwinger equations will also be 
discussed. 

4.1. Gauss-BRST symmetry 

Regardless of the specific gauge, the gauge fixing term in the QCD action violates the 
invariance under the local gauge transformation Eq. (A. 3). Except for the axial gauge, the 
invariance of the full (gauge-fixed) action can be recovered by performing a local gauge 
transformation similar to Eq. (A. 3), but with the following ansatz [69, 70] 



where both c and A are Grassman- valued quantities and A is a constant. By demanding 
that the gluon, quark and ghost fields transform as 



(4.1) 



5A 



a 



9 



6q, 



'a 



-i[T-Uc-6Xq^ 
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dc" = -^/-^^^c^c^JA (4.2) 

(c and c are independent Grassmann fields, and thus we have independent transforma- 
tions), we find that the gauge- fixing term remains also invariant. This new symmetry is 
called BRST (Becchi-Rouet-Stora-Tyutin) and is the basis for deriving the generalized 
Ward-Takahashi identities (also called Slavnov-Taylor identities). 

While the "standard" BRST symmetry is independent of the gauge considered, in 
Coulomb gauge the time derivatives are absent and hence one can further define a time- 
dependent — so-called Gauss-BRST symmetry [12]. As before, the spacetime-dependent 
paramether 9^ is factorized into two Grassmann components 

r(x) = c"(x)5At, (4.3) 

but in this case 8Xt (not to be confused with the colored Lagrange multiplier A°) is a time 
dependent infinitesimal variation. Within the standard, second order formalism, the fields 
transform as: 

9 9 

5c^ = -X^SXt, Sc" = --r^^c^c^5Xt, SX" = 0. (4.4) 
9 2 

The Gauss-BRST transform, Eq. (4.3), is the starting point for the derivation of the 
Slavnov-Taylor identities in Coulomb gauge. These identities, together with the peculiar 
features introduced by the time dependent transformation, will be discussed in the next 
section. 



4.2. Formal derivation 

As discussed in the previous section, the full (gauge fixed) QCD action in the stand- 
ard, second order formalism, Eq. (1.21), is invariant under a Gauss-BRST transform, 
Eq. (4.3). The Coulomb gauge Slavnov-Taylor identities arise from regarding the Gauss- 
BRST transform as a change of integration variables under which the generating functional 
Eq. (2.1) (together with the source terms Eq. (2.2)) is invariant. Since the Jacobian factor 
is trivial [29] and only the source term varies, we deduce that 



S\t=0 



^^^li^ exp {iSqcd + i^FP + i5, + iSSs] 

= jv^ exp {iSQCD + -^Sfp + \Ss] j d^xdit - ^o) I - ^ + F'^p'^Ac''. 

1 1 1 1 

_ jo. \7 _|_ fO.bc ja Ab c ■— „arpa \n'^7! T"^ a/ _l _ i _ fohc—a h c I 

'Jix^tx'^x^J '^ix^ix'^x ^Xax'^x-'-af^'il^x '■'^x'tpx-'' l3aXax ~r ^x'lx ' 'Ix'-x'^xf- 

9 5 ^ J 

(4.5) 

Notice the 6{t — XQ) constraint, which arises because the time-dependent variation 6Xt and 
is characteristic to the Gauss-BRST transform. It leads eventually to a non-trivial energy 
injection in the Slavnov-Taylor identities which is not present in the covariant gauge case. 
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The above identity is best expressed in terms of proper functions. Using the definitions 
Eq. (2.17) for the derivatives of W with respect to the sources and of T with respect to 
the classical fields, we arrive at the identity 



d!^x5{t — xq) 



x<i(a2<i<>)cS-r'"=<ia> 



9 



(-V2)^ 



<1C^> 



+\T^p <iqax> [<iX(5xK> -c%qpx] + '-^^a [<iX/3xi??> +clqpx\ <^Qax> 



(4.6) 



Note that functional derivatives involving the Lagrange multiplier result merely in a trivial 
identity such that the classical field can be set to zero [27]. Further, to derive the quark 
Slavnov-Taylor identities, one functional derivative with respect to icf is needed and then 
the ghost fields/sources can be set to zero. Implementing this then gives 



= / d'^x6{t - xo) 



X |-- (52 <iai>) 6{z -x)- r''^ <ia> 
1 

+- 

g 



b xdc 1 



6ic'i ^ 



<U%yrg> -iA)j'"^5{z-x) 



6 



6idi 



<^Qax> 



(4.7) 



Two further functional derivatives with respect to iq^;^ and iqg^ are taken and all remaining 
fields/sources set to zero. Now we use the following relation, derived in Ref. [27] 



<iPx^Vx> 



did 



j=o 



(4.8) 



j=o 



and obtain the Slavnov-Taylor identity for the quark-gluon vertices in configuration space: 



d'^x5{t — Xq) 



52 



+ir^« <iqsJqax> 



> 



^_^2'J ^^9Sv^Q'yi^^A^^> 

00 — X] 



<i^ici> 



6'^''6{z - x)6s^6{v - x) 



j=o 



<^Qax^Q-yuj> 



(4.9) 



47 



Chapter 4. Slavnov-Taylor identity for the quark-gluon vertices 



We now introduce the following definitions: 



J=0 



^9T^I3 [- <^Xpx^Xe><iqe^q^uj'i'^X><Jx^JK><Vx^VT><^Crid 

+ <^Xfix'^XK><^qJq-fi^iCrici><irf^ir]r>^ , 



i^'K> 



^q;ccq5a i^y ^i^) — 



J=0 



<jjlirir><iCric'^i^K: 



= I <rf^\7]r> <ic^icf i$K> <i J^i Ja> <iqs^iqei^\> <\Xe^xpx> 

->]\gTl^. (4.10) 



In the above, the sources/fields J and $ refer to either A or a, and the internal indices 
refer to all attributes of the object in question (summed or integrated over). A few more 
comments concerning the above definitions and notations are in order. Recall that in 
the "standard" Dyson-Schwinger equations, the loop expressions contain a propagator 
(< 'yXaz^Xpx > or any other propagator) and an associated tee-level vertex, stemming 
directly from the interaction terms in the Lagrangian. However, as seen in Eq. (4.10), in 
the case of the (nonabelian) Slavnov-Taylor identities one has a different structure which 
arises from the Gauss-BRST transform (or generally from the BRST transform): namely, 
one has to functionally differentiate objects such as <ix^2.ifj">, but in this case the vertex 
arising from a direct interaction is missing. The resulting expressions still have a (partial) 
meaning as loop integrals, even though the tree-level vertex is absent. These types of 
would-be loop expressions, denoted with a tilde, are called quark-ghost scattering kernels 
and are common to the Slavnov-Taylor identities of nonabelian theories.^ 

With these definitions, we return to the configuration space Slavnov-Taylor identity 
Eq. (4.9), and rewritte it as 



d!^x5{t — xq) 



<i95i.ig7^io-x> ) 5{z- x) + 



V, 



+ <iq^^\qax> 



q;ccqa'y V 



(-V2) 



{x,z,uj) + igT^^6{z - x)6{uj - x) 



+ 



^'!,;ccq5a(.x, z, v) - igT^J{z - x)6{v - x) 



i;i'?7CJ>j' • 



(4.11) 



Using the formula Eq. (2.39) for the Fourier transform of the vertex functions, one 
can write the Slavnov-Taylor identity in momentum space (momentum conservation is 



^In Coulomb gauge Yang-Mills theory these objects have been studied in detail in [27], and in linear 
covariant gauges, some well-known examples are the ghost-gluon and quark-ghost scattering-like kernels, 
appearing in the identities for the three-gluon and quark-gluon vertices (a nice presentation of this 
subject can be found in Ref [71]). 
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Figure 4.1.: Quark-ghost scattering-like kernels. Solid lines denote quarks, dashed lines denote 
ghosts and spring lines denote gluons. All internal propagators are dressed. The cross 
indicates the absence of a tree- level vertex (see text for details). 



assumed) 



\^^,AaP^{klMM)Tt{-h) 

Kn 



~^^qqa5{kl 



f-- ( 

q;ccq\ 



^1 + qo,h -qo;h) +igT° 



5/3 



+ 



rg;ccg(^2 + qo, k-i - qo; ki) - igT'^ Tqqspi-k2) (4.12) 



In the above, Tec is the proper ghost two-point function and qq is the (arbitrary) energy 
injection scale that arises from the time-dependence of the Gauss-BRST transform. This 
expression is very similar to the one derived in Ref. [27] for the Yang-Mills case. 

We now write the quark-ghost kernels Eq. (4.10) in momentum space. After performing 
the Fourier transform, we arrive at the following expressions (represented diagrammatic- 
ally in Fig. 4.1): 



^■,ccqiPUP2,P3) 



^q;ccqiPl^P2,P3) 



igT^ I duW^{u)Wqq{pi-u:) 



ccqq 



r^^^(^,P2, -P2 - '^)W:;l{p2 + U:)T'qqx{Pl - ^,P3,P2 + Uj) 



ccqq 



r|;?^(a;,p2, -P2 - ^^)W^lip2 + u})T^g^{p3,Pi -u},p2 + uj) 



{UJ,P2,P3,P1 - W) 



W-^^iu;)Wqqiu;-piyigT^, 



(4.13) 



where the indices k. A, refer to the gluonic field types a or A (with the associated spatial 
index) . 

Just as in the Yang-Mills case, the above Slavnov-Taylor identity, Eq. (4.12) is used 
as a constraint, i.e. it helps to find meaningful truncations for the vertex functions of 
the theory, which are then used to solve the Dyson-Schwinger equations. In particular, 
in the limit of the heavy quark mass, this equation will be used to write the temporal 
quark-gluon vertex Tqq^ in terms of purely spatial, ghost or quark propagators and proper 
functions. This will then be used to obtain an exact solution of the quark gap equation 
and of the Bethe-Salpeter equation for quark-antiquark bound states. 

Using the Feynman rules presented in chapter 3, it is easy to show that the identity 
Eq. (4.12) is trivially satisfied at tree-level. Also, it can be verified that this identity is 
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fulfilled at one-loop perturbatively. In the next section, we will present this derivation in 
detail. 



4.3. Perturbative analysis 

Let us now briefly show that the Slavnov-Taylor identity is fulfilled at one-loop perturbat- 
ive level. The proof is based on the translational invariance of the loop- integrals, without 
explicitly evaluating them. At one-loop perturbative level, Eq. (4.12) reads: 

/C3 



15 



(^l)f ^;ccg(^l + QO, h - qo, k2) + f q;ccg(fc2 + 90, ^3 " ?0, ki)T\^{-k2). 

(4.14) 

In the above, we have inserted the tree-level spatial quark-gluon vertex, Eq. (3.11), and 
the tree-level ghost two-point function from Ref. [28] 



(4.15) 



Using the expressions Eq. (2.40), Eq. (2.41) for the quark-gluon vertices, we obtain for 
the (nonperturbative) combination of quark-gluon vertex functions, on the left-hand side 
of Eq. (4.14): 

klT^lil^{kiM, h) + kiTf^i{ki,k2, k3) 

'klrfl\ks,U, -k3 + k3,rf^i^{k3,U, -ks - UJ) 



xWqgik-s + w)r^g^j(A:3 + w, /ca, h - w)VFJ^j.fc(A;i 
-'^QqA-ki +uj,ki,-uj)Wgq{uj) k^rf^^{k3,uj,-ks -u) + k3ir'i^^^{k3,uj,-k3 -u) 



xWqqiks + io)rl,{k3 + to, k2, ki - u:)wt{ki - io 



-^qAk{-ki +^,ki, -u)Wqq{uj)rlg^{uj, k^, " ^2 " u:)WZ{k2 + ^)WTA^f,{-k^ + a;) 
^ kt^^aAajik2 + ks, ki-uj) + k3ir^°\^^..{k2 + uj, ks, ki - u) 
-T%A-ki +u;,k,, -u;)Wqq{u;)rl^A,{u;, fes, -k2 - u;)WX%kik2 + u;)W^''A-ki + a;) 
k3^aAaj(.k2 + UJ, k^, ki - Uj) + k3iT^^]^^l-{k2 + UJ, ks, ki - w) 



-^qAl{-kl + W, fcl, -Uj)Wqq{uj)T\^A^X^, /cs , " ^2 " ^)WXA^k{k2 + w) W^Iaj/ (" ^1 + 



UJ) 



kl^^aAAjkik2 + UJ, ks, ki-uj) + A:3ir3XJfe(fc2 + uj, ks, ki - u 

T%,{-k^ +UJ,ki, -Uj)Wqq{uj)r1(uj, , -k2 - Uj)WZ{k2 + Uj)W',l{-ki + Uj) 



{0)dbc, 



X A:3ir^^f.^(A;2 + uj, k^, ki - lo) 



(4.16) 
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The lengthy expression above can be simplified if we make use of the tree-level Yang- Mills 
Slavnov- Taylor identities derived in Ref. [27].^ For completeness we list them all here 
(together with the tree- level quark Slavnov- Taylor identity, which is trivial): 



(4.17) 



klTfl'{h,k2, k,) + k,,Tf^i{ki,k2, ks) = rfj {ki)igT'' - igT^TfJ{-k2) 

Further, in order to derive the one-loop expression, in the full nonperturbative equation 
Eq. (4.16) we insert the tree-level vertices and propagators given by Eq. (3.4), Eq. (3.7), 
Eq. (3.10), Eq. (3.11). Putting all these together, we arrive at the following simplified 
expression: 

k'ir^il'ik^, k2, k,) + kiTf^i^{k^,k2, k,) 

duo iy-i^Wqq{hi + uj)-imAAjk{ki - cj) - -t^Wqq{u:)-imAAjk{ki - w) 

+l''Wqq{k^ + u)-l''W„„{kl - W) - l''Wqq{uj)^''W„„{kl - Uj) 

kzjtij{uj) 



]:T''N,g^ 



duj ■ 



^2(fc3 



-f°{k3-Uj)o-Y{k3-Uj)k 



[7°(u; + k2)o - j'^iuj + k2)k -m] . 



{U! + k2 



+1 



.[j%oj + ki)o-7''iuj + ki)k + 



m\ 



7°(A;3-w)o-7''(A;3-w)fc 



(4.18) 



(cj + kiY — m? 

We first compare the first integral of Eq. (4.18) with the combination of one-loop quark 
two-point functions (without ghost factors) on the r.h.s. of Eq. (4.14). Using the expres- 
sion Eq. (2.34) for the gap equation and recalling that the vf and (j) fields do not give a 
contribution at one-loop order in perturbation theory, we deduce that these terms are 
equal and hence they cancel in the Slavnov-Taylor identity. For the rest of the terms (the 
combination of quark- ghost kernels on the r.h.s. and the remaining ghost term on the 
l.h.s. of Eq. (4.14)) we get: 



Tfhk^)ti-(kr + q^,k3 



<k2) + f g.gcg(^2 + go, ks - qo, ki)rf^ ( 



-k. 



.igT-f^ri')-'i(-k,) = '-T'^N.g^ f d 

Kin J 



k3jtij{u) 



-(7°fcio - -f'^kik - m) 



[7°(w + k2)o - -i^{uj + k2)k - m] 



+ 1 



(cu + k2Y ~ nT-"^ 
. [-f^uj + ki)o - + ki)k + m] ^ 



Tksk 



{u} + ki)"^ — m? 



(-7^20 + 7^2fc - m) - 2/i:3i-^2 



(4.19) 



In principle, the Yang-Mills identities also contain ghost-gluon scattering-like kernels, but at tree-level 
these factors do not give a contribution. 
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After rearranging the terms, we find that the second integral in the formula Eq. (4.18) 
equals the above combination of quark-ghost kernels and the ghost term. Collecting 
all the results, we deduce that the l.h.s. of Eq. (4.14) equals the r.h.s., and hence the 
Slavnov-Taylor identity for the quark-gluon vertex is satisfied at one-loop perturbative 
level. 

As a useful check, we can also show that the divergent parts of the quantities entering 
the formula Eq. (4.14) combine such that the Slavnov-Taylor identity is again satisfied. 
For this we use the results derived in Chapter 3, i.e. the divergent part of the quark 
gluon vertices Eq. (3.99), Eq. (3.100), and the divergent part of the quark gap equation, 
obtained by inverting the result Eq. (3.89) for the quark propagator 

C 1 

Vqq{k) = i7j-T2-(7°^o " l^h - 4m) + finite terms, (4.20) 

together with the ghost two-point function (previously calculated in [30]) 

Tf(k) = -i(5"''P--^- + finite terms. (4.21) 
^ 3 (47r)2 e 

Further, writing out the explicit form of the quark-ghost kernels Eq. (4.13) and making a 
simple power analysis, we notice that these factors do not contain divergences.^ Collecting 
the divergent factors Eq. (3.99), Eq. (3.100), Eq. (4.20) and Eq. (4.21), it is straightforward 
to verify that that the Slavnov-Taylor identity is again satisfied. 



■^The convergence of the quark-ghost kernels has been also shown in Ref. [72], in the framework of split 
dimensional regularization. 
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Chapter 5. 
Heavy quarks 

So far, we have investigated the Dyson-Schwinger equations (along with the Slavnov- 
Taylor identities) for two- and three-point Green's functions at one-loop perturbative 
level, where the coupling between quarks and gluons is small. In the second part of this 
thesis, we are concerned with the nonperturbative (i.e., strong coupling) region of the 
theory, with the aim to describe bound states of quarks, and interpret them in relation 
to the linearly rising potential which confines them. 

Traditionally, most of the nonperturbative studies are performed in covariant gauges. 
They concentrate on the light quark sector, where one is concerned with dynamical chiral 
symmetry breaking (a review of this subject can be found in Ref. [73], see also Ref. [74] 
for the phenomenological implications for meson and baryon states). So far, the Landau 
gauge studies of the heavy quark sector are unfortunately not conclusive, since the heavy 
quark propagator contains spurious poles which prevent a solution of the bound state 
Bethe-Salpeter equation [75]. The difficulty arises from the fact that nonperturbative 
Dyson-Schwinger equations and heavy quark limit are not automatically compatible, i.e., 
both the heavy quark mass and the ultraviolet-cutoff scale of the loop integrals (supposing 
that a ultraviolet-cutoff regularization is employed) tend to be the largest scale in the 
problem. 

An alternative approach is to reformulate the heavy quark sector of QCD as an effective 
theory - the Heavy Quark Effective Theory [HQET]^. In this approach, the main sim- 
plification stems from the fact that the effective quark Lagrange function is rewritten as 
an expansion in powers of 1/m, where m is the heavy quark mass. Naturally, when the 
quarks are very heavy one can keep only the first few terms in the mass expansion. The 
effective Lagrange function also enjoys two additional symmetries, which turn to be very 
useful for practical applications [77]: the heavy flavor symmetry^ and the so called spin 
symmetry, which arises because in the heavy quark limit the spin degrees of freedom do 
not couple to the gluon field. Moreover, the apparent incompatibility of the heavy quark 
mass and the ultraviolet cutoff has been recognized. To handle this problem, the so-called 
matching procedure has been developed [81]: QCD operators are written as series in 1/m 
via HQET operators and the coefficients in these series are determined by matching on- 
shell matrix elements in both theories. These calculations have been performed up to 
three loops in perturbation theory [84-86]. 

In the present chapter, the aim is to explore the heavy quark sector of the theory 
by combining nonperturbative Dyson-Schwinger equations with the heavy quark mass 
expansion [87]. We will employ the full nonperturbative QCD functional formalism, in- 

^ In [76-80] we collected the original papers; extended reviews can be found, among others, in Refs. [81, 82] 
and in the textbook [83]. 

However, as we consider identical quarks, this is not of particular interest for the present work. 
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volving the complete quark fields, rather than HQET expressions that refer to the heavy 
quark degrees of freedom. The effect will be that, due to the complexity of the equations, 
we will be obliged to restrict ourselves to leading order in the mass expansion, and this 
means that we will not be able to describe real quarks (which are not infinitely heavy). 
However, this does not alter our goal to study the confining potential, since all quarks 
should be confined, regardless of their mass. 

5.1. Heavy quark mass expansion 

Let us start by writing out the explicit quark contribution to the full QCD generating 
functional, Eq. (2.1), 



j P^expji j d!^xq^{x) i-f°Do + i^-D - m ^qp{x)^ 

X exp |i j (fx ^a{^)qa{x) + qa{x)Xa{x)] + iSyM 



(5.1) 



with the temporal and spatial components of the covariant derivative (in the fundamental 
color representation) given by Eq. (1.23). The Yang-Mills contribution to the generating 
functional is in the standard, second order formalism [27, 28], where the auxiliary fields 
TT and (j) do not appear. 

Now consider the following decomposition of the quark and antiquark fields: 

g„(x) = e-'^-^[h{x) + H{x)]^, 

qjx) = e'^^'^[h{x)+H{x)]^, (5.2) 



such that 



K{x) = e'--o [P+q{x)]^ , Ha{x) = e''"^" [P-q{x)]^ 
K{x) = e-'--° [q{x)P+U > Ha{x) = e"''"^" [q{x)P-U ■ (5-3) 



In the above, the (spinor) projection operators are given by^ 

P± = i(l±70), P+ + P_ = l, P+P-=0, Pl = P±. (5.4) 

This decomposition is a particular case of the heavy quark transform underlying HQET 
[81]. There, the starting point is the observation that a heavy quark within a hadron is 
almost on-shell and moves with the hadron velocity v. Its 4- momentum can be written 

p/^ = mv^" + k^" (5.5) 

where |A:| ^ m\v\ and v"^ = 1 (such that when = 0, = m?). One then uses the 
general projectors P± = (1 it f)/2 and the exponential terms are generalized to e^^™^'^. 
The case used here corresponds to the rest frame of the quark, v^^ = (1,0). Intuitively, this 
corresponds to a shift in the position of the "zero energy" level, such that the energy of the 



^Recall that in the free Dirac theory the operators P± project on the positive and negative energy 
eigenstates (they select the upper and lower components of the quark spinor). 
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state m (the heavy quark at rest) is the new zero level. In other words, instead of the true 
energy pQ = m+ko of the heavy quark, the residual energy /cq is used. However, within the 
context of the generating functional Eq. (5.1), the decomposition Eq. (5.2) can be initially 
regarded simply as and arbitrary decomposition, which will later on prove to be useful in 
Coulomb gauge. In fact, this choice will result in an important simplification: the spatial 
components of the Yang-Mills Green's functions are absent at leading order in the mass 
expansion and they apear only in next-to-leading order. The projection operators satisfy 
the following further relations 







such that the following relations hold for the components of the quark field: 



hj^h = hh, H-f^H 



HH, h-f^H = H-f'^h = h-f'h = HfH = 0. 



(5.6) 



(5.7) 



Inserting the decomposition of the quark fields given by Eq. (5.2) into the generating 
functional Eq. (5.1) and using these relationships one obtains 



z[x,x] 



D$exp 



d X 



ha{x) [i-Do]„/3 h(s{x) 



+ ha{x) 



al3 



Hp{x) + H^{x) 



al3 



hp{x) + H^{x) [-2m - iL>o]„^ Hp{ 



exp i / d^x [e-'^'^^Xaix) [h{x) + H{x)]^ + ^^-<^ [h{x) + H{x)\^Xa{x)] +iSYM 



(5.8) 

At this point, it is illuminating to discuss the difference between the mass expansion (as 
used here) and HQET. As already emphasized in the introduction of this chapter, in this 
approach we work with the full quark fields, i.e. we retain the source terms for the quark 
fields g, q in order to derive the full gap, Bethe-Salpeter and Faddeev equations from 
x] (see below). Hence, the source term expression in the above is slightly modified 
by the appearance of the exponential factors. Since the Jacobian of the transformation is 
field independent (and thus trivial), our generating functional has not been altered, but 
merely rewritten in terms of different integration variables. This is in contrast to HQET 
where the quark sources are replaced with sources for the projected fields h and H. In 
principle, differentiation with respect to the sources of the generating functional Eq. (5.8) 
would then lead to the Green's functions of the projections /i, H of the spinor q. In order 
to derive a Green's function with an external heavy quark, the components H (multiplied 
with twice the quark mass) are integrated out (as below) and in addition their sources 
are set to zero [82]. 

Noticing that for the /i-field ('large') components, the quark mass parameter m does 
not appear directly, we integrate out the -ff-fields and get the following expression 



'[x,x] = j' 



P$Det [iDo + 2m] exp <^ i / d'^x K{x)[\DQ\^nhp{x) 



+ 



h{x)i^-D + e-'"'''''xix) [iDo + 2m]~^ i^-Dhix) + e'"'''°xix 



X exp i i / d^x [e-'"^^«x„(x)/i,(x) + e'"'^" ha{x)xa{x)] + iSyM 



(5.9) 
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Obviously, since we have integrated out a nontrivial component of the original quark 
field, our expression is nonlocal and this is where the heavy mass expansion is necessary. 
The determinant arising from the fermion functional integration depends on the temporal 
component of the gauge field and it can be written in the following way: 



Det [iDq + 2m] 



exp In [iDq + 2m] = N exp Tr In 
- TV + O (l/m^) 



1 + ^TV 



1 



ido + 2m 



(5.10) 



where is a (noninteracting) constant. If we expand the logarithm in power series we 
see that the first term, of 0{l/m), vanishes due to the fact that the generators T"" are 
traceless matrices, and the term at 0(l/m^) is proportional to Tr{T°'T^) = Nf (but, 
since we restrict to leading order in the mass expansion, this term will not appear in our 
calculations). 

Leaving the e^™^" factors as they are, we can thus write 



Z[X,X] = / 2?^exp <^ i / d^x [ K{x) [iDo]^^ hp{x) 



+ 



2m 



xexpji J d^x [e-^"^^OxJx)/i„(x) + e''"^°/i„(x)xa(x)] + iSyM^ + O {l/m''). (5. 11) 

Our generating functional is now local in the fields and arranged in an expansion in the 
parameter 1 /m (this will be referred to as the mass expansion although strictly speaking, 
it is an expansion in the inverse mass). However, locality does not mean that the above 
expression can be directly applied. Let us consider the classical (full) quark field in the 
presence of sources: 



I j V<^q^ix)exp{iS} = I 



SZ 



1 

Z 



-imxo 



haix) + 



2m 



ij-Dhix) + e'"''">x{x) }>exp{icS} + (l/m^) . 

(5.12) 



One sees immediately that even at 0{l/m), the classical quark field has components that 
involve interaction type terms {Dh ~ Ah) brought about by the truncation of the nonloc- 
ality. Of course, this is nothing more than the statement that the /i-field is nontrivially 
(and dynamically) related to the full g-field. It also means that if we want to use the 
nonperturbative gap and Bethe-Salpeter equations (i.e., those equations derived from the 
action for the full quark fields q and their sources x s-nd which are the equations of QCD 
as opposed to HQET) then we cannot expect that the mass expansion can realistically 
be extended far beyond the leading order in order to do practical calculations. As stated 
previously though, the aim is to investigate the connection between the Yang-Mills sector 
and the physical world made of quarks; the string tension that represents our goal is not 
dependent on the quark mass (both light and heavy quarks are confined in the same way, 
as far as we know). Therefore, we restrict our attention to the leading order in the mass 
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expansion as follows (and writing Dq explicitly): 

Z[x,x] = j P$exp|i j d^xK{x)[ido^+ gT^a''{x)\^php{x) 

X exp |i j (fx [e-''""«x„(x)/i,(x) + e"^^'K{x)Xc.{x)\ + i5yA/| + O (1/m) . 

(5.13) 

The standard machinery of functional methods is now employed. Just as for the usual 
QCD gap equation derived in Chapter 2, we start with the observation that up to boundary 
terms (which are assumed to vanish) the integral of a total derivative vanishes, and obtain 
for the quark field equation of motion: 

= / P$^^exp{i5} 
J dma[x) 

= j P${[i9o.+5rV(x)]„^/i/3(2;) + e^"^"°Xa(^)}exp{i5} + 0(l/m).(5.14) 

The field equation of motion for the antiquark gives equivalent results and can be neg- 
lected. From the generating functional of connected Green's functions VK[x, xli Eq. (2.7), 
and using the bracket notation, Eq. (2.9), for the derivatives of W with respect to sources, 
we obtain the classical quark fields, Eq. (2.14). 

Using the equation Eq. (2.11) and noticing that from Eq. (5.13), the following relations 
hold: 

^^^^ = [ P$e-''"^o/i„(x) exp {icS} + O (1/m) , 

diXaix) J 

{^^^l^} = - /"p$ei™^°7Ia(x)exp{i5} + 0(l/m) (5.15) 

(recalling the earlier discussion of neglecting the 0{l/m) terms), the field equation of 
motion, Eq. (5.14), can be written in terms of derivatives of W: 

= [ido^U <iX/3(x)> + [gTXp e^™"° [<ip''{x)ixp{x)> + <ip''{x)><ixp{x)>] 
+e'"'''°Xa{x) + 0{l/m). (5.16) 

Factoring out the exponential terms gives then 

= [i5ox - m]„^ <ix^{x)> + [gTXis [<^p''{x)iX^{x)> + <ip^{x)><iXf}{x)>] 

+Xa{x) + 0{l/m). (5.17) 

To continue, we employ the Legendre transform Eq. (2.15) in order to construct the 
effective action for the full quark fields (with the notations and conventions introduced in 
Chapter 2). The field equation of motion can then be rewritten in terms of derivatives of 
the effective action (which are the proper Green's functions when sources are set to zero): 

<iq^{x)>= [ido^ - m]^^qp{x) + [gTX^ [<ip%x)ix^{x)> +a^{x)qp{x)] + O (1/m) . 

(5.18) 
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With the field equation of motion written in the above forms, we can now derive the 
Feynman rules for the quark components of the theory (the Yang-Mills parts are already 
known [28]). Following the derivation presented in Chapter 3, we functionally differen- 
tiate Eq. (5.17), ignoring the interaction terms, and we get the tree-level propagator in 
configuration space 

= [ido, - m]^p <ix^(z)ix^(x)>(0) -iS^o.6{z - x) + O (1/m) . (5.19) 

Naively, one would write the solutions as 

<l(^) = TI^%:T + 0(l/m), W^^Jk) = -^ + Oil/m). (5.20) 



These propagators will turn to have a couple of striking features, which will be extensively 
discussed after deriving the nonperturbative solution for the gap equation. At this stage, 
we only emphasize that the quark and antiquark propagators must be treated separately, 
and this will turn to have important consequences for the bound state equations. Due to 
the mass expansion, we only have a single pole in the complex feo-plane, as opposed to the 
conventional quark propagator, which possesses a pair of simple poles. Hence, in order to 
define the Fourier transform, one must first define the Feynman prescription for handling 
the poles in the energy integral. For the quark propagator, we write 

Let us state for the moment that for the antiquark propagator, the following Feynman pre- 
scription is assigned (this will be explained in the context of the Bethe-Salpeter equation 
for bound states, see below): 



For the proper two-point function, we use functional derivatives of Eq. (5.18) and we 
get directly in momentum space 

= i[ko-m]5^^ + 0{l/m), 
rSl/3(^) = ■^[ko + m]6^p + 0{l/m), (5.23) 

Since the two-point function requires no Feynman prescription and is diagonal in the outer 
product of the fundamental color, flavor and spinor spaces, we have the following relation 
between the quark and antiquark two-point functions: 

In addition, WgqVgg = 1 as usual. For the three-point functions we obtain 



qqcraP ^ 



^TaB^kiMM) = -[gT\^ + 0{l/m). (5.25) 
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Notice the ordering of the indices for the Tgqo- vertex. Importantly, the tree-level spatial 
quark-gluon vertex does not appear at leading order in the mass expansion: 



Let us again stress that, because of our insistence of using the full quark sources, all 
the nonperturbative equations involving the quarks (the Dyson-Schwinger equations for 
two-point and three-point functions, the Slavnov-Taylor identities, the Bethe-Salpeter 
equation and the Faddeev equation) will not alter their form at leading order in the mass 
expansion, the only alterations being at the level of the tree-level factors. 

5.2. Truncation scheme 

In order to solve the nonperturbative system we must further specify our truncation 
scheme. In the context of the heavy mass expansion, we propose to consider only the 
dressed two-point functions of the Yang-Mills sector (i.e., the nonperturbative gluon 
propagators) and to set all the pure Yang-Mills vertices and higher n-point functions 
occurring in the quark equations to zero. 

Since the tree-level spatial quark-gluon vertex is suppressed by the mass expansion, 
in our approximation any diagram containing this vertex will not contribute at one- 
loop order in perturbation theory. Further, the fully temporal gluon Green's functions 
Fo-cro- 5 Fo-o-o-cr 5 • • • are zero at tree-level [28], and this means that the leading order perturb- 
ative corrections containing this vertices again vanish. This implies that that the number 
of loop diagrams arising because of the Yang-Mills vertices is heavily restricted and they 
first contribute to the next to leading order in perturbation theory. Physically, the most 
important point that will emerge is that when we set the Yang-Mills vertices to zero, we 
exclude the non-Abelian part of the charge screening mechanism of the quark color charge 
and any potential glueball states. On the other hand, the charge screening mechanism 
and glueball contributions of the gluon field (i.e., the color string) is implicitly encoded 
in the nonperturbative form of the temporal gluon propagator. We will come back to the 
physical implications of this truncation after deriving the confining potential between a 
quark and an antiquark from the Bethe-Salpeter equation. 

With the truncation scheme as outlined above, the Yang-Mills sector reduces to the 
inclusion of a single object: the temporal gluon propagator which is written as [28] 



There are three important features to this propagator. Firstly, there are indications on 
the lattice that the dressing function L'o-o- is largely independent of energy [49], justifying 
the energy independence of the above form. As we already emphasized in the first chapter, 
within the first order formalism one can justify that the temporal gluon propagator must 
have some part that is constant in the energy in order to cancel closed ghost loops and 
resolve the Coulomb gauge energy divergences'^. Second, the lattice analysis indicates 
that the dressing function Dg-o- is infrared divergent and is likely to behave as 1/}? for 
vanishing k"^. Since we are interested mainly in the relationship between D^^ (as the input 

''See also the more formal considerations of Ref. [34]. 




(5.26) 




(5.27) 
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of the Yang-Mills sector) and the string tension, we will not need the specific form until 
towards the end. Third, the product g^Da-a is a renormalization group invariant quantity 
and thus a good candidate for being related to the physical string tension [12, 27]. 



5.3. Gap equation 

5.3.1. Nonperturbative treatment 

Let us begin by considering the Dyson-Schwinger equation for the quark two-point proper 
function. In full, second order formalism, QCD it is given by^: 

{k, -u, u - k)W^g^^{u;)r'^^^^sji^, -k, k - uj)WtA,^ik - uj)} 

(5.28) 



, p(0)a 
qqAajSi 



(the spatial gluon propagator Waa will be unimportant here because of the suppression 
Eq. (5.26) of the spatial quark-gluon vertex). The gap equation is supplemented by the 
Slavnov-Taylor identity, Eq. (4.12), derived in Chapter 4 . 

In order to use the Slavnov-Taylor identity as input for solving the gap equation, we 
first apply our truncation scheme in the context of the heavy mass expansion at leading 
order. Starting with the dressed spatial quark-gluon vertex, consider the terms that 
contribute to the Dyson-Schwinger equation shown schematically in Fig. 2.3. According 
to the truncation scheme, we set all Yang-Mills vertices to zero, meaning that diagrams 
(c) and (e-i) are excluded. This then leaves us with the tree-level term (a) and the quark 
contributions (b), (d). However, all of these involve at least one tree-level spatial quark- 
gluon vertex, which is not present at leading order in the mass expansion. Thus, we obtain 
the nonperturbative result that 

^,Aap^{kl^k2M^O{l/m). (5.29) 

Similarly, the ghost-quark kernels of the Slavnov-Taylor identity, given their definition, 
Eq. (4.13), involve Yang-Mills vertices and can be neglected. Thus, in our truncation 
scheme and at leading order in the mass expansion, the Slavnov-Taylor identity reduces 
to 

kmq.ap{ki,k2, ks) = Tgg^siki) [i^T^l - [i^T^l Tggspi-ki) + O (l/m) . (5.30) 

L J op L J ao 

Clearly, the Slavnov-Taylor identity under truncations results in an Abelian type Ward 
identity. Moreover, since the temporal quark-gluon vertex is simply multiplied by the 
temporal gluon energy (an essential feature of Coulomb gauge Slavnov-Taylor identities, 
as opposed to covariant gauges) and the quark proper two-point function is color diagonal, 
we can immediately write the solution: 

^qaa^ih, k^M = ^ [Tqq{kl) " V-^.i-k^)]] + O (l/m) . (5.31) 



^This can be directly inferred from the expression Eq. (2.34), in the first order formalism, by setting the 
auxiliary fields n, to zero. 
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The above solution is trivially satisfied at tree-level. The apparent singularity arising 
for vanishing gluon energy (^3 = 0, but 7^ 0) must somehow be canceled by the 
difference of proper quark two-point functions. Since the spatial momentum configuration 
is arbitrary, this leads to the requirement that Tqq{k) — )■ T-qq^ko) + O (1/m). Later on we 
shall see that our solution does indeed satisfy this condition. The demand that the 
nonperturbative vertex solution to the Coulomb gauge Slavnov-Taylor identity be free of 
kinematic divergences (here, simply the factor) is a variation of the familiar covariant 
gauge situation considered in Ref. [88]. 

Inserting the results, Eq. (5.31) and Eq. (5.29), for the vertices, using the Feynman rules 
given by Eq. (5.23), Eq. (5.25), and the temporal gluon propagator given by Eq. (5.27) 
and resolving the color structure, the nonperturbative gap equation, Eq. (5.28), under 
truncation and at leading order in the mass expansion thus reads 

^qqaSih) = 1 [^0 " m] 6^5 " O^Cp [ '^^^'"'^^ ^Wqqo,(i{uJo) [Tqq{uJo) - Tqq{kQ)\i^^ 

J [ko - uJo){k - uiy 
+0{l/m). (5.32) 

There exists one particularly simple solution to this equation, given by 

^qqaisik) = iSa/B [ko - m - Ir] + O (1/m) (5.33) 



with the constant \duj = cfiuj/ (2Tr] 



31 



Ir = \9'C, [ ^^2^ + o (1/m) . (5.35) 

A short comment regarding the ordering of the limits in the spatial and temporal integ- 
rals and the potential divergences in the constant is in order. Here, we state that 
the temporal integral is performed first under the condition that the spatial integral is 
somehow regularized, i.e. finite (the implicit regularization is signaled by the subscript 
"r")^. Extracting the expression for the constant X^. from the gap equation under trun- 
cation, Eq. (5.32), we insert the solution given by Eq. (5.33), Eq. (5.34) (assuming that 
the spatial integral is regularized), and obtain 



{k-Cdy 27r R^oc J _ji [uo - m - Ir + ie] 
= / ^i^^^^i^ + O (1/m) . (5.36) 

Clearly, the effect of performing the temporal integral first is that the regularized constant 
Xr in the denominator factor can be shifted away and hence becomes irrelevant. Further, 
by shifting the integration variable in the spatial integral we find that this integral involves 
no external scale, and thus we arrive at our above result, Eq. (5.35). 



^This regularization will be assumed throughout the rest of this work. 
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Inserting the solution Eq. (5.33) into the Slavnov-Taylor identity, we also have that for 
the vertex 



and this means that the dressed temporal quark-gluon vertex is trivial and the gap equa- 
tion reduces to the rainbow truncation. 

Let us now briefly discuss the physical interpretation of these results. Firstly, it might 
be the case that there exist other solutions to the truncated gap equation. However, as 
will be shown in the next section, the above solution can also be derived from a semi- 
perturbative type of expansion. One has to bear in mind that in principle the fully 
nonperturbative solution might not be the same as the resummed perturbative solution. 

Secondly, as already outlined when defining the Feynman prescription, the quark propag- 
ator has a single pole, so cannot represent physical propagation (which requires a covariant 
double pole) and this arises obviously from the truncation of the mass expansion. From 
Eq. (5.34) it then follows that the closed quark loops (virtual quark-antiquark pairs) van- 
ish due to the energy integration, which implies that the theory is quenched in the heavy 
mass limit (see also Ref. [81] for an alternate discussion on this topic): 



The above result can be interpreted as being a consequence of the breaking of time reversal 
symmetry of the full Dirac equation (quark and antiquark moving either forwards or 
backwards in time according to causality), as a result of the mass expansion. Time 
reversal symmetry breaking means that there is only the quark (or only the antiquark) 
moving forward in time, corresponding to the above Feynman prescription, and, in this 
sense, the corresponding antiquark (or quark) is not present. On the other hand, the 
closed quark loop involves a quark going backwards in time (similarly for an antiquark 
loop) , which is prohibited according to our argumentation, so that such closed quark loop 
integrals vanish at leading order. Clearly, the breaking of the time reversal symmetry is 
also reflected in the presence of only a single pole in the heavy quark propagator. 

Thirdly, the quark propagator Eq. (5.34) is diagonal in the outer product of the funda- 
mental color, flavor and spinor spaces as a consequence of the mass expansion - physically 
this corresponds to the decoupling of the spin from the heavy quark system. In fact, W^'^^ 
is identical to the heavy quark tree-level propagator [81] up to the appearance of the mass 
term, and this is due to the fact that in HQET one uses the sources for the large /i-fields 
directly, while we retain the sources of the full quark fields. Also, note that the kinetic 
term of the (tree- level) propagator would read /2m in the denominator factor and 
hence appears at higher order in the mass expansion. Such terms are obviously important 
to the UV properties of the loop integrals but do not play any role in the infrared limit 
considered here. 

Finally, the fact that the solution involves potentially divergent constants is not a 
comfortable situation but does not necessarily contradict the physics, since the position 
of the pole has no physical meaning (the quark can never be on-shell). The poles in the 
quark propagator are situated at infinity (thanks to as the regularization is removed) 
meaning that either one requires infinite energy to create a quark from the vacuum or, 
if a a hadronic system is considered, only the relative energy is important. Indeed, it 
was shown some time ago [89] that the divergence of the absolute energy has no physical 




(5.37) 




(5.38) 
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meaning and only the relative energy (derived from the Bethe-Salpeter equation, see 
below) must be considered. Further, the divergences appearing in the quark propagator 
have no interpretation with regards to renormalization, at least within the context of the 
mass expansion to leading order. The mass parameter cannot be renormalized simply 
because one cannot construct an appropriate counterterm in the action. Also, the quark 
field renormalization is trivial at leading order, as one sees from the explicit form of the 
temporal quark-gluon vertex, Eq. (5.37). 

Having discussed the quark propagator, let us now discuss the antiquark propagator. 
Recall that at tree-level, we used a different Feynman prescription for the two denominator 
factors and this gives rise to some rather interesting physical consequences. As previously 
discussed, the heavy mass expansion employed here breaks the charge conjugation sym- 
metry relating particle and antiparticle, so we cannot expect that the two propagators are 
necessarily equivalent. Starting with the gap equation for full QCD, Eq. (5.28), we reverse 
the ordering of the quark and antiquark functional derivatives that form the quark Green's 
functions (still within the context of the full quark fields and sources) and rearrange the 
ordering to get the gap equation for the antiquark propagator: 

-r,,,.(-^^) = -rJ°L(-/c) 

du; {r'^g,s,i-k, u;, k - a;)H^,5^^(-a;)r(J;^ J-a;, k,u;- k)Wf„{k - u) 
+^gUS,j(-k, ^, k - uj)W,-,,p{-u)T^;^lp^^{-u:, k,u - k)WX%{k -co)}. (5.39) 
Applying our truncation scheme reduces the above to 

+ J duTl-^^s,i-Ku;,k-u)W,-,,p{-u;)r^^^^^i^^ 

(5.40) 

In similar fashion, we have the Slavnov-Taylor identity for the antiquark-gluon vertex: 



klr<'g^,p^{k2,ki,k3) = +rg-gps{k2) igT\ - igT" rg-gSa{-ki) + Oil/m). (5.41) 



nd 



T 



T 
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Before we give the solution to Eq. (5.40), let us remark that in order to construct a quark- 
antiquark pair in the Bethe-Salpeter equation (which has a physical interpretation of a 
bound state equation), one is not considering a virtual quark-antiquark pair but rather 
a system composed of two separate unphysical particles.'' Since the quark and antiquark 
lines of the Bethe-Salpeter equation are never connected by a primitive vertex (unlike 
the closed quark loop discussed above), we are allowed to assign the following Feynman 
prescription for the antiquark propagator: 

= [to+„""t + ..| ^°"-'""' 
r,5«M*:) = iSc0\h, + m-I,]+O{l/,n) (5.43) 



^It is known that in Coulomb gauge, Gauss' law forbids the creation of a colored state in isolation (the 
total color charge is conserved and vanishing [90]) and so, the Feynman prescription for the quark (or 
the antiquark) propagator has no physical meaning in isolation. 
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(a) (b) 



Figure 5.1.: Diagrams that contribute under our truncation scheme to the Dyson-Schwinger equa- 
tion for the temporal quark-gluon vertex (without prefactors or signs). Internal 
propagators are fuhy dressed and blobs represent dressed proper vertex and (re- 
ducible) kernels. Internal propagators represented by solid lines represent the quark 
propagator. 



with the corresponding solution for the vertex 



r^M(^i'^2,A;3) = - gT'^ ^ +0{l/m). (5.44) 

L J pa 

In the above, notice that the sign of the loop correction has changed and this will result (in 
the context of the bound state studies for mesons and diquarks) in a physical interpretation 
for the quark-antiquark Bethe-Salpeter equation as a whole. 

5.3.2. Semiperturbative treatment 

After solving the gap (and anti-gap) equation, we have seen that the solutions for the 
proper two-point function leads to a (nonperturbatively) bare temporal quark-gluon (and 
antiquark-gluon) vertex. In this section, we will reconsider this vertex within the context 
of a semiperturbative analysis. This will then introduce a technical feature crucial for 
considering the Bethe-Salpeter equation nonperturbatively. 

Under our truncation scheme, the nonperturbative Dyson-Schwinger equation for the 
temporal quark-gluon vertex involves the diagrams shown in Fig. 5.1. Semiperturbative 
expansion means in this case that in the loop expansion all internal propagators are taken 
to be dressed, but all internal vertices are tree-level. Clearly, the goal is to show that all 
the loop corrections (contained in diagram (b) of Fig. 5.1) vanish. For this, it suffices to 
consider two types of diagram, given in Figs. 5.2 and 5.3. 

The easiest to start with is Fig. 5.2, where a single ladder exchange correction to the 
temporal quark-gluon vertex is considered. This diagram (neglecting the overall color and 
prefactors) gives rise to the following scalar integral: 

dujD„„{k-Cj) 1 [°° duo ^^^^^ 



(k - ujY 27r y„oo [ooo-m-Zr + ie] [uq + qq - m - 1,. + le] ' 
Now we apply the following identity (for finite, real a, 6; the case a = h \s trivial): 

' ^°°d4, ' I .,Uo. (5.46) 



[z — a + ie\[z — b + le] {a — h) J_^ \[z — a + ie] [z — b + ie] 
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Figure 5.2.: Ladder type loop correction to the temporal quark-gluon vertex. Internal propagators 
are fully dressed: solid lines represent the quark propagator and springs denote the 
temporal gluon propagator. 




Figure 5.3.: Generic crossed box type loop correction to the temporal quark-gluon vertex. In- 
ternal propagators are fully dressed: solid lines represent the quark propagator and 
springs denote the temporal gluon propagator. The box represents any combination 
of interactions allowed under our truncation. 



Thus we see that where the spatial integral is regularized, the temporal integral vanishes. 
It is simple to see that the planar one-loop diagrams with two or more external temporal 
gluon legs (which under the truncation scheme considered here connect only to the internal 
quark line) and one internal temporal gluon will also vanish. 

Now let us consider a generic crossed box (nonplanar) type of diagram, illustrated in 
Fig. 5.3. Considering only the temporal double integral components of the explicit internal 
quark propagators, we have the following form 

^ diOQ dvQ 

[vo -ai + ie] [uq - a2 + ie] [uq + g'o - ^3 + ie] [^o - - Po - 04 + ie] 

dioo 



00 I'^o 

00 



a2 + ie] [ujQ + go - as + ie] ['^o - Po - - 04 + 2ie] 
1 1 



uio+Po + a.'i- ie] 



00 [[vo- ai + ie\ [Vo 

/oo 
,00 [wo - 02 + ie] [wo + go - as + ie] [wq - Po - ai - 04 + 2ie] 



= (5.47) 
where in the last line, we have used a variation of the identity Eq. (5.46). Thus we have 
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the result that the generic crossed box type of diagram shown in Fig. 5.3 also vanishes. 

With the result that both the single ladder type exchange diagram and the generic 
crossed box diagrams considered so far vanish, it is easy to see that any vertex dressing 
diagram will vanish (including all subdiagrams such as internal vertex corrections and 
so on), since all associated diagrams are merely variations or combinations of these two 
under our truncation scheme. We stress that this result is a consequence of the fact that 
the energy and Feynman prescription of the denominator factors follow the quark line 
through the diagram so that eventually the identity, Eq. (5.46), can be used. It is also 
precisely the reason why all closed quark loops vanish, according to the relation Eq. (5.38). 

Thus, the semiperturbative expansion confirms our previous result that the temporal 
quark-gluon vertex remains bare to all orders. Clearly, the result also applies to the 
antiquark-gluon vertex. With the corresponding simple forms for the self-energy integrals 
of the gap and anti-gap equations (which as we recall, reduce to the rainbow truncation), 
the results for the quark propagator functions are also confirmed. Notice though that 
whilst the all orders semiperturbative result must match the nonperturbative result, the 
converse is not necessarily true. It remains the case that there may exist further solutions, 
but these must be purely nonperturbative in character if they exist. 

To summarize, in this section we have introduced a set of valuable identities for the 
energy integrals, which confirm the nonperturbative case previously studied. This is very 
useful for further investigations, since we are allowed to apply them with confidence to 
study the bound state equations in the next section. 

5.4. Bound state equations 

Bound states of quarks appear as free-particle poles in their respective n-point Green's 
functions. For example, the quark-antiquark, and three-quark Green's functions exhibit 
meson and baryon poles, respectively. In this section, we will consider the correspond- 
ing bound state equations in the limit of the heavy quark mass - named Bethe-Salpeter 
equation for mesons and Faddeev equation for baryons - and interpret them in connection 
with the linearly rising potential which confines the quarks. As has been emphasized, due 
to the fact that we include the full quark sources and fields in our generating functional, 
we are allowed to use the full functional (nonperturbative) equations as a starting point 
and then subsequently apply our mass expansion and truncation scheme. In this con- 
text, we will also seek for a relationship between Yang-Mills sector of the theory and 
nonperturbative external physical scale (i.e., the string tension). 

5.4.1. Bethe-Salpeter equation for mesons and diquarks 

Let us now consider the homogeneous Bethe-Salpeter equation for quark-antiquark bound 
states. In full QCD this equation reads 

rap{p;P) = - I dkK^^.s^{p,k;P) [W^g{k+)r{k;P)W^,{k_)l^s (5-48) 

In the above, the minus sign arises from the definitions of the Legendre transform and 
Green's functions in Coulomb gauge. This will be shown in Chapter 6, where the 4- 
point Green's function (which includes the homogeneous Bethe-Salpeter equation) will be 
explicitly derived. The momenta of the quarks are given by k^ = k+£^P, A;_ = k + {S^ — l)P 
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Figure 5.4.: Homogeneous Betlie-Salpeter equation for quark-antiquark bound states. Internal 
propagators are fully dressed and solid lines represent the quark propagator. The 
box represents the Bethe-Salpeter kernel K and filled blobs represent the Bethe- 
Salpeter vertex function T with the (external) bound state leg given by a dashed line. 
See text for details. 



(similarly for p±), P is the pole 4-momentum of the bound state (assuming that a solution 
exists), ^ = [0, 1] is the so-called momentum sharing fraction that dictates how much of 
the total meson momentum is carried by each quark constituent. K represents the Bethe- 
Salpeter kernel and T is the Bethe-Salpeter vertex function for the particular bound state 
that one is considering and whose indices explicitly denote only its quark content. Later 
on, we will explicitly investigate what color, flavor and spin structure the solutions may 
have. We also mention that physically, the results should be independent of ^ and this has 
been numerically observed in phenomenological studies [91]. The Bethe-Salpeter equation 
is shown pictorially in Fig. 5.4. 

To solve the Bethe-Salpeter equation, we firstly address the problem of constructing 
the kernel K. For technical reasons the most widely studied system is based on the ladder 
kernel which is either constructed via the interchange of a single gluon (for example [92] ) 
or as a phenomenological potential (see for example Ref. [89]). However, there has been 
much recent attention focused on the construction of more sophisticated kernels. One 
key element of the construction is the axialvector Ward-Takahashi identity [AXWTI], 
which relates the gap equation to the Bethe-Salpeter kernel and which ensures that chiral 
symmetry and its spontaneous breaking are consistently implemented (e.g., Refs. [89, 93, 
94]). Here, we shall show that the ladder Bethe-Salpeter kernel is exact at leading order 
in the heavy mass expansion and under our truncation scheme. 

To construct the kernel we follow the semiperturbative analysis of the previous sec- 
tion. As before, at leading order in the mass expansion and under our truncation scheme 
we only have the temporal quark-gluon and antiquark-gluon vertices, both of which have 
been shown to be given by their tree-level forms. Consider now the generic semiperturbat- 
ive crossed box contribution, which contains all possible nontrivial contributions allowed 
within our truncation scheme. This diagram is depicted in Fig. 5.5. Such a diagram 
has at least the following terms in the temporal integral (as before, we assume that the 
spatial integral is regularized and finite so that we are able to firstly perform the temporal 
integral without complication) 

/ T -a P^y T- (5-49) 

J p\ - m -Zr + ie\ . . .[ujQ - k}t +m -Zr + ie\ 

The first factor corresponds to the explicit quark (upper) propagator, the last factor 
to the explicit antiquark (lower) propagator. The dots represent the multiple internal 
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-k. 

Figure 5.5.: Generic crossed box type of diagram that contributes to the Bethe-Salpeter ker- 
neL Internal propagators are fuUy dressed, whereas vertices are tree-leveL The 
upper (lower) solid line denotes the quark (antiquark) propagator; springs denote 
the temporal gluon propagator and the box represents any combination of nontrivial 
interactions allowed under our truncation scheme. See text for details. 




propagator factors which carry the same dependence on the integration energy and 
have the same relative sign for the Feynman prescription term, i.e., ujq + ie, regardless 
of whether they originate from internal quark or antiquark propagators. Therefore, this 
type of integral can always be reduced to the difference of integrals over a simple pole 
and with the same sign for carrying out the analytic integration, just as in Eq. (5.46). 
Thus, all generic crossed box diagrams in the Bethe-Salpeter kernel are zero and one is 
left with simply the ladder contribution to the kernel. In fact, this result can also be 
inferred from the AXWTI: this identity connects the self-energy term of the gap equation 
and the Bethe-Salpeter kernel and since it has been shown explicitly that the self-energy 
integral reduces to the rainbow truncation, the corresponding Bethe-Salpeter kernel is 
simply given by ladder exchange. 

Having analyzed the kernel K, let us now inspect the quark and antiquark pieces of 
the Bethe-Salpeter equation. Recalling that the antiquark propagator must be treated 
as distinct from the quark propagator, in Fig. 5.4 we must explicitly specify which line is 
assigned to the quark propagator and which one, to the antiquark propagator. For the 
quark-antiquark system under consideration (later on we will analyze the quark-quark, or 
diquark system), the Bethe-Salpeter equation more properly reads 

T^p{p;P) = - I dkK^p.Sj{p,k;P)[W^,{k+)T{k-P){-)Wj^{-k.)]^^ (5.50) 

where we have explicitly identified the antiquark propagator contribution (it corresponds 
to the lower line of Fig. 5.4) by reordering the functional derivatives, i.e. Wqq{k-) = 
—W^{—k-). The above equation is still valid in full QCD. Further, we explicitly replace 
the tree-level form for the kernel, and we arrive at the following expression Bethe-Salpeter 
equation for the quark-antiquark system, at leading order in the mass expansion and 
within our truncation scheme: 



rapip;P) = - j dk[V^^^{p+,-k+,k-p)W-qq{k+)] 



{-)W'^^{-k_){-)T%^{-p_,k.,p-k) Wf„{p-k)Ty^{k-P) + 0{l/m). 

(5.51) 
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Notice the antiquark contribution within the kernel which absorbs the expUcit minus sign 
from Eq. (5.50). Also, recall that we are implicitly considering only the flavor non-singlet 
case. The above equation can be trivially rewritten as: 



-k+,k- p)Wt{p - k)T%,si-P^,k^,P- k) 



xW^g,x{k+)W^gS.i-k-)rx.{k; P) + O (1/ 



m] 



(5.52) 



Inserting the nonperturbative results for the propagators and vertices so far, Eqs. (5.34), 
(5.43), (5.37), (5.44) and taking the form, Eq. (5.27), for the temporal gluon propagator, 
we obtain the equation 



dkD„„{p- k) 



{p- kf 



2tt 



dkn 



[kl 



m 



Cr + ie] [fcO - m + Xr- - ie] 



[r'^r(A;;P)r«^ + 0(l/m). (5.53) 



We see immediately that the flavor and spin structure of the meson decouples from the 
problem — this is well a known property of the heavy mass expansion. The color structure 
will be discussed shortly. Since the external energy po does not enter the right-hand side 
of the above equation, we further have that the Bethe-Salpeter vertex r(^^(p;P) must 
be independent of the quark energy (and only implicitly dependent on the bound state 
energy Pq). Thus we can write 



dkD„^{p- k) 



(p- ky 



-9 



T''T{k-P)T"' 
dko 



al3 



■ Ir + is] [fe- — m + Ir — ie] 



m 



+ 0{l/m) 



dkD„„{p- k) 



TT{k;P)T'' 



al3 



{p-kf 



[Pq - 2Zr + 2ie] 



+ 0{l/m) . 



(5.54) 



Thus, at leading order in the mass expansion, inserting the expression Eq. (5.35) for Ty. 
it is now clear that 



dCoDaa{uj) 



rai3{p;P) 



-g 



dkD^„{p- k) 



{P- kf 



T''T{k]P)T'' 



+0{l/'m) . 



a/3 

(5.55) 



We notice that the explicit quark mass contributions of the self-energy expressions cancel. 
This is a feature of the quark-antiquark Bethe-Salpeter equation — it does not make any 
reference to the origins of its constituents (this explains, for example, why the pion can 
be a massless bound state of massive constituents). Physically, one can visualize that the 
quark and antiquark are moving with equal and opposite velocities such that the center 
of mass is stationary. This is related explicitly to the choice of Feynman prescription for 
the constituent quark and antiquark. Were the Feynman prescription for the antiquark 
chosen to coincide with that of the quark, the right-hand side of Eq. (5.54) would simply 
vanish and there would be certainly no physical quark-antiquark state. The Feynman 
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prescription for the antiquark corresponds precisely to a particle moving with the opposite 
velocity. Also, at leading order in the mass expansion, the momentum sharing parameter, 
^, has dropped out, and therefore the results are independent of ^. Shifting the integration 
momenta, we can write 

PoTapip; P)=9^ j '^"'^r^'^^ [CpT^pip-, P) - [T'^Yip- Pm^p] + O (1/m) . 

(5.56) 

To see the physical meaning of this equation, we rewrite the Bethe-Salpeter vertex func- 
tion as a Fourier transform: 

V^p{p-P) = I dye~'P yr^^{y) (5.57) 

(in the homogeneous Bethe-Salpeter equation, the total momentum P denotes the solution 
and is not a variable). We also assign the following color structure 

[T^mTXp = CMTapiy) (5.58) 

where Cm is yet to be identified. Then, the Bethe-Salpeter equation reduces to 

J dye~'P-yPoT^p{y) = J dye-'^ yg^ J ^^^g-^^) {CpT^p{y) - e'^-yCMTaM} 

+0{l/m) (5.59) 

with the simple solution 

P^ = g2l^ d^^B. |c7^ _ e^^-yCM} + O (1/m) . (5.60) 

As already mentioned, because the total color charge of the system is conserved and 
vanishing [90], neither the quark or antiquark can exist as an independent asymptotic 
physical state. Thus, the qq system is either confined, such that the bound state energy 
Pq increases linearly as the separation between the quark and antiquark increases, or 
the system cannot be physically created, such that the energy Pq is infinite when the 
hypothetical regularization is removed. Whether the system is confining or disallowed 
can only depend on the color structure, since the temporal gluon propagator dressing 
function would be common to both situations. 

An infrared confining solution is characterized in configuration space by the solution 

Pq = for large (5.61) 

where a ~ IGeV/ fm is called the string tension, such that as the separation between the 
quark and antiquark increases, the energy of the system should increase linearly without 
bound and infinite energy input is required to fully separate them^ (at least in the absence 
of unquenching) . The Fourier transform integral needed to obtain the above form of the 
solution is 

If - ='^'1 - (-^) 

*The small \y\ (and large |cl5|) properties are of no concern here. 
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This implies that the temporal gluon propagator dressing function diverges like 1/uP and 
in addition the condition 

Cf = Cm (5.63) 

must be satisfied. Moreover, with the these conditions, the spatial integral in Eq. (5.60) 
becomes automatically convergent and hence the energy of the system is well-defined, as 
the regularization is removed (since we are interested in the low regime, it becomes 
clear that the regularization here would be infrared in character). Using the Fierz identity 
for the generators, Eq. (A. 12), we get the condition 

Cp^a-iiy) = CM^a^iv) = I3f}{y) - l^^a-^iv), (5-64) 

or with the definition Eq. (A.ll), 

^a-i{y) = 5a^T{y). (5.65) 

In other words, the quark-antiquark Bethe-Salpeter equation can only have a finite solu- 
tion for color singlet states where the divergent constant integral coming from the un- 
physical quark self-energy cancels; otherwise the energy of the system is divergent. 

Assuming that in the infrared (as is indicated by the lattice data [49] or by the above 
argument about the non-existence of asymptotic quark states), = X/u"^ where X is 
some combination of constants (and further knowing that g^X is a renormalization group 
invariant [12, 27]), then 

Po = a\y\ = + O (1/m) . (5.66) 

The above result shows that there exists a direct connection between the string tension 
and the nonperturbative Yang-Mills sector of QCD (i.e., the temporal gluon propagator) 
at least under the truncation scheme considered here. 

Let us now shortly discuss what are the possible consequences of including the pure 
Yang-Mills vertices in our approach. Recall that the confining potential follows from the 
rainbow-ladder Bethe-Salpeter kernel (i.e., dressed temporal gluon propagator and tree- 
level quark-gluon vertices). This suggests that the dressing function D^ja does implicitly 
contain all nonperturbative effects associated with the dynamical dressing of the color 
charge (including, for example, potential glueball states), whereas the quark-gluon vertices 
correspond to the naked quark color charge. Pictorially, one can visualize this as a dressed 
color string confining two naked color sources. Consequently, we anticipate that the effect 
of including the non-Abelian corrections to the formalism presented here would not result 
in the removal of the linearly rising bound state energy (this would correspond to the 
cancellation of the ladder exchange). Instead, one expects a lowering of the string tension 
a by shifting the pole position by some finite amount. Physically, this corresponds to the 
screening of the quark color charge — this statement has been phrased by Zwanziger as 
"no confinement without Coulomb confinement". 

Let us now turn to the diquark Bethe-Salpeter equation. From a technical point of 
view, the difference between this and the previously considered quark-antiquark system is 
rather simple, but as we shall see it leads to a completely different physical result. Since 
the quark and the antiquark propagators share the same Feynman prescription relative 
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to their energy, the result that the crossed box contributions to the Bethe-Salpeter kernel 
vanish extends to the diquark case. This means that we can immediately write down the 
Bethe-Salpeter equation for diquarks, at leading order in the mass expansion and within 
our truncation scheme: 

xW^g^x{k+)W^gS.{-k^)rx.{k;P) + O (1/m) . (5.67) 

Again, the indices of the Bethe-Salpeter vertex function correspond to the quark content 
of the diquark and since the flavor and spin content decouple from the system, we shall 
only be interested in the color content of the diquark. Expanding this out as before, we 
get the analogous result 

Po-2m-g Cf r„^(p; P) = g ^2 [T ] J Txk{p - P) 

+0{l/m). (5.68) 



Fourier transforming as previously, and writing 



[^I/Jk rA«(y) = CoTaM (5.69) 

gives the solution 

Po = 2m + g^J^ ^^^7^^^ [Cf + e'^'^^C^} + O (1/m) . (5.70) 



The dependence of the solution on the quark mass simply indicates that in contrast to 
the quark-antiquark system, there are now two co-moving quarks. For the anti-diquark 
system the solution is identical to the above, but with minus twice the mass - their 
velocities are simply reversed. The diquark is antisymmetric under interchange of the two 
quark legs and this means that the color structure must be antisymmetric. Similar to the 
quark-antiquark system, the system can only have a confining (finite) energy solution, 
i.e., if Cd = —Cp, or no finite solution at all. Using the Fierz identity Eq. (A. 12), the 
color condition then reads 

- CpTafjiy) = CoTafiiy) = ^Tfi^iy) - ^Tapiy). (5.71) 

Demanding the diquark color antisymmetry and with the definition Eq. (A. 11) this be- 
comes 

-Nc-2 = Nc = -1, 2. (5.72) 

This means that in SU{Nc = 2) there exists a confined, antisymmetric bound state of 
two quarks - the SU{2) baryon - and otherwise no physical states are allowed. In the 
next section, we will consider the SU{3) baryon, and explicitly demonstrate how the 
confinement potential is obtained in this case. 
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Figure 5.6.: Faddeev equation for three quark bound states. Solid lines represent the quark 
propagator, the box represents the diquark kernel K and the ellipse represents the 
Faddeev vertex function with the bound state leg depicted by a triple-line. See text 
for details. 

5.4.2. Faddeev equation for baryons 

Baryons appear as three-quark bound-states in the Faddeev equation, which generally can 
be written as: 

r = K^^^r (5.73) 
where T is the quark-baryon vertex and A'(^) is the three-body kernel, containing an 

(3) 

irreducible term K^^ and the sum of permuted diquark kernels: 

i=l 

The Faddeev equation [95] and its subsequent developments [96, 97] (for an extended 
review see [98] and the textbook [99]) provide a general formulation of the relativistic 
three-body problem. It is a bound state equation (the direct analogue of the homogen- 
eous two-body Bethe-Salpeter equation) and it has been efficiently applied in QCD to 
study baryon states, via the Green's functions of the theory. Typically, these studies are 
performed in Landau gauge and, due to the complexity of the equations, they have been 
mainly restricted to rainbow-ladder truncation, where the kernel is reduced to the single 
exchange of a dressed gluon. Within this approximation and by employing phenomeno- 
logical ansatze for the Yang-Mills part of the theory, the nucleon and A properties have 
been studied [100-103]. Other simplifications include the three-body spectator formal- 
ism [98], a Salpeter-type equation with instantaneous interaction [104] or the diquark 
correlations [105]. 

In the following, we investigate the Faddeev equation for three-quark bound states, 
by employing only the permuted two quark kernels K\ (which coincide with kernel 
appearing in the Bethe-Salpeter equation for diquark states) and neglecting the three- 
quark irreducible diagrams, i.e., genuine three-body forces [106]. This approximation is 
also motivated by the fact that in the quark-diquark model the binding energy is assumed 
to be mainly provided by the two-quark correlations [107] . In this truncation, the Faddeev 
equation reads (see also Fig. 5.6): 

ra/37(Pl,P2,P3;^) = 

- j dk [^K^^.^,p,{k)W-^^i^n{pi + k)W-^i^,pn{p2 - k)T^npn^{pi + k,p2 - k,p3;P) 
+K^p.p>^> {k^)W-^l^,p» {P2 + k)W-^^>y, {p3 - k)T^i^„y, {pi,p2 + k,p3- k; P) 

+ ^a7;ya'(^)^,-,yy'b3 + k)W^ga'a"iPl - k)r^^>>y^{pi - k,p2,P3 + k;P)] (5.75) 
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where Pi,P2,P3 are the momenta of the quarks, P = P1+P2+P3 is the pole 4-momentum 
of the bound baryon state and T is the so-cahed quark-baryon Faddeev vertex for the 
particular bound state under consideration and whose indices denote explicitly only its 
quark content. Due to the fact that in the heavy mass limit the spin degrees of freedom 
decouple from the system, at leading order in the mass expansion the Faddeev baryon 
amplitude Ta/s-y becomes a Dirac scalar, similar to the heavy quark propagator Eq. (5.34). 
The explicit momentum dependence of the kernels K is abbreviated for notational con- 
venience. As in the homogeneous Bethe-Salpeter equation, the integral equation depends 
only parametrically on the total four momentum P. 

As discussed in the case of the Bethe-Salpeter equation, the kernel K reduces to the 
ladder approximation (constructed via gluon exchange) and it reads 

K,. ,^,(A;) =r^ ,W^^(k)T''. = g^T'' ,WfJk)Tl„, (5.76) 

with the temporal quark-gluon vertex and the temporal gluon propagator given by Eq. (5.27) 
and Eq. (5.23), respectively. Similar to the Bethe-Salpeter equation for meson bound 
states, the energy independence of this propagator will turn to be crucial in the deriva- 
tion of the confining potential. 

Let us now investigate the energy dependence of the equation Eq. (5.75). As shown in 
the previous section, the Bethe-Salpeter kernel was energy independent, and thus it was 
straightforward to show that the Bethe-Salpeter vertex itself did not contain an energy 
dependent part. This observation was then used to calculate the confining potential from 
the Bethe-Salpeter equation, via a simple analytical integration over the relative energy 
variable. Unfortunately this approach cannot be extended to baryon states: despite the 
instantaneous kernel, a relative energy dependence still remains and thus one cannot 
assume an energy-independent Faddeev vertex. Therefore, in order to proceed, we make 
the following separable ansatz for the Faddeev vertex: 

ra/37(pi,P2,P3;^) = ^a/37rt(p?,P^,pll;P)r.(#i,p2,P3;^) (5.77) 

where we have introduced a purely antisymmetric (in the quark legs) color factor ^ (the 
possible baryon color index is omitted) and the symmetric (Dirac scalar) temporal and 
spatial components Fj and F^, respectively. 

Inserting the explicit form of the kernel Eq. (5.76) and the quark-baryon vertex ansatz 
Eq. (5.77), the Faddeev equation Eq. (5.75) can be explicitly written as (for simplicity we 
drop the label P in the arguments of the vertex functions): 

ra/3^(pi,P2,^53) = -g'^T^rT^^^TKj j dkW,,{k) 

X Wg5(pi + k)Wqq{p2 - k)Vt{p\ + ko,p9, - A;o,P3)F,(pi + k,p2- k,p^) + c.p., (5.78) 

where the explicit color structure has been extracted {W^^ = 5°'^Waa, Wqqap = SapWqq). 

Analogous to the Bethe-Salpeter equation, we use the Fierz identity for the generators, 
Eq. (A. 12), to write the color structure as 

with Cb = {Nc + 1)/2A'^C5 where is the number of colors, yet to be identified (i.e., the 
baryon is not assumed to be a color singlet). 
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In the next step we perform the Fourier transform for the spatial part of the equation, 
recalhng that the heavy quark propagator is only a function of energy. We define the 
coordinate space vertex function via its Fourier transform 

rs(pi,P2,P3) = J dfidf2rf^3e-^Pi-"^-^^^^-^^-^^^^-"^r,(fi,f2,X3) (5.80) 



such that 



J dkWaa{k)Tsipi + k,P2-k,p3) = 

J dxidx2dx3 e~'P'-^'~'P'-'^'~'P^-^-m^^{x2 - f i)r,(f 1, f2, xs). (5.81) 

Clearly, the component Tg trivially simplifies (as before, we have separated the temporal 
and spatial integrals, under the assumption the spatial integral is regularized and finite) 
and the equation Eq. (5.78) reduces to [dko = dkQ/{2TT)]: 



rt{p'i,p'i,p'i) = g^CBW^^{x2 - xi) j dko Wg.ip'i + ko)Wggipl - ko)Tt{p'i + A;o,p^ - ko,pl) 
+ c.p.. (5.82) 

At this point we make a further simplification, motivated by the symmetry of the three- 
quark system: we restrict to a particular geometry, namely to equal quark separations, 
i.e. |f| = 1^2 — a^il = 1 3^3 — 2^2! = |xi — X3I. By inserting the explicit form of the quark 
propagators, Eq. (5.34), we have 

^t{Pi,P2,P3) = -9 CBW^^{\i^) / dko ^ , • 1 r — 1 ^ , . 1 

J [p\ + ko - m -Ir + ie\ [p)^ - ko - m - Ir + ie\ 

+ c.p.. (5.83) 

With the assumption that the vertex Fj is symmetric under permutation of quark legs, 
an ansatz that satisfies this equation is: 

r.(AApg)-_£ ,^„_3(^„;„^,^)^,^ . (5.84) 

Since the explicit derivation is rather technical, we only give here the solution and defer 
the details to the Appendix E. 

Notice that in the expression Eq. (5.84) there are simple poles (in the energy) present. 
These poles however do not occur for finite energies and cannot be physical. As discussed, 
this is also the case for the quark propagator. Intuitively, when a single heavy quark is 
pulled apart from the system, the qqq state becomes equivalent (i.e., it has the same 
color quantum numbers) to the qq system in the sense that the remaining two quarks 
form a diquark which for A'^c = 3 would be a color antitriplet configuration, and hence 
the physical interpretation of the vertex Eq. (5.84) can be directly related to the heavy 
quark propagator Eq. (5.34): the presence of the single pole in Eq. (5.84) simply means 
that this cannot have the meaning of physical propagation (this would require a covariant 
double pole). Moreover, the divergent constant Ir appearing in the absolute energy does 
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not contradict the physics - the only relevant quantity is the relative energy of the three 
quark system. 

With this ansatz at hand, we return to the formula Eq. (5.83), insert the definitions 
Eq. (5.27) and Eq. (5.35) for W^a^x) and X^, and arrive at the following solution for the 
bound state energy Pq, in the case of equal quark separations: 



Fq = 3m + -g I auj 



Cf - 2CBe" 



(5.85) 



The following reasoning is similar to our discussion from the case of the bound states of a 
meson or diquark system. Since the quarks can not be prepared as isolated states, the only 
possibilities for the qqq state are either that the system is confined (i.e., the bound state 
energy Pq increases with the separation), or the system is physically not allowed (i.e., the 
energy Pq is infinite). From the formula Eq. (5.85) and knowing that Dcrcr(w) is infrared 
enhanced, it is clear that in order to have an infrared confining solution (corresponding 
to a convergent three- momentum integral), the condition 



(5.86) 



must be satisfied. This is fulfilled for Nc = 3 colors, implying that ^'^/j^ = Ea/j^ and that 
the baryon is a color singlet (confined) bound state of three quarks; otherwise, for A'^c 7^ 3 
the energy of the the system is infinite for any separation \f\. 

As in the preceding section, with the assumption that in the infrared Da-a{(^) = X/io'^ 
(as indicated by the lattice data [47-49, 108, 109] and by the variational calculations in 
the continuum [41]), where X is some combination of constants, it is straightforward to 
perform the integration on the right hand side of Eq. (5.85), with the result that for large 
separation \r\: 

Sg^CpX. 



Pn 



3m H — - 
2 



87r 



\r\. 



(5.87) 



This mimics the previous findings for qq and qq systems, namely that there exists a 
direct connection between the string tension and the nonperturbative Yang-Mills Green's 
functions (at least under truncation). In this case, the standard term "string tension" 
refers to the coefficient of the three-body linear confinement term cT3q|r|. Also, comparing 
with the result of the previous section, we find that the string tension corresponding to the 
qqq system is 3/2 times that of the qq. To our knowledge, no direct comparison between 
the string tensions of the two systems has been made and hence this relation would be 
interesting to investigate on the lattice. The appearance of three times the quark mass 
stems from the presence of the mass term in the heavy quark propagator Eq. (5.34) which 
enters the Faddeev equation. 
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We have seen that in general, the underlying equation for the description of meson bound 
states is the two-body homogeneous Bethe-Salpeter equation. In the rainbow-ladder 
approximation, this equation has been successfully used to describe the properties of 
light mesons (see, for example [91, 92] and for a recent review [43]), where the driving 
mechanism is the chiral symmetry breaking. Beyond this approximation, models with 
dressed vertex contributions [93, 94, 110-113] and unquenching effects [114-116] have been 
considered, and more sophisticated numerical methods to solve both the homogeneous and 
the inhomogeneous Bethe-Salpeter equation have been recently developed [117]. In spite 
of this success, an exact derivation of the meson or diquark bound state energies via 
Green's functions techniques has not been yet reported. The difficulty stems from the 
fact that the (irreducible) interaction kernel contains higher order vertex functions which 
in general can not be calculated exactly. 

In this chapter, we continue our investigations of the quark-antiquark and diquark 
states by using Green's functions techniques. As before, our study is based on the heavy 
mass expansion underlining HQET and with the truncation of the Yang-Mills sector to 
include only dressed two-point functions. By means of functional methods, we explicitly 
derive the (fully amputated) quark 4-point Green's functions and give an exact, analytical 
solution. This will enable us to verify that bound states are related to the occurrence of 
the poles in the Green's functions and hence we will be able to provide a direct connection 
between the homogeneous Bethe-Salpeter equation considered in the previous chapter and 
the singularities of the Green's function, at least within the scheme considered here. 



6.1. 4-point Green's functions for quark-antiquark systems 

In this section we derive the Dyson-Schwinger equations for the one-particle irreducible 
and for the fully amputated 4-point quark-antiquark Green's functions (from which later 
on the bound state energy for the heavy quark system will be extracted). The presentation 
follows [118]. 

Let us start by deriving the Dyson-Schwinger equation for the one-particle irreducible 
(IPI) 4-point function. As an illustration of the functional differentiation techniques, we 
present the explicit derivation of the first term in this function and notice that the rest 
of the terms follow from an identical calculation. Just as for the gap equation, we first 
take the functional derivative of the quark Eq. (2.19) with respect to iQ'yz, and obtain the 
result Eq. (2.29); further, we functionally differentiate with respect to iqrw and, using the 
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(f) 



(g) 



(h) 



(i) 



(j) 



Figure 6.1.: Diagrammatic representation of the one particle irreducible 4-point quark-antiquark 
Green's fmiction. Blobs represent dressed proper (IPI) 4-point vertex, solid lines rep- 
resent the quark propagator, springs denote cither spatial (A) or temporal (cr) gluon 
propagator and cross denotes the tree level quark-gluon vertex. Internal propagators 
arc fully dressed. 



product rule, we obtain: 

<iPj/lX/3x- 



-S[X,n] 



<iPyUx> 



<i$Aig72i^'K><i^KiX/3x> 



+ <iplUx><i^xiq-yzi^K> 



6iqr 



(6.1) 



As explained above, we only retain the first term in the product (and denote the rest with 
dots). Again, we make use of the formula Eq. (2.27) and obtain: 



<iPyiX/3x> = S[X,K,fl,u] 



< <ippJf,><i<^f,iqrwi'^u><iJu^Jx><i'^xiqjzi^K><iJ^iXi3x> +■■■■ (6.2) 
A last functional derivative with respect to the quark field qr^t gives 

J3 



6iqrwSiq^zSiqr,t 
5 



<iPyiX/3x>= S[X, K, p, 



<i<I>^,ig™i$^> <i J^i Ja> <i<^xiqyzi'^K> <i J«iX/3x> 



+ <\pylJf,> <\^f,iqr^Mrw'^^v> <U^Ux> <i^xiq-fzi'^K> <iJKiX(Sx> 



+ <i/9^iJ^><i«>p,i,?™i$^> 



Siqrit 



<\J?;Ux> 



<\^xSq-iz'^^K><^JKi-xpx> 



+ <\py\J^> <\^^iqrwi'^u> <i J;ji JA><i^ AlQr^t ig72i$K> <^Jk^XI3x> 
+ <iPyU^,> <i^^i^™i^z^> <i J"i JA><i^ Aig72i^K> 



Siqrjt 



<i-^KiX/3x> 



■+....(6.3) 
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Again, we take only the first term from the above sum and as before we use the formula 
Eq. (2.27). We obtain: 



5^ 



diqrw^iqjjiqrit 



<i/5^iX/3x> 



y^A.px^ - -S[X,K,fi,u,e,6] <ipyiJs><i^ekvt^^5><^Js^Jfi> 
X <i^>^i^™i$j.><iJ^iJA><i$Ai97zi^«><iJKiX/3x> + • • • • (6.4) 



Identifying the various fields and sources, we arrive at the the first term in the expres- 
sion below for the 4-point Green's function (all the other terms are derived by a similar 
calculation and we omit the terms which will vanish when the sources are set to zero): 



+ 



+ 
+ 



<iX/3iXK> <\q^iq^\qx^qr,> <ig^igi.icr|^> <iXi.iXA> <\p^^\Py> 
<iX;3iX5> <i^5i^r?io-> <ip£ipf> <^Mu'^(^li> <iX!.iXA><%i97io"f> <ip|^iPy > 
<iX/3iXK> <i9Kig7i^r i9A> <iXAiX5> <ig5ig,,icT^> <ip£ipy> 
<iX/3iXK> <\q^\q^\qr^qr,<y\> <ipiipy> 
<iX/3iX5> <iqsiqvK> <^Pe^pi> <i^r i97i^Ki^A> <^Px^Pt> 
<^Xl3^Xiy> <^uktMMr,> <ix^iXK> <i^Ki97i'^A> <W\^Py> 



<1X/31X5> <lQ5l9rylO"£> <^Pe'^Pu> <'^qT^qfi^Cru> <1X^1Xk> <^qK^q^^Crx> <Wx^Py> 



+ 



(6.5) 



where the dots represent the A vertex terms which are not considered here and we have 
already replaced the tree-level temporal quark-gluon vertex with its expression Eq. (3.10). 
Also, notice the minus sign in the fourth term, which corresponds to the minus sign arising 
in the Bethe-Salpeter equation considered in Chapter 5. This equation is diagrammatic- 
ally represented in Fig. 6.1. 

The next step is to derive the Dyson-Schwinger equation for the connected 4-point 
quark-antiquark Green's function, which is related to the IPX Green's function via Le- 
gendre transform. Starting with the identity 



<ig/3ig^><iX7iXa>= Spa, (6.6) 

we first take a functional derivative with respect to the source \Js and, with the help of 
the formula Eq. (2.22b), we derive the relation: 

=<iqpiqp><iXpyXa^J5> +\S[k] <U sU ^,\qpiq-y> <\x^yXa> ■ (6-7) 
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Figure 6.2.: Relation between the one particle irreducible (filled blob) and amputated (dashed 
blob) 4-point Greens function for the quark- antiquark system. Internal propagators 
are fully dressed. 

Identifying Js with a gluon source, we have that 

<iXaiXpi«^5>= (-i) <iXQiX7><i^7i9/3i'/'K><iX/3iXp><i«^Ki«^5> • (6.8) 

Taking a further functional derivative of Eq. (6.7) with respect to the quark source vxy^ 
(in this case, Js is identified as quark source), we arrive at the following expression 

<iXaiX5iXAiXr?> = -i <iXaiX5i</p><i9eii?/3i<Ap><iXAiX£><iX/3iX»?> 
+i <iXaiX7> <i^7i9/3i'/'K> <iXAiX<5i-^K> <iX/3iXr,> 
- <iXaiX7><i^7i9Ki9ei9/3><iXAiXe><iXKiX5><iX/3iX»7> • (6.9) 

We now replace the the quark-gluon vertices with the expressions Eq. (6.8) and obtain for 
the 4-point quark-antiquark connected Green's functions, written in terms of IPI Green's 
functions: 

<iXaiX<5iXAiX»7>= 

- [<^Xa^XT><^MiA(t>y><yX^m>] [<iXAiXe><i^£i9/3i</'p><iX/3iXr?>] <iJpiJi.> 
+ [<iXaiX7><i97i9/3i<^K><iX/3iXr,>] [<iXAiXT><igri9Mi'A!^><iX^iX5>] <\J,,\Ju> 
+ [<iXaiX7><iXAiX£><i^7i9Kigeig/3><iXKiX5><iX;3iX7,>] • (6.10) 

The relation Eq. (6.10) can be further simplified by introducing the fully amputated 
Green's function, i.e. dividing by the quark propagators (cut the quark legs), as shown 
in Fig. 6.2. The resulting expression for the IPI Green's function (as function of the 
amputated Green's function), is then replaced in Eq. (6.5). 

6.2. Solution in the heavy mass limit 

In the heavy mass limit and under our truncation, introduced at the beginning of the 
previous chapter, Eqs. (6.5, 6.10) simplify dramatically, allowing us to derive exact solu- 
tions for both the IPI and the fully amputated 4-point quark Green's functions. In this 
section, the consecutive steps required to arrive at the simplified form of these equations 
will be discussed in detail. Furthermore, the physical implications of the corresponding 
solutions will be discussed. 

6.2.1. One particle irreducible Green's function 

Let us start with the IPI Green's function, Eq. (6.5), and the corresponding diagrammatic 
representation. Fig. 6.1, and apply our truncation scheme in the heavy quark limit, at 
leading order in the mass expansion. 
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Figure 6.3.: Crossed ladder diagram that contributes to the IPI 4-point Green's function. The 
upper line denotes the quark propagator, the lower one, the antiquark propagator 
and springs denote the temporal gluon propagator. 



For the quark-antiquark system, we consider the flavor non-singlet Green's function in 
the s-channel, recalling that the quark and the antiquark are regarded as two distinct 
flavors. Hence, the diagrams (a), (c) and (i) are excluded. The diagram (b) (crossed 
ladder type exchange diagram) explicitly reads (see also Fig. 6.3): 



r^glasbl' - ^' ^Wqq5cl>{pi + Uj)T\g„^^{pi + UJ,Pi, -pi - Pi - Uj) 

rggaT/.(P3, W - P3, -w)Wgg;,A(P3 " ^)^qqa\-i{Pi - '^,P2,Pl + P4 + Uj) 

xWZ{-u})W^i{pi+P4 + u}) (6.11) 

Isolating the energy integral over the quark propagators and using the fact that Wqq{k) = 
—W^q{—k) (before the truncation), we find that the energy integral vanishes, due to the 
fact that both quark and antiquark propagators have the same Feynman prescription, just 
as the crossed box contributions from the kernel of the Bethe-Salpeter equation discussed 
in Section 5.4 of Chapter 5: 



J du}Q Wqq{pi + L0)Wqq{u - ps) 



{Pi + P3) + 2m J \ujQ + p^ — m — Ir + is wq — + w, — + ie 



(6.12) 



In the above, the first factor corresponds to the explicit quark propagator, and the second 
factor, to the explicit antiquark propagator (upper and lower line in the diagram Fig. 6.3, 
respectively) . 

Let us now analyze the diagram (d), containing a quark-2 gluon vertex. This vertex 
can be obtained as a solution of the corresponding Slavnov-Taylor identity. The deriv- 
ation is identical to the Slavnov-Taylor identity for the quark-gluon vertex presented in 
Chapter 4. Starting with Eq. (4.7), we functionally differentiate with respect to the quark, 
antiquark and gluon fields. In principle, the resulting equation contains a large number 
of terms, however most of them simplify in the heavy mass limit and under truncation. 
To be specific, after taking the functional derivatives there are four categories of terms 
entering the equation. Firstly, the terms multiplied by a spatial quark-gluon vertex do 
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not contribute, since this vertex is suppressed at leading order in the mass expansion, 
according to Eq. (5.29) (see also the related discussion). Secondly, the terms containing 
a 4-point function Tgqj^u are also of 0{l/m), due to the fact that in the corresponding 
Dyson-Schwinger equation at least one vertex must be at tree-level, and we are at liberty 
to chose this to be TqqA, which is truncated out by the mass expansion. Thirdly, the 
ghost kernels arising from the functional derivatives also vanish, since they only interact 
with the spatial Yang-Mills sector of the theory. Hence, under truncation and in the 
heavy mass limit, only the terms that involve a temporal quark-gluon vertex will survive. 
Explicitly, the equation Eq. (4.7), from which the Slavnov-Taylor identity for the quark-2 
gluon vertex is derived, reduces to: 

= y d^x6{t - xo) 

X <ic7^>) 6{z -x) + <i<> ia'jiz - x) 

<ncx> nnxKz -x)- iT^^ig^^ <nax> ^{z - a;)} . 

(6.13) 

We now functionally differentiate with respect to iqsy,iqpt and icr^ and arrive at the 
following expression: 



= y dxo6{t — xq) 

X |-^^° <^Qptiqeyicr'^^icrl> 6{z - x) - if'"^ <\qptiqeyi(Tl,> 5{z - x)5{x - w) 

+r4 <\qpMaywl,> 5{z - x)5{x -y)- T^^ <iqaz''^qey^crtj> S{z - x)6{x - t)| 

(6.14) 

Replacing the temporal quark-gluon vertex with the expression Eq. (5.23), and using the 
identity Eq. (A. 4) for the generators, it is straightforward to see that the color structure 
cancels and hence Tggo-a is zero (and correspondingly, the diagram (d) vanishes). 

Let us for the moment discard the diagrams (f) and (g), which include the IPI 4-point 
quark Green's function, and the diagram (e), containing a 4 quark-gluon vertex. Then we 
are only left with the diagram (h) and the rainbow-ladder term (j). This simplification 
enables us to derive a solution for the corresponding (truncated) equation for the IPI 
4-point quark Green's function. With this result at hand, we will then return to the 
diagrams (f), (g) and (e), and explicitly show that they cancel (and hence our assumption 
is justified). With these observations, Eq. (6.5) reduces to: 



<iq^iq-.,iq^iq^r,>= - [gT^j^]^^ |<ix/3iXi.><i?i.i?Mi?ri9r,><iX/,iXK><igKig7io-A><ipAip"> 

}■ 

(6.15) 



+ 



It is convenient to express the resulting equation in momentum space. We define the 
momentum space Green's functions via their respective Fourier transform (in order to 
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avoid proliferation of indices, we introduce the convention that xi,ki correspond to the 
index a etc.) 

<iq^iq^iq,iq^>= J dh dk2 dh e-''=i"^-''=^"2-^'=3"3-^'=*"^ri4^,^(fei, ^2, A^g, ^4) 

(6.16) 

and arrive at the following Dyson-Schwinger equation for the IPI 4-point quark Green's 
function in the s-channel (shown diagrammatically in Fig. 6.4): 



du! 



qqarfj, 
Tad 



^tquccsiP^^ -P^ - ^' ^Wggs^ipi + a;)r^g^^^(pi + uj,pi, -pi -pi-uj 



(P3 , P2 - W , Pi + P4 + (^)Wqqf,x (w - P2 )T'lgaXj (w - P2 , P2 , " W ) 



xW^^i-iIj)W^UPi+p^ + u} 



xWf^{-uj)r''^l^^{pi + U,P2- UJ,ps,P4). (6.17) 
In order to proceed, we make the following assumption for the function T^^^: 

r^^\pi,P2,P3,PA) = r(4)(Po;Pl +P4), (6.18) 

with Pq = Pi+P2- This implies that in the above equation the 4-point function ^^^^^j^{pi + 
U),p2—u},p3,p4) does not depend on the integration variable ujq, and hence we can separate 
the energy and three-momentum integrals and perform the energy integration over the 
quark propagators. Moreover, with this ansatz we are allowed to Fourier transform the 
resulting spatial integral back to coordinate space, as shall be explained shortly below. 

After inserting the expressions Eq. (5.34), Eq. (5.43), for the quark and antiquark 
propagators and completing the energy integration, Eq. (6.17) simplifies to: 

[p? + - 2Ir + 2ie] r(,4^^^(Po;Pi + Pi) 



+i 



.9 
1 — 

4 
9' 



Nc- 



dujW^„{ij)W„^{k + i;j) 
j dujW„„{Gj)T^^l^^^{P^-pi + P4 + ^). 



(6.19) 
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Figure 6.5.: Momentum routing for the diagram (g). See text for details. 
Let us now make the following color decomposition for the function T^^^: 

Ta^rr? = ^a-y^m^i + 6anS'yT'^2- (6.20) 

where (for a given flavor structure) F^^^ and T^2^ are Dirac scalar functions. 

At this point, it is convenient to Fourier transform back to coordinate space. In general, 
since Eq. (6.17) might in principle contain momentum-dependent vertex functions, as 
well as mixing of energy and three- momentum variables, this transformation could not be 
carried out. However, in our case momentum-dependent vertices are absent and moreover, 
with the ansatz Eq. (6.18), the energy and tree-momentum integrals have separated such 
that the spatial integral is performed only over two functions (spatial gluon propagator 
and the spatial component of the quark 4-point function). Hence the Fourier transform 
simplifies to the usual convolution product: 

J dCdW^^{u)T^^\q + uj) = j dxe-'"^-'^W„„{x)V^^\-x), (6.21) 

Using the Fierz identity for the generators, Eq. (A. 12), and sorting out the color factors, 
it is straightforward to obtain for the components f[^\ Tj^^: 

Tf\p,,-x) = -^f^y ^^^^ 



2 A^c J Po - 2X, + i^ W^^ (x) + 2\£ 

A^e [(Po - 2Xr){N^, - 2) + ig^CpW^^jx)] 
Po- 2Ir - ig^CpW^^ (x) + 2i£ 

where x is the separation associated with the momentum pi + P4. Inserting the above 
results into the decomposition Eq. (6.20), we find the final formula for the IPI quark 
Green's function: 



F(^) (Pn- 



Po - 2Ir + i^W^^ix) + 2ie 

e . {Po-2Ir)N,{N^-2)+ig^N,CFW,^ix) , 
^''^^^^ Po - 2Ir - ig^CpW^^ix) + 2i6 + '^"'^'^^^ 



where we recall that Pq = Pi + P2 
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m 





n 




m 






n 

















Figure 6.6.: Perturbative expansion of the diagram (e). Boxes comprise m and n gluon legs, 
respectively, with m,n> 1. See text for details. 



Having derived the solution Eq. (6.23) for the IPI Green's function, we return to the 
diagrams (f), (g) and (e) and show that they do not contribute to the final result. To see 
this, we use our usual trick and consider the energy integral. In the case of the diagram 
(g), this reads (with the momentum routing from Fig. 6.5): 



j dujQWqq{pi + w)Wg<y(w -P4)r^^(pi + U ,P2,PZ,UJ - Pa) 



duo -r-Q 



[p5 + ujq + 771 — Zr + ie] [uj^ — p\ + m — Xr + ie 



[p^ + uq + m - Ir + ie] [co^ -p\ + m-Zr + ie] 



[Pq + ujQ + Ci + 2ie] [Pq + Wo + C2 + 2ie] 



(6.24) 



where a, Ci and C2 are energy independent constants that appear in the expression 
Eq. (6.23) for the IPI Green's function. The above formula can be trivially rewritten as 



/ 



dojQ 



[p\ + ojQ + m -Xr + ie] [uj^ - p1 + m - Xr + ie] 
1 a-Ci 



+ 



Pq + u)q + Ci + 2\£ [Pq + Wo + Ci + 2ie] [Pq + wq + C2 + 2ie] 



(6.25) 



Clearly, both terms in the above sum (having the same Feynman prescription) can be 
reduced to differences of integrals over a simple pole, and with the same sign for performing 
the integration in the complex plane and this vanishes (same as for the kernel of the 
Bethe-Salpeter equation and the diagram (b) from above). An identical calculation for 
the integral (f) leads us to the fact that this integral is also vanishing. 

Finally, we are now in the position to show that the diagram (e), containing the 4 quark- 
gluon vertex, is also vanishing. The argumentation is based on our previous findings, 
namely that the diagrams (f) and (g), containing the IPI quark Green's functions, are 
zero. We first observe that the diagram (e) can be written as a combination of diagrams 
of the form shown in Fig. 6.6, where the boxes contain an arbitrary number of gluon legs 
(ladder resummation).^ On the other hand, as a result of the Dyson-Schwinger equation 
for the IPI quark Green's functions, the diagrams (f) and (g) can also be written as a 
ladder resummation, which exactly coincide with the two terms in the diagram Fig. 6.6. 



^To see this, it is enough to analyze the first few terms in perturbation theory of the 4 quark-gluon 
vertex, which are then included into the diagram (e). 
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Figure 6.7.: Dyson-Schwinger equation for the fully amputated quark-antiquark 4-point Green's 
function in the s-channel. 

Hence, the perturbative series of diagram (e) has been reorganized such that although 
the function Tqgqqu itself does not vanish, this 5-point interaction vertex and and the 
gluon line on top of it provide the cancellation at every order perturbatively. In turn, 
this implies that our original assumption is correct and the solution Eq. (6.23) is valid 
nonperturbatively. 



6.2.2. Amputated Green's function 

In the following, we return to the Dyson-Schwinger equation for the fully amputated 
4-point quark-antiquark Green's function in the s-channel and, with the simplifications 
outlined in the previous section, we will derive a solution to this equation. We will verify 
that this is consistent with the IPI Green's function obtained in the previous section and 
moreover, we will analyze the position of the poles and compare them with the Bethe- 
Salpeter equation for physical states. 

The Dyson-Schwinger equation for the fully amputated 4-point quark-antiquark Green's 
function is obtained from the formula Eq. (6.15), by replacing the IPI Green's function 
with the expression Eq. (6.10) and cutting the legs. Thus this equation (shown diagram- 
matically in Fig. 6.7) is equivalent to the Dyson-Schwinger equation Eq. (6.17), derived 
for IPI Green's functions. It is given by: 



G':>.^,^{p+,p.; k.) = Wt{p- k) [ry 

- j dq [Tlq^Wqq{q+)]^^ [t^^^-( -g_ )r^|,] I^^^ (# - 0^^^,^ (g+ , g_ ; A;+ , ) . 

(6.26) 

In the above, the momenta of the quarks are given by p+ = p + ^P, P- = p — {1 — £,)P 
(similarly for k and q), and ^ is the momentum sharing fraction. P indicates the depend- 
ence on the total four momentum, which will become important for the investigation of 
the bound-state contributions to the Green's function. Eq. (6.17) is an inhomogeneous 
integral equation which - unlike the Bethe-Salpeter equation - contains both a resonant 
component (later on used to reproduce the bound state confining energy of the qq system), 
and a nonresonant term. 

The right hand side of equation Eq. (6.17) does not depend on the external energy pQ, 
implying that the 4-point function G*-^^ has to be independent on the relative energy qo, 
and we further assume that G^^^ depends only on the relative momentum p— k. Accord- 
ingly, we replace the quark and antiquark propagators with the expressions Eq. (5.34), 
Eq. (5.43), and as before we perform the energy integration. We arrive at the following 
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expression: 

.(4) 



WT^nT^p- — 



2Zr + 2ie 



(6.27) 



In the above, we have also replaced the vertex functions with their tree-level expressions 
Eq. (5.37), Eq. (5.44). Having integrated out the energy, it is convenient to rewrite the 
above formula back into coordinate space. Using the definition Eq. (6.16) and the relation 
Eq. (6.21), the equation Eq. (6.27) simplifies to: 



Pq - 2Xr + 2ie 



Waa{x)G^^l,,^{P^:S). (6.28) 



Again, we decompose the function G^^^ 



G 



(4) 



(6.29) 



(for a given flavor structure, G^f^ and Gg^^ are Dirac scalar functions), use the Fierz 
identity, Eq. (A. 12), to sort out the color factors, and obtain the following results for the 
components G^^\ ^2^^: 



Gf\P^;x) 



{Pq - 2Zr)W„,{x) 



9_ 

2 y Po - 2X^ + {^W„,{x) + 2i£ 

(Pq - 2X,) 

Po - 2X. - i^ (iVc - ^e) ^-(^) + 



(6.30) 



g(')(Po;x) 



9 



2Nr 



(Po - 2Ir)W^^{x) 



Po-2X, + i^W,,(f) + 2ie 



(6.31) 



Replacing these results in the formula Eq. (6.29), we get the final result for the function 
G(4): 



G(^)(Pc 



{P0-2Ir)W^^ix) 



2 Po - 2X, + i^W^^ix) + 2ie 



(Po - 2X,) 



Po - 2X^ - i3^W^^{x)CF + 2ie 



(6.32) 



A few comments regarding the structure of the above equation are here in order. Firstly, 
a direct calculation shows that our result for the amputated 4-point function is related 
to the result Eq. (6.23) for the IPX Green's function, via the formula Eq. (6.10) (or 
alternatively. Fig. 6.2). 

Also, notice that despite the truncation, the poles in the IPI and amputated Green's 
functions are identical. Moreover, in this approach the physical and nonphysical poles 
disentangle automatically, as opposite to the "standard" QCD setting, where this separ- 
ation does not occur. Using the form Eq. (5.27) for the temporal gluon propagator, we 
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notice that the bound state (infrared confining) energy Pq = a\x\ emerges as the pole of 
the resonant component (first term in the bracket of Eq. (6.32)), for arbitrary number 
of colors. Hence, this provides an explicit analytical dependence of the 4-point Green's 
function on the qq bound state energy, which results from the Bethe-Salpeter equation 
presented in Chapter 5. The term multiplying the bracket has a pole for Nc = and 
hence this cannot represent a physical state. This (common) pole in the nonresonant 
term can be shifted to infinity (as the regularization of Ir is removed) and, as discussed in 
the previous chapter, is nonphysical just like the poles in the quark propagator or in the 
baryon vertex. In the case of the 4-point Green's function, this can be simply absorbed 
in the normalization. Also, the appearance of this spurious pole does not contradict the 
physical results, since the bound state energy, stemming from the first pole of Eq. (6.32), 
is the only relevant quantity. 

6.3. 4-point Green's functions for diquarks 



Let us now consider the diquark 4-point Green's function. This can be easily obtained 
from the equation Eq. (6.5) for quark-antiquark systems, by interchanging the quark 
legs and inserting the appropriate minus signs. We obtain (see also the diagrammatic 
representation from Eq. (6.8)): 



<iqaiq^iq^iq^>= [gT^^la/? j dy5{x-y) 



+ 



+ 



+ 
+ 



<iX/3iXK> <^Jqx'^<li'^(l'n> <^qMu'^<^\,> <iXviXA> <ip|^i/'^J> 

<iX/3iXK> <^K^qMim><^Xx'^X5> <%i9»7io"e> <iPeiPy> 
<iX/3iXK> <^tMMl'^<lri<yx> <^p]SPy> 
<iX/3iX5> <i^5i9r,io-> <wl'^pi> <ig^ig7icjf io-^> <ipii/9y> 
<iX/3iXv> <\qy\qMpi^(lr,> <^Xi^Xk> <iqjqjicr\> <\p\\pp- 



<\q^\qfj\af,> <iXu}XK> <iqJqji(Ti> <ipiipy> 



(6.33) 



As in the qq case, the dots represent the A vertex terms which are not considered here and 
the tree-level temporal quark-gluon vertex has been replaced with its expression Eq. (3.10). 
Also, notice that the diquark is antisymmetric under the exchange of two quark legs. In 
this case we explicitly take into account the flavor structure, i.e. we consider equal (heavy) 
mass quarks but with different flavors. 
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(a) 



(b) 



(c) 



(d) 



(e) 





(f) 



(g) 



(h) 



(i) 



(i) 



Figure 6.8.: Diagrammatic representation of the one particle irreducible 4-point diquark. Green's 
function. Same conventions as in Fig. 6.1 apply. 




p,+p,-(o 



Figure 6.9.: Crossed ladder diagram that contributes to the IPI 4-point Green's function in the 
diquark channel. 



As before, we analyze the diagrammatic representation from Fig. 6.8 and show that the 
same type of cancellations occur. Starting with the diagram (a), we notice that this is a 
crossed ladder type exchange diagram (see also Fig. 6.9): 



■.(O)a 



fj{pi,UJ - pi, -Uj)Wqqfj^{pi - Uj)r'lg^^^{pi -U,ps,U}-pi - ps) 



^qqaT^liP2,P4+^^,Pl + P3 " w)Wgg^,5(-p4 " Uj)T ^g^g^{-p4, - UJ,Pi,Uj) 



xWf,iu))W^',ipi+P3-L3). 



(6.34) 



It has been already shown that the integral over the quark propagators (with the same 
Feynman prescription) vanishes and thus the diagram (a) is zero. A similar type of 
integral arises in the diagram (j), and hence this term is also not giving a contribution. 
Further, the diagrams (b) and (h) are zero, since the corresponding quark-2 gluon vertex 
vanishes according to the Slavnov-Taylor identity, Eq. (6.14). As in the case of the qq 
system, let us for the moment assume that the integrals (f), (i), containing the diquark 
4-point vertex, and (g), containing the 5-point functions ^gqqqa are also zero, and solve 
the Dyson-Schwinger equation with the remaining terms. Having derived the solution, 
we will then return to the diagrams (f),(i) and (g) and show that they vanish. 

Putting all these together, the Dyson-Schwinger equation for the diquark 4-point func- 
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tion with the remaining diagrams, i.e. diagrams (c), (d) and (e) explicitly reads: 

^qqqqcraT-yr] {Pl, P2, P3, Pi) = 



qqcFTV 



n(0)a 
qquafi 



{P2,UJ-P2, -Uj)Wqque{p2 - ^^)^qqae-n{P2 " W, P4, Pi + P3 + ^ 



xWZ{-u)Wt{p^+p^ + u: 



+ I du 



r^gl^lPli - Pl,^Wgqf}5{pi + Uj)r1q^Sr,iPl + ^^P4, - Pl - P4) 



r 



qqaTuiP2,^^ - P2, -l^)Wqqyx{P2 " w)r^g^;^^(p2 " U},P3,Pl + P4 + Uj) 
dc I 



xWfA-^)W^'APi+P4 + c3 



n(0)a 
qqaa 



r 



qqcFTV 



{P2,UJ -P2, -Uj)Wqqux{p2 " 



"qqqqaKX-yrtiPl +UJ,P2- , P3, P4)W^^{(J) 

(6.35) 



After sorting out the color factors and applying the usual separation of the energy and 
three-momentum integrals, we arrive at the following formula: 



P2,P3,Pi) 



1 



-2m- 2Ir + 2ie 



+1 



9 



2N, 

2 N 2 



2Nr 



2Nc5a,jSr^] J dujWaa{i2)Waa{Pl+P3 + 

[{Nl + 1) 2N,5 ] / dLdW„a{i^)Waa{pi+Pi + '^ 



UJ 



9 f f 



X dujQ 



ujQ + Pi — m — Ir + — p% + m + Zr- - ie 

J divW„a{i^)T^xl-fri{Pl+^^,P2-(^,P3,PA)^ 



(6.36) 



Since the diquarks we are not restricted to flavor singlet, we make the following flavor 
decomposition: 

T^(4) ^ r/ r/ -p(c) , r/ r/ -p(p) 

(the superscripts c and p stand for crossed and parallel configurations, respectively). As 
before, in order to carry on the energy integration, and in the view of Fourier transforming 
to coordinate space the resulting spatial integral, we make the following ansatz for the 
components V^'^ and F^^'^: 



f('')(pi +a;,p2 - a;,P3,P4) = F('')(Po;Pi +P3 + '^) 
F(P)(pi +a;,p2 -^^,^3,^4) = F('=)(Po;Pi +P4 + W), 



(6.38) 
(6.39) 



with Pq = Pi + ^2- 
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With these notations, Eq. (6.36) decouples to: 

f 1 

( T,; 4- no 1 = 

.2 



X 



g 



2Nr 



+ [Nc6 



(6.40) 



r(f) (Po;m+P4) = i 



5 



1 



2A^c P° - 2m - 2X,. + 2ie 



2Nr 



- [NcSaxSrt. - SaJrx] j duW„„ {^)^^nX'yr, i^O IPl+PA+l^)^ 

Further, we make the color decomposition: 

p(c,p) _ r r r(c,pl) _i_ X X p(c,p2) 

where r(^'Pi), r(^'P2) 

are Dirac scalar functions. Inserting this into Eqs. (6.40, 
obtain the following set of equations: 

r(-)(Po;Pi +P.) = H)^ po_2..-2r +2i. 



(6.41) 

(6.42) 
6.41), we 



^[N^^+l) j d^W„„{u)W„^{pi+p:i + 



UJ] 



+ / dcjw„„{0j) \N,T'^'^^\Po;pi+P'i + uj) - r(^i)(Po;m +P3 + 



UJ] 



(6.43) 



r(^'i)(Po;m+P4) = (-i): 



2Nc PO - 2?n - 21r + 2ie 
X js'^y duW„a{Cj)W„„{pi +pi + uj) 

+ / dwt^<,,(La) \Nj:'^P^\P^-pi+pi + Lo)-T^P^\PQ-pi+p^ + oj) 



(6.44) 



2iVc -2m- 2Xr + 2ie 
^^9^ j duW„„{Cj)W„„{pi+p:i + uj) 

d uj W„„ {Co) [iV,r(^i) (Po ; m + P3 + ^) - r(^2) ; pi + p3 + , 



(6.45) 



2Afc P^ -2m- 2Zr + 2ie 



^(Af^^l) j dujW„„{^)W„„{pi+pi + Cj) 

dujw„„{Lo) [iv,r(pi)(Po;Pi +P4 + ^) - r(f2)(Po;p', + 



(6.46) 
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In the next step, we transform to coordinate space and rearrange the terms, such that 
we obtain: 



r(^i)(Po;y)-r(^')(Po;y) = (-i) 



T(P'\Po;x)-T(P^\Po;x) = i-i) 



9 



3 \2 



2NJ _ 2m - 2Ir - i^(l - N,)Wa^{y) + 2ie 



(6.47) 



.2 \ 2 



,-r^2 



2NJ _ 2m - 2Ir - i^(l + N,)W„,iy) + 2ie 



(6.48) 



.2 \ 2 



(A^,-l)2H^,,(f)2 



2A^cy Po - 2m - - i^(l - Nc)W„^{x) + 2\e 

(6.49) 



.2 \ 2 



(7Vc + l)'W^.a(i' 



^^2 



2A^cy Po - 2m - 2X, - i^(l + N^)W„„{x) + 2ie 

(6.50) 



In the above, x and y represent the separations corresponding to the momentum pi + p4 
and pi +P3, respectively. The final expression for the diquark 4-point Green's function 
reads: 

r„.^^(Po; y) = 5i^5l^ [ \ {5^.,5rr, + 5^^5r^,) [r(^i) (Po; y) + r^^^) (p^. ^) 



+ I {Sa-ySrr, " .^^r^^r^) [t^^^HPo; x) - r^P^\Po; x) 



:x 



(6.51) 



with the components given by Eqs. (6.47 - 6.50). 

As in the case of the qq systems, with this solution we return to the diagrams (f), (i) 
and (g) . Writing out the explicit form of the energy integrals we notice that their form is 
identical to the quark-antiquark case, since the e prescription is similar, regardless of the 
internal quark (or antiquark) propagator. Thus, these diagrams are also vanishing. 

The above equation can be rewritten such that the pole structure becomes manifest. 
Introducing the notations 



urn = (-i) 

Uy) = (-i) 

we have that 



,2 \ 2 



,-r^2 



{N,-iyW^^{y) 



^NcJ Po - 2m - 2Ir - i^(l - N,)W^^iy) + 2ie 



.2 \ 2 



,-r^2 



{N, + iyW^^{y) 



■ ctT'yri 



{Po;x,y) 



2NJ Po _ 2m - 2X, - i^(l + N,)W^,{y) + 2ie 



{5a^6rn + Sa^6r^) 6^ f+{y) - 6l Jl U{x) 



(6.52) 
(6.53) 



-|- ( Sd'y^T-rj ^ari^T'Y 



5i^5y.{y) + 5i^5l.J^{x)\ } (6.54) 
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Analyzing the pole structure of the above equation we notice that, as in the case of the 
qq systems, we have two different pole conditions: 



From the first equation (corresponding to the color antisymmetric term in the above 
formula), we find that Nc = 2, —1, and from the second equation, we obtain Nc = —2, 1. 
Hence, the only physical solution, with Nc = 2, corresponds to color antisymmetric and 
fiavor symmetric configuration, in agreement with our findings from the previous chapter 
that bound states exist only for color antisymmetric SU{2) baryons (in that case fiavor 
symmetry was implicit). Notice also that the poles in the case of the diquark system 
disconnect, as opposite to the meson case, where they multiply. 




c 



(6.55) 
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Chapter 7. 

Summary and conclusions 



This thesis has been constructed from three major building blocks: the general derivation 
of the quark Dyson-Schwinger equations in Coulomb gauge first order formalism, the 
perturbative studies at one-loop order, and the nonperturbative investigations in the limit 
of the heavy quark mass. We will present our summary and conclusions separately for 
each part. 



7.1. Functional derivation of Dyson-Schwinger equations 

Starting with the QCD Lagrangian, we have introduced the gauge fixing technique and 
the Faddeev-Popov method in Coulomb gauge. Further, we have converted to first or- 
der formalism. There are two reasons motivating this choice: firstly, the (unphysical) 
ghost degrees of freedom can be formally eliminated and secondly, the system can be 
reduced to the "would-be physical" degrees of freedom. Formal here refers to the fact 
that the resulting equations are nonlocal and hence very difficult to approach in practical 
calculations. 

Using functional derivation techniques, in Chapter 2 we have explicitly derived the 
quark field equations of motion. Based on these equations, we have obtained the Feyn- 
man rules and the relevant tree-level proper two-point and vertex functions. We have also 
discussed the general decomposition of the quark (proper and connected) two-point func- 
tions. Based on the Lagrange transformation, we have derived exact relations between 
the corresponding dressing functions. Given that at one- loop order in perturbation theory 
the fourth Dirac structure ^'^7* of the quark gap equation vanishes, the set of equations 
relating these functions is significantly reduced. Starting with the quark field equation of 
motion, we have derived the Dyson-Schwinger equations for the quark proper two-point 
function, quark contributions to the gluon proper two-point functions and quark-gluon 
vertex functions. Moreover, in Chapter 6 we have presented the explicit derivation of the 
four-point Green's function for quark-antiquark and diquark systems. 

Since the Dyson-Schwinger equations build an infinite tower of coupled integral equa- 
tions, they cannot be solved exactly, and hence approximation schemes have to be em- 
ployed. These have to preserve the symmetries of the theory, which are reflected in the 
Slavnov-Taylor identities. Thus, starting with the time- dependent BRST transform we 
have explicitly derived the Slavnov-Taylor identities for the quark-gluon vertices and, in 
the limit of heavy quark mass, for the four-quark-gluon vertex. 



Chapter 7. Summary and conclusions 



7.2. One-loop perturbative results 

In the second part of this thesis, a one- loop perturbative analysis of the quark gap equa- 
tion, quark contribution to the gluonic two-point functions and quark-gluon vertex func- 
tions has been undertaken. The various propagator and two-point dressing functions have 
been evaluated at this order. To this end, the (dimensionally regularized) results for the 
required noncovariant massive integrals have been obtained, using differential equations 
and integration by part techniques. 

The results for the various two-point functions are rather illuminating. Since the sin- 
gularities are absent in both Euclidean and spacelike Minkowski regions, the analytic 
continuation between Euclidian and Minkowski is justified. The second important phys- 
ical results is the renormalization of the quark mass and propagator. Namely, we have 
verified that the one-loop renormalized quark mass agrees explicitly with the calculation 
performed in linear covariant gauges. Also, up to color factors the results for the quark 
propagator dressing function agree with the corresponding results obtained in Quantum 
Electrodynamics. Finally, the correct one-loop coefficient of the perturbative /3-function 
has been obtained. 

Turning to the quark-gluon vertex functions, we have explicitly evaluated their divergent 
parts and considered them in conjunction with the corresponding Slavnov-Taylor identity. 
This identity contains some peculiar objects, the quark-ghost scattering kernels (analogous 
to ghost-gluon kernels known from Yang-Mills theory). In order to gain some more 
information about them, it has been helpful to analyze the Slavnov-Taylor identity at 
one-loop perturbative order. This has been done in two different ways: firstly by using 
the translation invariance of the loop integrals, without explicitly evaluating them, and 
secondly by considering only the divergent parts of the integrals entering the identity. In 
this later investigation, we have employed the results for the divergent parts of the quark- 
gluon vertex functions previously derived and have found that the equation is satisfied. 
Also, the quark-ghost kernels appear not to contain divergences, in agreement with the 
calculations performed using the method of split dimensional regularization. 

7.3. Heavy quarks 

In the third part of the thesis, we have considered the limit of the heavy quark mass. After 
performing a heavy quark transform of the quark field (adapted for Coulomb gauge), we 
have expanded the generating functional of the theory in the mass parameter and retained 
only the leading order. In this case, the system has simplified dramatically: we have 
found that only the temporal gluon propagator contributes at leading order, whereas the 
spatial gluon is suppressed by the mass expansion. Further, we have truncated the Yang- 
Mills sector to include only the non-perturbative (energy-independent) temporal gluon 
propagator and neglect all the higher order Yang-Mills Green's functions. 

In this setting, we have used the full nonperturbative quark equations of QCD in order to 
study the confinement properties. The gap equation, supplemented by the Slavnov-Taylor 
identity, has been solved and we have shown that the rainbow-ladder approximation to the 
quark (and antiquark) propagators is exact in this case. These equations have been used 
along with the Bethe-Salpeter equation for mesons and Faddeev equation for baryons, 
and we have demonstrated that the ladder approximation to the Bethe-Salpeter equation 
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is also exact. 

Using the assumption that the temporal gluon propagator is energy independent, it was 
then straightforward to find solutions for these equations. In the meson case, from the 
Bethe-Salpeter equation we have found that only color singlet mesons and SU{2) baryons 
have finite energy that increases linearly with the distance, i.e. confining solution (and 
otherwise the system is physically not allowed). Further, we have found that there exist a 
direct connection between the temporal gluon propagator and the string tension (at least 
within the approximation scheme considered here). 

Turning to the baryon case, we have considered the Faddeev equation for three-quark 
states in a symmetric configuration. As in the case of the qq systems, we have found that 
the bound state energy between three quarks increases linearly with the separation and 
we have provided a direct connection between the temporal gluon propagator and the 
physical string tension. Furthermore, we predict that the string tension for three-quark 
states is 3/2 times that of the qq system. 

At this point, we make a short comment regarding our approximations. Firstly, we 
have assumed that the temporal gluon propagator is energy independent. As explained 
in the text, this assumption is supported by the lattice results and moreover, from formal 
arguments it can be inferred that this propagator must have at least a part that is en- 
ergy independent, in order to cancel the ghost loops in the Yang-Mills expressions and to 
solve the energy divergence problem of Coulomb gauge. The second approximation was to 
neglect the Yang-Mills vertices and, given that the spatial quark-gluon vertices are sup- 
pressed by the mass, we have found that the temporal vertex remains nonperturbatively 
bare. Since the temporal gluon propagator is dressed, we have the situation of a gluon 
string connecting two naked color charges, and thus the inclusion the Yang-Mills vertices 
would represent the dressing of these naked charges, i.e. the screening mechanism. Thus, 
we anticipate that this would only lower the value of the string tension, and would not 
alter the linear behavior of the confining potential. The fact that the three-gluon vertex is 
irrelevant for the infrared properties of the theory has been also shown in the Hamiltonian 
approach [36]. 

Finally, in the heavy mass limit and within the same truncation scheme, we have also 
considered the 4-point quark Green's functions, and in particular we have studied the role 
of singularities of these functions. Firstly, we have found that both for qq and diquark 
systems the poles (physical and unphysical) naturally separate — as is well-known, this 
would not be the case in the usual QCD calculations. Also, we have found that the bound 
state energies emerge as a pole of the resonant component, thus providing and explicit 
connection between the (nonperturbative) 4-point Green's function and the bound state 
energy resulting from the Bethe-Salpeter equation. 
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Conventions and useful formulae for the 
gauge group SU{Nc 



c 

In this work, natural units are used: 

h = c=l. (A.l) 

Initially we work in Minkowski space, with the following metric: 

5;,. = diag(l,-r). (A.2) 

Where perturbative integrals are to be explicitly evaluated, we analytically continue to 
Euclidean space, i.e. — )• ik^, as explained in the text. 

The group elements of SU{Nc) are unitary Nc x Nc matrices with determinant one and 
can be written in the form 

U{x) = exp{-ir(x)r'^} (A.3) 

where 6°'{x) specify the angle of rotation in color space. 

The A'^^ — 1 (Hermitian) generators of SU (Nc) obey the Lie algebra 

where the numbers f"'^'^ are the completely antisymmetric structure constants of the group. 
Prom the relation 

det U = exp {Tr In U} = exp{-ie"(2;)Tr T"} (A.5) 
it follows that the generators are traceless matrices: 

Tr T" = 0. (A.6) 

In addition, we have the normalization condition 

Tr(T"r'') = ^(^"^ (A.7) 

The sum rules used in the text are: 

T'^T" = Cpt (A.8) 

j,bj,cjdbc ^ It'^Nc (A.9) 

T^T^T" = (cf - ^ J (A. 10) 
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where the Casimir invariant is given by 



iV2 - 1 

Cp = (A.ll) 



2N, 



We note also the Fierz identity: 

2 [TX(S insj = ^o.,6sf, - ^6^p5s^. (A.12) 
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Appendix B. 

Standard massive integrals 



In this Appendix, we present the derivation of some standard Coulomb gauge massive 
integrals used in Chapter 3. Consider the integral: 

/d oj 

In the case /U = = 1 this gives the scalar integral associated with, for example, the 
fermion loop in quantum electrodynamics [3]. 

We present here a method to evaluate such integrals for arbitrary denominator powers 
(developed originally in Ref. [119]) and generalize to the various additional noncovariant 
integrals. We start by writing the Taylor expansion of the massive propagator in terms 
of a hypergeometric function in the following way: 

:iFoU;--^ . (B.2) 



Now, the idea is to use the Mellin-Barnes representation of the hypergeometric function 
ii^o(/^*;^): 



^°^'''^^ = f^2^ / ^^(-^)'r(-^)r(/^ + ^), (B.3) 



where the contour in the complex plane separates the left poles of the F functions from 
the right poles. A first advantage of this representation is that the "mass term" gets 
separated from the massless propagator and the remaining integrals can be calculated 
with the Cauchy residue theorem, as we shall see below. Also, in order to study various 
momentum regimes, the results can be written as a function of either fc^/m^, or m} j}^ . 
This we do by using the standard formulas of analytic continuation of the hypergeometric 
function (for an extended discussion, see [119]). 

Applying Eq. (B.3) to the massive propagator we can rewrite the integral Jmik"^) as: 

ioo 

Jm{k^) = tAiFT^vt// dsdt{m^y+'r{-s)T{-t)r{fi + s)r{u + t) 



(27ri)2r(/.)r(z. 

du} 



(w2)M+^[(A;-a;)2]^+«' 
Inserting the general result for the massless integral, derived in Ref. [30]: 

doj _ [A:2]2-M-^-e r(/i + + e - 2) r(2 - /i - e)r(2 - u - e) 



I 



[a;2]Mp-a;)2]^ (47r)2-^ r(//)r(i/) r(4 - // - i/ - 2e) 



(B.4) 



(B.5) 
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we get for the integral Jm'- 



m 



2\ s+t 



— ioo 

With the change of variable t = 2 — e — fi — v — u — s (for such a replacement, the left 
and right poles of the T function are simply interchanged and therefore the condition of 
separating the poles is not contradicted) we obtain: 

ioo _ 

[m2]2-M-^-£ 1 1 f (w?\ T{-u) 



'^'^^^ ^ (47r)2-^ (2^i) r(^)r(z.) y "^'^ V r(2 - e + n) 

—ioo 

ioo 

X- / dsT{-s)T{2- e - iJL- s)T{-2 + e + V + ^i + u + s)T{^i + u + s). 

(27ri) J 

—ioo 

(B.7) 

To evaluate the integral over s we use the Barnes Lemma: 

ioo 

1 /■ , x^/ , N T(a + c)T{a + d)T{h + c)T{b + d) ,^ , 

—ioo 

and for the integral Eq. (B.7) it follows immediately that 

7 (P) - ["^']'"""''' ^ 1 
'"^ ^ " (47r)2-- (27ri) r{fi)T{i^) 

ioo _ 

/ rfn fH^] r(-n)r(/i + n)r(z. + u)r{^L + u-2 + e + u) 

—ioo 

Closing the integration contour on the right we have: 



(47r)2-s (27ri) r(^)r(z^) 
m2\~^' 1 r(^4 

fc^y j! r(/z + i/ + 2j) 



j=0 

With the help of the duplication formula 



r{2z) = 22--V-i/2r(z)r + ^) , (b.h) 
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we can rewrite Jm. as: 
Jm{k 



2, _ [^2]2-M-v-er(^ + t,_2 + e) 



E 



1 iT{fi + j)T{u + j)T{f, + u-2 + e + j) r(^) r(ii±^- 



2 



.^^ , A;V 22ij! r(//) r(r/) r{fi + i^-2 + e) r{i^ + j) T{ii±^ + j)' 

(B.12) 

The sum is clearly a representation of the hyper geometric 3^2(0, 6, c; d, e; z) (see, for ex- 
ample, Ref. [65]) and we finally obtain: 

"^"^ ^ " (47r)2-- r(/i + z/) 

X3F2 ^^,,u,f, + u-2 + e;^, ^ ; -^J . (B.13) 

A trivial computation shows that the result Eq. (B.13) is consistent with the known results 
in the limit m = 0. All we have to do is to invert the argument of the hyper geometric 
according to the formula (found, for example, in Ref. [120]) 

p. , , ^ r(6i)r(62) 



r(ai)r(a2)r(a3) 



^ [ r(ai)r(a2 - ai)r(a3 - ai) ^ 
\ r(6i - ai)r(62 - ai] 

X3F2 ( ai,ai - 61 + 1, ai - 62 + 1; ai - 02 + 1, ai - 03 + 1; ^ 

V z 

r(a2)r(ai - a2)r(a3 - 02) 
r(fei - a2)r(62 - 02) ^ ' 

X3F2 ( a2, 02 - 61 + 1, 02 - 62 + 1; -ai + 02 + 1, 02 - ^3 + 1; - 

V z 

^ r(a3)r(ai - a3)r(a2 - 03) ^ _^ya^ 
T{bi - as)T{b2 - 03) 



X3i^2 (^a3,a3 - &i + 1,03 - 62 + 1; -ai + 03 + 1,-02 + 03 + 1; -j j- (B.14) 

The same method can be applied to the other noncovariant integrals with different de- 
nominator structures that appear in the text and the general expressions for these read: 

duj [m2]2-M-^-^r(2-^-e)r(/i + i/ + e-2) 



/ 



[u^^]^'[{k - oj)^ + m^]'' (47r)2-^ r(z^)r(2 - e 

du _ [^2]2-M-i^-^ r(| - /i - £)T{fi + iy + e-2) 

[cj2]M[(A:-w)2 + m2]^ ~ (47r)2-£ T{u)r{3/2 - e) 

X2F1I ^i,fi + u + e-2;3/2-e;--^] . (B.16) 
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Appendix B. Standard massive integrals 



Differentiation with respect to or ki gives rise to expressions for integrals with more 
compHcated numerator structure. For completeness, we list here the first order e expansion 
of the integrals arising into the one-loop perturbative expressions considered in this work: 



/ 



(cj^ + m'^)[{k — w)^ + m^] 



/ 



duj 



I 



uj'^[{k — co)"^ + m^] 



aujuji 



duj 

uj'^[{k — w)^ + m^] 



+ V,(B.17) 
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Appendix C. 

Checking the nonstandard integrals 



One way to check analytically the results for the nonstandard integrals introduced in 
Chapter 3, Am and Af^, Eq. (3.65) and Eq. (3.76), respectively, is to make an expansion 
around = and evaluate the resulting integrals with the help of the Schwinger para- 
metrization. Let us first consider the integral A^, Eq. (3.30). Using Schwinger parameters 
[121] (see also [122]), we can rewrite the denominator factors as exponential functions to 
give: 



Ar 



I dad^d'y / doj exp { — (a + /3)tJ4 + 2/3/c4a;4 
+2pk-oj - 131? - Pm^Y 



{a + 13 + 7)w^ 



(C.l) 



Applying similar reasoning as in Ref. [30], we come to the following parametric form of 
the integral (recall that x = k\, y = }?): 



An 



{x + y + vn?)' 



-l-e 



(4^)2-^ 

-1-/3 



d(5 



-T{l + e) 



da 







(a + /3)V2 



a/3 



X 



{a + j3) {x + y + m?] 



+ 



/3(1 - /3)y + /3m2 
X + y + m? 



2-1 



.(C.2) 



For general values of x, the integral above cannot be solved because of the highly nontrivial 
denominator factor. Since there can be no singularities at x = (this would invalidate 
the Wick rotation which, as discussed in Chapter 3, does hold in this case), we can safely 
make an expansion around this point and then integrate. To first order in powers of x we 
have: 



Ar 



>o {x + y + m 



2-e 



(47r) 

+ y(l -/3) 
m? + y 



-V{l + e) 

-2-e 



dl3 
a/3 



da (a + /3)-^/2 



{m? + y){a + l3) 



/3 



m? + y(l - /3) 



{m? + yY 



m2 + y(l-/3) 1"^"" 
w? + y 

{l + e)x +0(x2)| . 

(C.3) 



After performing the integration we get: 

x^o {x + y + m'^)~^~^ 



Amix,y) 



1. 



r(l + e)2Fi { -e,2 + e;l- e; 




y 



+ y 



+ y 



+ Oix^) + 0{e) 



(C.4) 
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In the above formula, we isolated the hyper geometric term in order to evaluate the e 
expansion separately. For this purpose, we differentiate the hypergeometric function with 
respect to the parameters. In general, differentiation of 2-^1(0, b; c; z) with respect to, e.g. 
the parameter 6, gives (similar expressions are obtained for differentiation with respect to 



a, c : 



iF^ ' ' ' ' {a,b;c;z) = 2^ — — - ^'(6)2^1 (a, 6; c; z) , (C.5) 



k=0 



(c)fc 



k\ 



where ^{k) is the digamma function and (a)^ is the so-called Pochhammer symbol 
(defined, for instance, in the standard textbook [120]): 



{a)k 



T{a + k) 



With the help of formula Eq. (C.5), we obtain: 

2F1 {-e,2 + e;l-£;z) = 1 + eln(l - z) + 0{£) 



(C.6) 



(C.7) 



Inserting this back into Eq. (C.4), we can write down the result for the integral Am (to 
first order in powers of x): 



x^o {x + y + m 



2\-l-e 




In 



m 



+ 



+ y J m? + y 



+ 0{x^) + 0{e) 



(C.8) 



which agrees explicitly with the corresponding expansion of the result given in Eq. (3.65). 

We now turn to the integral A^, given by Eq. (3.31). The parametric form has the 
expression: 



(x + ?/ + m^)~^"^r(l + e) 



(4vr) 



2-s 



f-1 rl-fS 

X dp da 
Jo Jo 



/5 



(a + /3)3/2 



a/3 



+ 



/3{l - /3)y + Pm'^ 



(a + /3) {x + y + m?) x + y + m? 



-l-e 



(C.9) 



Calculations similar to the integral Am bring us to the following result (to first order in 
x): 



2x 
'3y 



(47r)^~^ 



m 

y 



2 In- 



m 



ni^ + y 



m' , 
-2 1 + 1 In 



m 



ni^ + y 




+ 0{x') + 0{e) , 



(C.IO) 
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again in agreement with the expansion of the result given in Eq. (3.76). 

Another useful check comes from the study of the mass differential equation, i.e., the 
mass is in this case regarded variable. 

With /" given by Eq. (3.32), the derivative with respect to the mass reads: 

m^= [ I 2 ] . (C.ll) 

From the relations Eq. (3.33), Eq. (3.34) and Eq. (C.ll) we get the following relation: 

Qjn Qjn Qjn 

+ + '"7^ = {d + n- 5)r. (C.12) 

ok4 okk am 

Using the same procedures as described in Chapter 3, we can then derive a differential 
equation for the integral in terms of the mass: 

dm? k"^ + k"^ + m'^ J [{k — uj)"^ + m^J^o;^ 

Starting with the case n = where = A^, we see that by inserting the solution, 
Eq. (3.65), we have that in the limit e — t- 

2dAm m" _ m2(x + y + m2)-2-^ 

^ + U + ^'^) 7,2 , ™2^' 




dm^ k'^ + (47r)2-^ 

+0{e). 

In terms of Schwinger parameters, the explicit integral of Eq. (C.13) reads: 
f duj m? r(l + e) ''^ 



k'^-\-m? J [(A; — 0;)^ + m^J^oj^ x + y + (47r)2 ^ 
and for m? 7^ indeed 



(C.15) 




k'^ + m'^J [(/c - a;)2 + m2]2cj2 (47r)2-s 

+0(e), (C.16) 

showing that the mass differential equation is satisfied. In fact, in the differential equation 
Eq. (C.13), there is a potential ambiguity arising from the ordering of the limits — )• 
and e — )• 0. Namely, for = 0, the right-hand side of Eq. (C.14) vanishes as m'^lnir?, 
whereas the parametric form of the integral in Eq. (C.15) goes like m? je. However, 
this problem is not manifest because of multiplication with the overall factor m? in the 
differential equation. Since the solution of the mass differential equation is in principle 
formally derived as the integral over rr? and m? = is the only the limit of this integral. 



^Recall that when deriving the differential equations presented in Chapter 3, the mass has been treated 
as a parameter, and only the differentiation with respect to and ki has been considered. 
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Appendix C. Checking the nonstandard integrals 



the ambiguity encountered may be regarded as an integrable singularity and presents no 
problem. 

Turning now to the case n = 1 where = Af^, we first extract the overall factor as 
before by defining = k^Am such that the differential equation is 

^2dAn ^ n? -^^ _ rr? f du 

dm? k"^ + m^ ™ fc^ + m^ J [{k — tj)^ + m^J^a;^ 

Notice that in the integral term we have used the identities 



/ 



dio LOi f dLo{ki-LOi) _ f /n 1e^ 



[(fc-u;)2+m2]2^2 J j^2 + ^2]2(^_^)' V [{k - u)^ + m^]^u}^ 



Now, for ^ 0, the integral term of Eq. (C.17) is finite as e — )• 0; however, the = 
limit is again ambiguous but as above this can be regarded as an integrable singularity. 
Also, when = 0, ^ is known to be e finite (it is the massless integral considered in 
Ref. [30]). This means that as e — >• we have the simple integral expression 

= 1 \^ ^ In I \ " I . (C.19) 

dm^ x + y + m^ (47r)2-£ ^ ' ~ ' ' ~ ' 

Knowing the solution, Eq. (3.76), it is straightforward to show that when = the 
original massless integral from Ref. [30] is reproduced and that the derivative of the 
massive solution satisfies the above. 
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Appendix D. 

One-loop non-covariant vertex integrals 



In order to extract the divergence of the non-covariant three-point integrals appearing in 
the quark-gluon vertex functions presented in Chapter 3, we use a method based on the 
Schwinger parametrization [121]. As an example, consider the divergent integral: 



I 



(D.l) 



In principle, one can also consider massive quark propagators, i.e. integrals of the form 
l/w^[(tJ -|- fc3)^][(c(; — /ci)^ — m?]. However, as will shortly become clear, the mass factor 
does not contribute to the divergent part and hence in the following we will set m = 0. 
Using the Schwinger parametrization, the formula Eq. (D.l) can be rewritten as: 



dad(3 d'-f / duo 



1 



2a + /3 + 7 



+ 



(a + /3 + 7)2 



X exp 



|-(a + 7)a;| - (a + /5 + j)uj'^ + 2-/uJiki - 2(2 ■ (/3A?3 - 7^1) - 



(D.2) 



We start by considering the first term (the factor 5jj/2 has been left aside): 

dadpd'y J duj — ^ ^ ^ — exp { — (a -|- 7)0;! — (a + /3 + 7)(i;^ + 27a;4A:4 
-2a;-(/3fc3-7^i)-/3^l'-7fcf} (D.3) 



(in fact, it turn out that the second term in Eq. (D.2) is convergent and therefore not 
interesting for our purpose). After making a shift 



a; — )• w ■ — — — , a;4 — )• a;4 + 



7/E 



14 



a + /3 + 7' ' ' a + 7 
and performing the momentum integration, we are left with the parametric integral 



(D.4) 



1 



(47r 



da df5 d'y 



(a + /3 + 7)5/2-£ (a + 7)1/2 



X exp 



(/34-7fci)2 ^ 72^2 



+ 



14 



a + /3 + 7 a + 7 



f^kl-jkj 



(D.5) 



Following the usual procedure (also used in the Appendix C to check the nonstandard 
two-point integrals), we insert the identity 1 = dX 5{X — a — /3 — 7), make the rescaling 
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a — )• aX, (3 — >• /3A, 7 — )• 7A and perform the integral over A. After rearranging the terms, 
we get: 

Tie) 1 
Jo = \J_ / dadpd-f6{l-a- ^ --/)- 



(47r)2--7o ''(a + 7)^/2 

x|/5(l-/5)fci + 7(1-7)^2 + 2/37^1- fc3 + ^2Mi\ (d.6) 
I a + 7 J 

In the above equation, we see that the integrand is convergent in the hmit e — t- 0, and the 
divergence of Iq is only given by the overall factor r(e) and thus in the exponent we can 
set e = 0. It also becomes clear that the mass term (which would enter the argument of 
the exponential) cannot not influence this result, and hence our initial setting with m = 
is justified. The result is: 

^o = 7T^|- + finite. (D.7) 
(47r)^ 3 e 

By a similar calculation, one can show that the second term in Eq. (D.2) is convergent 
and therefore the overall divergent structure of the original integral is 

r - ^ + finite. (D.8) 



(47r)2 3e 

A straightforward example is the covariant integral 

/duj UJjLOl , 
^77 yTFT y?- CD.9 

This can be calculated without explicitly writing out the parametric form, by using the 
following simple relation: 

doj uji{ijj + k2)^ _ ,2 [ duj uji I' duj ojiuj^ 



0J^{ki-0j)^{k2+Uj)^ '"V 0J^ki-u)^k2+u)^ ''J UJ^ki-Oj)^k2+0Oy 

/doj ujj ,^ 

where the first term and the temporal part of the second term on the right hand side are 
convergent. Using the fact that the result has the form dul* (where /* includes the factor 
1 /d) , it is simple to get the e coefficient by combining the terms on the right hand side of 
Eq. (D.IO). We obtain: 

Finally, let us examine the integral 

J oj'^lki — uj)'^{k2 + ujj^uj^ J J 

As before, we rearrange the terms and obtain 

dio bJiUjjOJi I' duj iOiUJjUJi 



2k2kI*{SijSki + Sii5jk + 6ik6ji 



U^{ki-Ujyuj^ J u'^[k^+ujf(k2-ujf 

UJiUJik2j +UJiUJjk2l+UJiUJjk2l , „ ■w-rv10^ 

dio ^-^ h hnite(D.13) 

a;2(fc3 + w)2(/c2 -a;)2 
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The last integral on the right hand side has been already evaluated and hence we only 
need to consider the integral 



a;2(A;3 +l5)2(A;i - a;)^ 

This we evaluate by differentiating the parametric form of the integral Eq. (D.2) (i.e., the 
integral with the factor WjWj in the numerator) with respect to kn. The result is 

J = J doj {[yidijkii + Sukij + 5jiku) - j3{5ijhii + Suk^j + Sjik^a)] 

^ {-fkij - Pk3j){-fku - Pk3i){-fkii - pk^i) \ ^ 

+ (a + /3 + 7)3 jexpi...| (U.ibj 

(the argument of the exponential function is identical to the one of Eq. (D.2)). After 
making the usual manipulations, we extract the e coefficient from the parametric integral 
and obtain: ^ ^ 

-X^ij^ki + ki^jk + ^ik^ji) + finite. (D.16) 



' (47r)2 20e^ 

We list here all the divergent (covariant and non-covariant) one-loop vertex functions 
used in the evaluation the Coulomb gauge quark-gluon vertices: 

^ + finite (D.17) 





UJ^iki 


-uy{k2+ujf 




doj oj\ 


OJ^ki 


-a;)2(/c2 +a;)2 




duj uJiUJi 


{ki- 


Ujf{k2+LoY0j^ 




doj uj\ 



(47r)2 4e 
1 1 
(47r)2 4e 
5ii 1 



+ finite (D.18) 
+ finite (D.19) 



(47r)2 3e 

,, ,,,, — + finite (D.20) 

{ki-ujY{k2+ioYuj^ (47r)2 e 

/ ^i^l^ji^fc ^ + ^.^ 5 +^.^5 ) + finite (D.21) 

J ijj^{ki — uj)^[k2 + ujYuj''' (47r)2 20e 

j' doj iJlijjjUk 6jk 1 



u^iki-uj)^{k2+uj)^Cd^ (47r)2l2e 



+ finite. (D.22) 
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Appendix E. 

Temporal component of the quark-baryon 
vertex 



In this Appendix we present the exphcit derivation of the energy-dependent part of the 
Faddeev vertex, Eq. (5.84). We start with Eq. (5.83) and consider the first of the per- 
mutations of the energy integrah 

dko r , ■ A — I ^ , ■ 1 ^t{Pi + ko,P2 - ko,pl)- 

[p^ + Kq — m — Ir + ie\ \_p)2 — Kq — m — Ir + \e\ 

(E.l) 

1 1 



Using 



fE 21 

[z -I- a -I- ie] [z -I- 6 -I- ie] {b — a) \ z + a + ie z + b + ie \ 

and shifting the integration variables, we find that the integral /, Eq. (E.l), depends 
only on the momentum (and implicitly on the bound state energy of the system Pq). 
Explicitly, it reads (using the symmetry of F^): 

/ = ^-i- I dk, r.(Po-pO + feo,-.o,p°) 

[Po-P3-2(m-fX^)] y ^ [ko + Po-pl-m-Ir + ie] 
Replacing this in the equation Eq. (5.83), we find: 

rM,plpl) = -2g^CBW„,{\r\) V 



™ [Po-p^-2{m + Xr 



[Po-Pi+kQ-m-Zr + lej 

The form of the equation Eq. (E.4) suggests that the function F^ can be expressed as a 
symmetric sum 

UpUIpI) = M) + fip'2) + fipl), (E.5) 

such that the integral equation for Ft (function of three variables) is reduced to an in- 
tegral equation for the function / (of only one variable). The function /(p^) should be 
chosen such that the integral on the right hand side of the equation Eq. (E.4) generates 
a factor proportional to [Pq — p^ — 2(m H-X^)] , to cancel the corresponding factor in the 
denominator. To examine this possibility, we impose the following condition: 

1 f^j^ /(Pq - p° + fco) + /(-fco) + /(pO) _ a i 

[Po-p°-2(m+X,)] y ° [Po-p° + A;o-m-X, + ie] Pq - 3(m + X,)^ 

(E.6) 
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where a is a (dimensionless) positive constant which remains to be determined. Rearran- 
ging the terms to factorize the function /(/cq), the above equation can be rewritten as 



/ 



dko f{ko 
= (-i) 



+ 



'io - m-Ir + ie Pq- pI -k 
(2a - l)Po - 2ap° + (3 - 4a)(m +X,.) 



m — Xr + ie 

0- 



2[Po -3(m + X,.)] 
Then the most obvious ansatz for the function / is 

1 



f{ko) 



{2a - l)Po - 2ako + (3 - 4a)(m + ie 



(E.7) 



(E.8) 



such that on the right hand side of the equation Eq. (E.7) the numerator is cancelled by 
f{pi)- The next step is to complete the integration on the on the left hand side, which 
gives (note that the e prescription is chosen such that only the first term in the bracket 
survives - the integration must not give rise to any new terms containing the energy p^): 



d/co /(/co) 



1 



kQ — m —Zr + \e 



{2a - l)Po + (3 - 6a){m+Ir) 



(E.9) 



It is then straightforward to compare Eq. (E.7) and Eq. (E.9) and find that the equality 
is satisfied for a = 3/2, leading to the expression for the vertex Ff used in the text. 
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